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Abstract 

Due to coherence, there are strong electromagnetic fields of short duration in very 
peripheral collisions. They give rise to photon-photon and photon-nucleus collisions 
with a high flux up to an invariant mass region hitherto unexplored experimentally. 
After a general survey of the field equivalent photon numbers and photon-photon 
luminosities, especially for relativistic heavy ion collisions, are discussed. Special 
care needs to be taken to include the effects of the strong interaction and nuclear size 
in this case. Photon-photon and photon-hadron physics at various invariant mass 
scales are then discussed. The maximum equivalent photon energy in the lab-system 
(collider frame) are typically of the order of 3 GeV for RHIC and 100 GeV for LHC. 
Diffractive processes are an important background process. Lepton-pair, especially 
electron-positron pair production is copious. Due to the strong fields there will be 
new phenomena, like multiple e^e“ pair production. The experimental techniques 
to select 77 -processes are finally discussed together with important background 
processes. 

Key words: photon-photon processes, photon-hadron processes, relativistic heavy 
ion collisions. 
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1 Introduction and Purpose 


With the first collisions at the relativistic heavy ion collider RHIC at Brook- 
haven (BNL) in Jnne 2000 heavy ion physics has entered a new stage. The 
interaction of nnclei at snch high energies has only been observed np to now in 
cosmic ray interactions. Many interesting physics topics can now be stndied in 
the laboratory. The main aim is to stndy the violent central collisions. One is 
looking for the formation and the signatnre of a new state of hadronic matter, 
the Qnark Glnon Plasma. The present statns of this held is exposed in the 
Proceedings of Qnark Matter QM2001 [1]. 

Relativistic heavy ions are also very important tools for other physics inves¬ 
tigations. We mention here projects (sometimes referred to as “non-QGP” 
physics, a term which we will avoid in the following) like the search for new 
physics at very high rapidities in the CASTOR snbproject at ALICE. This 
project is related to cosmic ray physics and searches for the so called Centanro 
events at the LHC [2]. Other “exotic” physics topics like the possible occnr- 
rence of CP violation or “disoriented chiral condensates” (DCC), see [3,4], 
have also been investigated in the past. 

It is the pnrpose of this Report to review the physics of very peripheral^ 
collisions. We pnt the main emphasis on the energy regions of the colliders 
RHIC (Relativistic Heavy Ion Collider) at Brookhaven and the forthcoming 
LHC (Large Hadron Collider) at CERN/Ceneva. Dne to the coherent action 
of all the protons in the nnclens, the electromagnetic held snrronnding the 
ions is very strong. It acts for a very short time. According to Fermi [5] “fhzs 
time-dependent electromagnetic field can be replaced by the field of radiation 
with a corresponding frequency distributiori\ see also Fig. 1. He called this 
“aqnivalente Strahlnng”. The spectrnm of these photons can be calcnlated 
from the kinematics of the process. The ’’eqnivalent photon method” (EPA) 
is often called also the ”Weizsacker-Williams method”. For a more popnlar 
introdnction see, e.g., [6], and also the corresponding remarks in [7]. 

A very nsefnl view on the electromagnetic processes is the parton pictnre 
with the photons seen as the partons. We give the basic argnment here in the 
introdnction, a detailed explanation will then follow in Chap. 2. In this pictnre 
the scattering is described as an incoherent snperposition of the scattering of 
the varions constitnents. For example, nnclei consist of nncleons which in tnrn 
consist of qnarks and glnons, photons consist of lepton pairs, electrons consist 

^ In the following we will always use “very peripheral” to denote the distant colli¬ 
sions with b > Ri R 2 we are interested in. Here b denotes the impact parameter 
and Ri and R 2 are the radii of the two nuclei. The term “ultraperipheral” collisions 
is also sometimes used for them. This is an attempt to distinguish it clearly from 
what is sometimes called “peripheral collisions”, where b rs Ri R 2 . 
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of photons, etc.. Relativistic nuclei have photons as an important constituent, 
especially for low virtuality and a low ratio of the photon energy compared 
to the one of the ion (see below for quantitative arguments). This is due to 
the coherent action of all the charges in the nucleus: for these conditions the 
wavelength of the photon is larger than the size of the nucleus, therefore it 
does not resolve the individual nucleons but sees the coherent action of them. 
This has some similarity to the “low-x physics” and the Weizsacker-Williams 
approach of [ 8 ] to gluons in the initial state of heavy ion collisions. 

The coherence condition limits the virtuality = —q^ of the photon to very 
low values 


< l/R^ (1) 

where the radius of a nucleus is approximately R = 1.2 fm with A the 
nucleon number. This is due to the rapid decrease of the nuclear electromag¬ 
netic form factor for high values. For most purposes these photons can 
therefore be considered as real (“quasireal”). From the kinematics of the pro¬ 
cess one has a photon four-momentum of = (uj, g±, q^ = u/v), where u and 
qj_ are energy and transverse momentum of the quasireal photon in a given 
frame, where the projectile moves with velocity v. This leads to an invariant 
four-momentum transfer of 


Q" 


UJ n 

— + 

7 


( 2 ) 


where the Lorentz factor is 7 = E/m = 1/Vl — The condition Eq. (1) 
limits the maximum energy of the quasireal photon to 


^ 7 UJmax ~ 


2 

R’ 


( 3 ) 


and the perpendicular component of its momentum to 




1 


( 4 ) 


At LHC energies this means a maximum photon energy of about 100 GeV in 
the laboratory system, at RHIC this number is about 3 GeV. We dehne the 
ratio X = oj/E, where E denotes the energy of the nucleus E = Mjq'-^A and 
Mtv is the nucleon mass. It is therefore smaller than 


_ 1 _ 

“ RMnA ~ R ’ 


( 5 ) 
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where Ac(^) is the Compton wavelength of the ion. Xmax is 4 x 10“^, 3 x 10“^, 
1.4 X 10“^ for O, Sn, Pb ions, respectively. Here and also throughout the rest 
of the paper we use natural units, i.e., h = c = 1. 



Fig. 1. A fast moving nucleus with charge Ze is surrounded by a strong electro¬ 
magnetic field. This can be viewed as a cloud of virtual photons. These photons can 
often be considered as real. They are called ^^equivalenf or ^^quasireal photons”. In 
the collision of two ions these quasireal photons can collide with each other and with 
the other nucleus. For very peripheral collisions with impact parameters b > 2R, 
this is useful for photon-photon as well as photon-nucleus collisions. 


The collisions of and e~ has been the traditional way to study 77 collisions. 
Similarly photon-photon collisions can also be observed in hadron-hadron col¬ 
lisions, see Fig. 1. Since the photon number scales with [Z being the charge 
number of the nucleus) such effects can be particularly large. This factor of 
> 6000 (corresponding to Au, Z = 79) is the reason why the name “gold 
flashlight” [9] has been used to describe very peripheral (AuAu) collisions at 
RHIC. 

Similarly the strong electromagnetic held can be used as a source of photons 
to induce electromagnetic reactions in the second ion, see Fig. 1. Since the 
ion, which is hit by these photons, is moving in the collider frame the photon- 
hadron invariant masses can become very high. In the rest frame of one of the 
ions (sometimes called the ’’target frame”) the Lorentz factor of the other ion 
is given by 7 * 0 ^ = ‘^ifab ~ 1) where ■jiab is the Lorentz factor in the collider (cm) 
frame. The maximum photon energy in this frame is 500 TeV for the LHC 
and 600 GeV for RHIC. 

This high equivalent photon flux has already found many useful applications in 
nuclear physics [10], nuclear astrophysics [11,12], particle physics [13] (some¬ 
times called the “Primakoff effect”), as well as, atomic physics [14]. Previous 
reviews of the present topic can be found in [15-18]. 

The theoretical tool to analyze very peripheral collisions is the equivalent 
photon method. This method is described in Chap. 2. The equivalent photon 
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method is often used for high energy reactions, e.g., in the description of 
e’''e“ and ep collisions. We put most emphasis here on the peculiarities which 
arise due to the hnite size and the strong interactions of the ions. These 
equivalent photons can collide with the other nucleus and with each other. This 
leads to the possibility of doing photon-hadron and photon-photon physics at 
relativistic heavy ion accelerators, with high fluxes of equivalent photons in 
hitherto inaccessible invariant mass regions. 

Thus there will be new possibilities to study photon-hadron interactions at 
RHIC and LHC. It extends the yp interaction studies at HERA/DESY to yA 
interactions. The photon-hadron invariant mass range at RHIC will be some¬ 
what below the one at HERA, whereas at LHC one will reach higher invariant 
masses than those possible at HERA. As was mentioned above, these photons 
can be regarded as quasireal and the freedom to vary the four-momentum 
transfer is not given in the heavy ion case. According to [19,20] relativis¬ 
tic heavy ion colliders can be regarded as “vector meson factories”. This is 
discussed in detail in chapter 5. 

This vector meson production is essentially a diffraction process: the (equiv¬ 
alent) photon emitted by one of the nuclei is diffractively excited to a vector 
meson on the other nucleus. This can be viewed as a photon-Pomeron in¬ 
teraction. In addition there are also Pomeron-Pomeron processes, which are 
interesting processes of their own [21]. Their importance as a possible back¬ 
ground will be studied in chapter 6. 

The physics potential of photon-photon collisions is discussed in chapters 3 
and 4. These chapters extend our previous studies for the CMS Heavy Ion 
Programme [22,23], see also [24]. It ranges from studies in QCD (Chap. 3) like 
meson production and the total yy —hadron cross section to the search for 
new particles (Chap. 4) like a light Higgs particle. 

Chapter 7 is devoted to lepton pair production. Due to their low mass es¬ 
pecially the electrons have a special status and EPA is not always a good 
approximation. Effects of the strong fields manifest themselves essentially in 
multiple pair production. Besides the free pair production we also discuss 
some important related processes like bremsstrahlung from these leptons, the 
production of muons with a large transverse momenta and bound-free pair 
production. 

Central collision events are characterized by a very high multiplicity. Therefore 
all major heavy ion detectors are tuned to deal with the large amount of data 
in this case. On the other hand, the multiplicity in very peripheral collisions 
is comparatively low (therefore the name “silent collisions” has been coined 
for them). The ions do not interact strongly with each other and move on 
essentially undisturbed in the beam direction. The only possible interactions 
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are due to the long range electromagnetic interaction and diffractive processes. 
Still some of these ’’silent events” are not so silent, e.g., quite a few particles 
will be produced in a 100 GeV on 100 GeV 77 -interaction leading to a hnal 
hadron state. The background coming from especially grazing collisions needs 
to be taken into account. Also a way to trigger on very peripheral collisions is 
needed. This will be discussed in chap. 8 . 


1.1 A Short History of Very Peripheral Collisions 


Relativistic heavy ion collisions have been suggested as a general tool for 
two-photon physics about twelve years ago and the field has grown rapidly 
since. Yet the study of a special case, the production of pairs in nucleus- 
nucleus collisions, goes back to Landau and Lifschitz in 1934 [25] and to Racah 
in 1937 [26]. In those days one thought more about high energy cosmic ray 
nuclei than about relativistic heavy ion colliders. The general possibilities and 
characteristic features of two-photon physics in relativistic heavy ion collisions 
have been discussed in [27]. The possibility to produce a Higgs boson via 77 
fusion was suggested in [28,29]. In these papers the effect of strong absorption 
in heavy ion collisions was not taken into account. Absorption is a feature, 
which is quite different from the two-photon physics at colliders. The 

problem of taking strong interactions into account was solved by using impact 
parameter space methods in [30-32]. Thus the calculation of 77 luminosities 
in heavy ion collisions was put on a hrm basis and rather definite conclusions 
were reached by many groups working in the held, as described, e.g., in [16,15]. 
This opens the way for many interesting applications. 

Up to now hadron-hadron collisions have rarely been used for two-photon 
physics. The work of Vannucci et ah [33] in 1980 may be regarded as its 
beginning, see Fig. 2. In this work, the production of pHpi~ pairs was studied 
at the ISR. The special class of events was selected, where no hadrons are seen 
associated with the muon pair in a large solid angle vertex detector. In this way 
one makes sure that the hadrons do not interact strongly with each other, i.e., 
one is dealing with very peripheral collisions; the photon-photon collisions 
manifest themselves as “silent events”, that is, with only a relatively small 
multiplicity. It seems interesting to quote from Vannucci [33]: “The topology 
of these events, their low transverse momentum, and the magnitude of the 
eross section can he most naturally be interpreted by the 77 process. This 
effect, which is still a small background at the ISR compared with the Drell-Yan 
mechanism, grows with energy and could trigger enough interest to be studied 
for itself, possibly at the pp collider, probably at the ISABELLE machine. ” 
Although ISABELLE has never been built, we have now the Relativistic Heavy 
Ion Gollider (RHIG) in this ring, where 77 processes will be studied with the 
“gold flashlight” in this way. Another example is the search for a magnetic 
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b) 


(a) (b) 

Fig. 2. The production of a muon pair at the ISR with (a) and without the produc¬ 
tion of hadrons (b) is shown. Typical characteristics of very peripheral collisions are 
seen: low particle multiplicity and a small sum of transverse momenta. Reproduced 
from Fig. 2, p. 244 of [33] with kind permission of Springer-Verlag and the author. 

monopole in photon-photon elastic scattering in pp collisions at the Tevatron 
following the idea in [34]. This is discussed in chapter 4. Dimuons with a very 
low sum of transverse momenta (i.e., coming from 77 -interactions) are also 
considered as a luminosity monitor for the ATLAS detector at LHC [35]. A 
more recent discussion can also be found in [36]. 

The physics potential of relativistic heavy ion collisions for photon-hadron 
studies was gradually realized during the last decade. Nuclear collisions with¬ 
out nuclear contact have been used very successfully over the past 50 years to 
study nuclear structure, especially low lying collective excitations. Coulomb 
excitation, mainly at energies below the Coulomb barrier has been reviewed by 
Alder et ah [37-39]. The theory of relativistic electromagnetic excitation was 
given by Winther and Alder in [40]. This seminal paper started a whole new 
era of investigations. It was recognized for some time that the cross section 
for the excitation of the giant dipole resonance is huge. This is the reason why 
also new phenomena like the double giant dipole resonance excitation could 
be studied by electromagnetic excitation. The excitation of the giant dipole 
resonance is also an important loss process in the relativistic colliders, as the 
giant dipole resonance decays predominantly by neutron emission. The decay 
neutrons are also a useful tool to measure the collider luminosity. 

In addition it was recognized in [41] that there is a sizeable cross section for 
photon-nucleus interactions beyond 1 GeV. Thus a held of studies similar to 
the one at HERA will be opened. However, as mentioned already above, in 
the heavy ion case the photon is restricted to be quasireal, and the study of 
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interactions of virtual photons with high with hadrons is not possible. On 
the other hand, the photon-nucleus interaction (rather than the photon-proton 
interaction) can be studied with RHIC and LHC. This is described in [19,20]. 
Since the vector meson production cross section is very big — of the order 
of the geometric cross section — relativistic heavy ion colliders may justly be 
called “vector meson factories” [19]. Photon-gluon processes were discussed in 
[42] and later in [43,44]. The very interesting possibility to produce ti in this 
way was recently described in [45]. 

At the STAR (Solenoidal Tracker At RHIC) detector at RHIC a program 
to study photon-photon and photon-nucleus (especially diffractive “photon- 
Pomeron”) interactions in very peripheral collisions exists [9,46-49]. First ex¬ 
perimental results are just becoming available [50,51]. At RHIC the equivalent 
photon spectrum goes up to about 3 GeV. Therefore the available 77 -invariant 
mass range is up to about the mass of the rj^. At the RHIC/INT workshop at 
the LBNL (Berkeley), the physics of very peripheral collisions was discussed 
by Klein [52] and Brodsky [53]. The experimental feasibility of making these 
measurements with STAR were described, expected backgrounds along with 
the techniques and triggering algorithms to reject these signals are discussed, 
see also [9,46-49]. 

When the “Large Hadron Collider” will be scheduled to begin taking data 
in 2006/2007, the study of these reactions can be extended to both higher 
luminosities but also to much higher invariant masses, hithero unexplored, 
see, e.g., [54,55,22]. 

To conclude this introduction we quote Bjorken [56]: It is an important portion 
(of the FELIX program at LHC [24]) to tag on Weizsdcker Williams photons 
(via the non-observation of completely undissociated forward ions) in ion-ion 
running, creating a high luminosity 77 collider. Although the FELIX detector 
will not be realized in this form at the LHC, these photon-photon and photon- 
ion collisions can and will be studied at other LHC detectors (ALICE and 
CMS). 
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2 Prom Impact-Parameter Dependent Equivalent Photon Spectra 
to 77 -Luminosities 


Electron-positron colliders have been and still are the basic tool to study 77 
physics, see, e.g., the rich physics program at LEP [57,58]. The lowest order 
graph of two-photon (TP) processes is given in Fig. 3(a). The plane-wave 
description is very accurate for e~^ and e~, which only interact electromagnet- 
ically. This is, for example, described in great detail in [59]. For lepton-hadron 
and lepton-ion collisions (e.g., for ep, or eA) and hadron-hadron, hadron-ion, 
and ion-ion collisions (like pp, pA, and AA) such kind of two-photon interac¬ 
tions occur as well. If both particles are nuclei with charge Z (symmetric AA 
collisions), the cross section of a TP process is enhanced by a factor of Z^. 
For eA and pA collisions the enhancement is a factor Z^. This has been called 
“the power of coherence” by Brodsky [53]. In these later cases there is also the 
interesting possibility of photon-ion processes, see Fig. 3(b), which we will call 
one-photon (OP) processes in the following. 

Under rather general circumstances these processes can be described by the 
concept of equivalent photons. In this way a great simplihcation is achieved. 
The cross section factorizes into an equivalent photon spectrum n{u) and the 
cross section for the photon-ion interaction process in the case of one- 

photon processes: 


O'OP 


duj 

u 




( 6 ) 


For the general case of two-photon collisions the cross section for the reaction 
A 1 A 2 A[A 2 ^f factorizes into the photon-photon luminosity and the 
cross section of the photon-photon interaction process '-py ^ Xf where Ai, A 2 
are the initial particles. A), A '2 their hnal state after the photon emission, and 
Xf the hnal state produced in the photon-photon collision, see also Fig. 4 and 
Eq. (14) below: 


O'TP — 


/—ni(ct;i) f ^^n2{u}2)o-^^^XfW4:u;iU}2) 

J Ui J OJ 2 


doj 


iOl 

dW dLry^ 
dW dL 


dWdY 




77 


W dW 




(7) 

( 8 ) 

(9) 


The cross sections for the corresponding subprocesses are given by and 
the equivalent photon spectra and the 77 luminosity are denoted by 
n{uj) and dL^^/dWdY, respectively, where Y is the rapidity of the produced 
system Xf. The invariant mass of the 77 system is denoted by IF = 
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The 77 luminosity is obtained by a folding of the equivalent photon spectra, 
as will be explained below. 

As compared to the point-like and non-strongly interacting electrons and po¬ 
sitrons there are in the case of incoming hadrons or nuclei additional initial 
and hnal state interaction effects. They must be taken into account in the cal¬ 
culation of the equivalent photon spectra and 77 luminosities. We distinguish 
three conceptually different kinds: 

(1) For extended objects one has an extended charge distribution, as well as, 
an excitation spectrum. This is taken care of by elastic and inelastic form 
factors, see Fig. 3(c). 

( 2 ) For strongly interacting particles (like in pp, pA, or AA systems) there 
are effects like absorption (nuclear interaction), which modify the 77 lu¬ 
minosities. They may be called “initial state interactions”, see Fig. 3(d) 
and (e). In the case of heavy ions the “initial state interaction” con¬ 
sists of both the nuclear absorption and the Coulomb interaction. For 
the Coulomb interaction the Sommerfeld parameter 7 ^ Z 1 Z 2 Q; is larger 
than one and the semiclassical description can be used. The main ini¬ 
tial state interaction is then the nuclear absorption for trajectories with 
b < 2R. This can be very conveniently taken into account by introduc¬ 
ing the appropriate cutoff for the impact parameter dependent equiva¬ 
lent photon spectrum. It would be much more cumbersome to take these 
strong interaction effects into account in a purely quantal description. 

(3) Particles produced in 77 interactions, which have strong interactions, can 
have hnal state interactions with the incident particles, see Figs. 3(f). 


2.1 Equivalent Photon Numbers 


Before looking in more detail at the equivalent photon approximation for the 
different cases of interest, let us discuss here hrst some of the main characteris¬ 
tics in a qualitative way. Here we make use of the plane wave approximation, 
neglecting the “initial state interaction” between the two ions. One should 
remember that therefore the statements made are only rather qualitative con¬ 
cerning the heavy ion cases, where the semiclassical straight line approxima¬ 
tion is the relevant one, as the Sommerfeld parameter is always larger than 
one, see below and should be used instead. 

For the elastic emission from a nucleus with spin 0, characterized by an elastic 
form factor Fel (g^) and using the plane wave approximation, the equivalent 
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Fig. 3. In the collision of either leptons with hadrons or hadrons with hadrons 
photon-photon (a) and photon-hadron collisions (b) can be studied. The principal 
diagrams are shown schematically here. In the collisions of hadrons additional effects 
need to be taken into account: inelastic photon emission processes (c), “initial state 
interaction” (d) and (e), as well as final state interaction (f). 

photon number of a nucleus with charge Z is given, see Eq. (24) below, by 



The elastic form factor F^i is characterized by the property, that it is ~ 1 up 
to values of = —q^ of about 1/R^ and quickly falls of for larger values 
of Q^. This leads to the two main characteristics already discussed in the 
introduction: for on > 'j/R the equivalent photon number quickly decreases 
due to the fall-off of the form factor, giving a maximum “usable” photon 
energy of Umax - l/R- 

In addition also the transverse momentum distribution of the photon is limited 
to rather small values. This is again due to the fall-off of the elastic form 
factor, restricting the transverse momenta essentially to g_L ^ 1/R. But in 
addition the integrand of the expression Eq. (10) peaks at values of gj_ ~ u/j. 
Most equivalent photons will therefore have even smaller transverse momenta. 
This will also apply to the TP case. Also here the produced system will be 
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characterized by rather small transverse momenta with Pj_ ^ l/P. 


Approximating rather crudely the form factor by 0(1//?^ — where = 
—q^, the integration can be done analytically and the equivalent photon spec¬ 
trum in the leading logarithmic approximation is given by 


n{u) 


TT \UiRJ 


( 11 ) 


There are more rehned formulae for the equivalent photon approximation than 
Eq. (11). 

For heavy ions the relevant equivalent photon number is the one obtained in 
the impact parameter approach given below, see Chap. 2.6. The approximate 
form of Eq. (11) can be obtained also in this case from Eqs. (44) and (46) 
below. This equation is useful for quick estimates and also to see in a simple 
way the general properties of the spectrum. We see that the spectrum falls off 
logarithmically with oj up to Umax = 1 /R- Note also that the l/u dependence, 
characteristic of the equivalent photon spectra, has been taken out in the 
dehnition of n{(jj) in Eq. (9). 

Equivalent photon spectra of point-like particles like the electron require a 
further discussion. The integral Eq. (10) is divergent, if the form factor is 
constant. The choice of the cutoff depends then on the specific process to be 
considered. This is thoroughly discussed in [59], and it is sufficient here to 
refer to this reference. 


The photon-photon luminosity can be found (again neglecting possible initial 
state effects)from Eq. (9). One obtains, see also Eq. (39) and Eq. (49) below. 


dWdY W V 2 J 


n 



( 12 ) 


Due to the decrease of the equivalent photon number with u, this luminosity 
has a maximum at rapidity Y = 0. It falls as a function of |E|. The width 
Ymax of the distribution can be estimated to be 

Ymax ln(2a;max/IE). (13) 


The width decreases with larger invariant mass W. 
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2.2 The Quantum Mechanical Plane Wave Formalism 


Let us first deal with the effects due to the hnite size of the particles. This 
leads to two effects: we need to introduce a form factor for the elastic pho¬ 
ton emission and in addition have to take the inelastic photon emission into 
account, where the particle makes a transition to an excited state. Both can 
be dealt within the general plane wave framework, we are showing here. We 
will in the following deal only with the derivation of the eqnivalent photon 
approximation in the two-photon case; the one-photon case then follows along 
the same line. 

Let us look at the most general case, where two incoming particle Ai and A 2 
nndergo a transition to A'^ and A '2 while emitting a photon each and where 
the collision of these photons then prodnces a hnal state Xf, see Fig. 4. We 
follow in onr derivation here [59], bnt are more general in treating both the 
elastic and inelastic case and nsing a somewhat different notation, see also 
[60]. The cross section of this two-photon process can be written as 


(27r)‘^(47rQ;)^ /■ c-m/ 1 \ imi2 

d(TA,A2^A[A'^Xf = -- J (?1 +d2- Pf) \T\ 

d'^qid^q2dRi{pi - gi)di?2(p2 - q2)dRf{qi + ^ 2 ), (14) 

where pi,p 2 are the fonr-momenta of the initial particles, Ei = p^, E 2 = 
P 2 their energies; qi,q 2 are the fonr-momenta of the exchanged photons and 
dRi{p'j) denotes the integration over the hnal state i (the snmmation and 
averaging over the hnal and initial spin states is assnmed to be implicitly 
included): 


dRM) = n 

k£i 


J(4) 

{2xf2Ek 


Pi 



(15) 


The square of the matrix element |Tp is given by 


\T? 


1 1 


■p/i*-p ly-p I/'* 
-'- 1^-1 2 -*- 2 


M, 


flfl' UV 1 


(16) 


where describes the prodnction of the system Wj in the photon-photon 

interaction and Tj are the electromagnetic transition cnrrents. Integrating over 
all possible hnal state with momentnm p'- = Pi — qi {i = 1, 2), we can reexpress 
these electromagnetic currents in terms of the electromagnetic tensor , 
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see, e.g., [61] 


»r' = 5^ / rffl.(p- ?)(2ir)‘rSTf. 


(17) 


Ai'.Pi’ 


Xf, Pf 


A2’,P2’ 

Fig. 4. The general form of a two-photon process is shown. The two incoming par¬ 
ticles Ai and A 2 either stay in their ground states or undergo a transition to exited 
states Ai and while each emitting a photon. The two photon fuse to a final 
state Xf. Also shown are the momenta of all particles involved 

From Lorentz invariance and current conservation considerations the general 
form of the tensor can be written as 







Pi -Qi A 
qf 


(18) 


where hFj^i and Wi ^2 are two scalar functions of the two independent invariants 
qf and z/j = —pi ■ qi/rrii, which characterize generally the electromagnetic 
structure of the electromagnetic currents. Please keep in mind that we are 
treating the photons as being emitted here, whereas in electron scattering, 
the photons are normally assumed to be absorbed by the particle. The cases 
differ only by the sign of q. 


At high energies the main contribution to this cross sections comes from small 
values of q^ and where the three-momentum of the photon is almost aligned to 
the beam axis. In this case only the transverse part (with respect to the mo¬ 
mentum of the photon) of the tensor is important. Assuming in addition 
that in this case essentially only depends on the energies of the two 

photons {uji = go,i) and not on the photon virtualities qj, it can be related to 
the cross section of the corresponding process for two real photons. For 
a more thorough discussion of the applicability of the EPA, we refer to [59]. 
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Under these assumptions the cross section can be written as 


da 


'AiA2^A\A'^Xf 


/— f ^^ni{u!i)n 2 {ui 2 )da^^^Xf{i^i,i^ 2 ), (19) 

J UJ\ J UJ2. 


doj' 

UJ2 


with the equivalent photon number given by 



a ufrui 


a 


= ^ dqi± dvi 


TT^ 


„2 p2 

21 Uu + 


LOtmt 


{q! 


2^2 


P? 


w^,l + qLw^,2 


( 20 ) 

( 21 ) 


where ^ Pi are energy and momentum of the incoming particle. This 
equation corresponds to Eq. (D.4) of [59]. 


One important feature of the description of the electromagnetic structure (and 
therefore of the equivalent photon spectrum) in terms of fUi and IU 2 is the 
fact that they are Lorentz scalars and can therefore be calculated in any frame 
— e.g., in the rest frame of the projectile (nucleus). 


2.3 Elastic Photon Emission 


Let us hrst discuss the most relevant case: the elastic photon emission. In this 
case the two invariants are not independent of each other but are related via 
2mu = —q^. Therefore the two structure function will be of the form 

1Ui,2(z2, q^) = Wi^2(q^)5(v + q^/2mi) (22) 


and the integration over v can be done trivially. In many cases — and espe¬ 
cially in the heavy ion case — we have = —q^ mU we can then neglect 
the recoil and use to a good approximation z/ ^ 0. We then get qz = oj/p., 
where f3 is the velocity of the particle and 


For the case of a nucleus with a = 0+ ground state (an even-even nu¬ 
cleus) we have Wi = 0 and IU 2 = Z‘^ |Fe;(—g^)|^, with the elastic form factor 
normalized to Fel (0) = 1. In this case the equivalent photon number is given 
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by 


( ^ /■j2 r^2/ 2^ 

= —J‘‘ )■ 


( 24 ) 


It should be clear that in this equation coming from the plane wave approxi¬ 
mation the important strong interaction effects (which correspond to a cutoff 
in impact parameter space, as will be explained below) are not taken into 
account. For a collection of some other relevant cases see, e.g.. Table 8 of [59] 
(with the relation of C and D there to Wi and W 2 given by C = p^hFi, 

D = 14 ^ 2 ). In many cases the dominant contribution comes from IF 2 alone and 
Wi can be neglected. 

For nuclei various degrees of sophistication for the elastic form factor can be 
applied: Gaussian form factors or form factors corresponding to a homoge¬ 
neously charged sphere have been used in the literature. A comparison of the 
different results can be found, e.g., in [16]. 

One can include also the effects of electromagnetic moments like Ml and 
E2, with their corresponding contributions to Wi and IF 2 , see also Eqs. (31) 
and (31) below. E.g., the projectile ^®^Au — which is used at RHIC — has a 
magnetic moment of p = 0.14485 n.m. and an electric quadrupole moment of 
q = -1-0.6 barn [62]. A classical version of an E2 equivalent photon spectrum 
is given in [63]. Perhaps at some level of accuracy of the RHIC experiments 
the inclusion of such effects will become necessary. 

Quite general spectra of the equivalent photons in the plane wave approxi¬ 
mation are given explicitly in [16] (with their dehnition of n{u!) including an 
extra l/cu): 



uj/'r 


for any elastic form factor E^i- The elastic form factor for extended charge 
distributions is in general characterized by a radius R. Dehning a maximum 
photon energy cuq = 7 /i?, the equivalent photon number can then be cast into 
the form 


n{(jj) 


2Z‘^a 

n 




(26) 
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where / is a universal function, independent of 7 and given by 



X 



( 27 ) 


Three different cases are given in [16] point particle, homogeneous charged 
sphere and a Gaussian form factor. The exact result for the equivalent photon 
spectrum of the proton can be found in [64], where the usual dipole parame¬ 
terization of the form factors is used. 


2.4 Inelastic Form Factors 


The structure of the electromagnetic current in Eq. (18) is equally valid for 
inelastic photon emission processes. In this case Wi{u,q‘^) and W 2 {n,q‘^) are 
functions of both invariants u and q^. At low excitation energies, discrete 
states will dominate — resonances of either the nucleus or the proton. These 
resonances will be at well dehned discrete energies (neglecting their hnite 
width), and — as we can again neglect recoil effects in most cases — u has 
to be equal to A, where A is the energy of the excited state in its rest frame. 
This can be seen by evaluating the relation {p — qY = (m-|-A)^ and neglecting 
terms quadratic in q and A. Therefore, similar to the elastic case, Wi and W 2 
will be of the form 


= n\2{q^)6{v - A), 


( 28 ) 


and the integration over du in Eq. (21) can be done trivially. The effect of the 
hnite excitation energies on can be taken into account [60], with in this 
case given by 



~ ^- 

7 7 


A' 2 

H-y + Q'u- 

7 


(29) 


One clearly sees that as long as cu > 7 A the modihcation of g^ is small and 
can be neglected in these cases. 

In the case of nuclei it is more convenient to express Wi and W 2 in terms 
of the more familiar (to nuclear physicists) Coulomb, transverse electric or 
transverse magnetic matrix elements [65-67,60]: 


IEi = 27r |T"p + |T" 


(30) 
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(31) 


Wq = 27r 


(A2 _ q2y 


2 |M^P - 


A 2 - 


(|T"| 2 + |T“p) 


Inelastic photon processes on nuclei are dominated by the giant dipole reso¬ 
nance (GDR). For them the Goldhaber-Teller model gives [60] 


IFi = 

W2 = 


NZ A 
A 2mN 
NZ 


A 2m AT A 


Fez(A2 - 
Fel{A^ 



(32) 

(33) 


where F^i is the elastic form factor of the nucleus, rriN the nucleon mass and 
the excitation energy A is given experimentally by 


80MeV 

A = EgdR = ^1/3 


(34) 


For protons the lowest excited state is the A-resonance. For this transition the 
contribution to the photon spectrum was estimated in [15], using the param¬ 
eters given in [68]. A contribution of the order of 10% was found. This is in 
contrast to the similar situation in heavy ion collisions, where the contribution 
of the most important excited state, the GDR, is in general less than 1% [60]. 

For higher excitation energies and momentum transfers on nuclei the quasielas¬ 
tic emission of nucleons is dominating and the quasifree approximation can be 
used. For an estimate the Fermi gas model can be used. One obtains a simple 
expression in terms of the elastic nucleon form factors [69,70]: 

H-r = C(() \ZW! + (A- Z)HT], (35) 


with 


C(t) 



Qrec > 2Pf 
otherwise. 


(36) 


with + the Fermi momentum Pp = 0.25 GeV. 

Finally at very large momentum transfer {\—q^\ 3> IGeV^) the parton model 
applies and Wi and can be expressed in terms of the quark structure 
functions fi{x) [71,72]: 
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( 37 ) 


I 

W2 = ^^elxfi{x), (38) 


that is, Wi and W 2 are related by 2MWi = UW 2 /X and the Bjorken variable is 
X = —q^/{2mv). Of course in this case the integration over dv can be replaced 
by one over dx. In the case of the proton, the equivalent photon spectra were 
calculated in [71,72], Whereas the integration over in the elastic case is 
limited due to the form factor to small values of here the integration goes 
up to the kinematical limit (—< 4i7^). 


2.5 Integrated Photon-Photon Luminosity 


Let us now calculate the luminosity function dL.y.y/dW as dehned in Eq. (9). 
Neglecting possible “initial state effects” between the collision partners, we 
can fold the equivalent photon spectra ni(a;i) and n 2 {uj 2 ) to obtain [32] 


dL 




(39) 


where r = W'^/s = A 001002 /S is the square of the fraction of the total energy 
carried away by the photon-photon system with s = {pi -\-p 2 )‘^- The photon 
distribution function / here is related to the equivalent photon number via 


fix) 


— nixE) 

OJ 


n{xE) 

X 


(40) 


and E is the total energy of the initial particle in a given reference frame, 
a; = cu/E is as before the ratio of the photon energy uj to the energy of the 
incoming particle E. 

Explicit calculations of 77 luminosity functions exist for the pp case. The 
strong interaction between the protons are not taken into account. Several 
approximate results are given in [72]. The simplest one (in leading order of 
ln(r)) is given by 


dL 




(41) 


A more detailed result is given by Drees et ah [73] based on the equivalent 
photon spectra of the proton given in [74], 
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Such a formula for the 77 -luminosity (Eq. (39)), in which the initial state 
interactions are neglected, can be applied to 77 interactions in ep and eA 
collisions. For eA collisions this was done by Levtchenko [75]. Eq. (25) was 
used for the equivalent photon spectrum of the nucleus. In ep and eAcollisions 
initial state interactions are weak and can be safely neglected. 

In the case of protons calculations have been made with either one or both pro¬ 
tons breaking up [71,72]. This leads to 77 luminosities which are larger than 
the ones in the elastic-elastic case. Since the charges of the quarks (2/3, —1/3) 
are comparable to the proton charge, this is not unexpected. In the case of a 
heavy nucleus on the other hand the charges of the constituents (protons) are 
mnch smaller compared to the total charge Z of the ion; therefore incoherent 
effects, althongh present, are generally less important. Possible strong interac¬ 
tion effects in the initial state are not taken into acconnt in these calculations. 

The interesting possibility to stndy photon-photon events in pp collisions at 
the LHC by tagging the two hnal protons was studied in [76]. Only protons 
emitting a photon with a snbstantial part of their energy (corresponding to 
a; > 0 . 1 ) can be tagged in this way, as their path after the interaction deviates 
snfhciently from the non-interacting ones. Photon-Photon luminosities nnder 
these conditions (and also inclnding transverse momentnm cnts, which are less 
important) were calculated. 


2.6 Semiclassical Impact Parameter Description 


For heavy ions (as opposed to, e.g., pp collisions) the semiclassical descrip¬ 
tion works very well, as the Sommerfeld parameter p = is mnch larger 

than one, see, e.g., the discussion in Chap. 2 of [15]. Therefore an impact 
parameter dependent eqnivalent photon nnmber N{u,b) can be dehned. It is 
most important to note that in this semiclassical impact parameter descrip¬ 
tion it is very easy to take strong absorption effects into acconnt. This is done 
by introdncing a cntoff for those impact parameters where the ions interact 
strongly with each other. The calculation of the impact parameter depen¬ 
dent eqnivalent photon spectra is explained in [7] for the case of El (electric 
dipole) excitations. The generalization of the eqnivalent photon spectrnm to 
all electromagnetic mnltipoles was given by [40]. For the El case 


b) 


Z'^a / cn 

TT^ \7/9 / 


KKu) + 

/yz, 


(42) 
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with (3 and 7 the velocity and Lorentz factor of the ion and where u = ub/{'yj3). 
Integrating from some Rmin to inhnity gives 


UEliuj) 






(43) 


where ^ = uRmin/lP- Using the expression for Kq and Ki for 0, it can 
easily be seen, that this expression agrees in logarithmic approximation with 
Eq. (11). Impact parameter dependent eqnivalent photon spectra N{u,b) for 
extended charge distribntions are calculated in [77]. They are 


N{uj,b) 


Z^a 1 

TT^ 


00 

J dv n^Ji(n) 

0 


J^el - —) 


(44) 


where u = ojb/'y(3 as above. In the case of a monopole or dipole form factor 
analytical results can be found [78,79]. 


2 .7 Effects of Strongly Interacting Particles 


Let us now turn to the effects of strong interactions between the colliding 
particles (“initial state interactions”). For the case of nucleus-nucleus (AA) 
collisions we can use the semiclassical method: to a very good approximation 
the nuclei move on classical straight line trajectories with an impact parameter 
b. Collisions with b < Rmin = R1+R2 are dominated by strong interactions and 
electromagnetic effects (even though they still occur) are generally completely 
swamped by those violent processes. As an example for the presence of these 
electromagnetic effects in central collisions we mention Ref. [80]. In this paper 
production due to the electric helds of the ions is calculated. Their effect 
was found to be much smaller than e'''e“-pair production due to hadronic 
processes and from meson decay [81]. Usually, the experimental conditions 
imposed on the study of electromagnetic processes ( 77 , 7/1 and 7 A) exclude 
these collisions, as they cannot be extracted from hadronic events. Of course 
one has to be careful to exclude them in the same way in the theoretical 
study of these processes as well. One generally imposes a (sharp) cutoff on the 
impact parameter 


b ^ Rmin Rl T R2 • 


(45) 
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In the one-photon case one then gets the equivalent photon number as 

OO 

n{ijj) = J 2TrbdbN{u, b). (46) 



Fig. 5. In the semiclassical picture, which is valid in the case of heavy ion collisions, 
initial state interactions between the two ions take place if the impact parameter 
between the two ions is smaller than Rmin = Ri + R 2 - Final state interaction can 
occur if the individual bi are smaller than R^. 

Let us look first at the one-photon case (“ 7 A collisions”) between two ions. 
It is a good assumption that the total nuclear charge is contained inside the 
sphere with radius i?*. As it is well known from electrostatics, the electric 
held outside of a spherical symmetric charge distribution is the same as if the 
total charge is concentrated in the center. The two ions cannot come closer 
than Rmin = Ri + R 2 - Impact parameter dependent equivalent photon spectra 
for extended charge distributions can in principle be calculated. But from the 
argument given above it is clear that the form of the charge distribution does 
not enter in this case in Eq. (46). It only matters that the nuclear charge Z is 
entirely contained inside the nuclear radius Ri. 

The corresponding electromagnetic process can be written as 


(Top — 


R, 


'TTlZn 

OO 


2'Kbdb ( —Niuj,b)aJuj) 
J u 

27rbdbPop{b), 


Rrr 


(47) 

(48) 


see Eq. ( 6 ), allowing to dehne the impact parameter dependent probability 
for the one-photon process in this case. In numerous studies of nuclear ex¬ 
citations of fast (relativistic) heavy ions this has been proven to be a good 
approximation. 


The sharp cutoff assumption may be relaxed using a smooth absorption prohle 
T( 6 ), see, e.g.. Chap. 2 of [15]. The importance of such uncertainties depends 
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on each case: for dipole excitations the dependence on Rmin is only logarithmic, 
for higher multipolarities it depends on an inverse power of the cutoff Rmin- So 
in the former case the uncertainties in the choice of a cutoff will be small. For 
RHIC the effect of the surface thickness was calculated and found to contribute 
to about 5-10% [49]. An interesting question is the choice of a cutoff in the case 
of deformed nuclei. Depending on the orientation of the nuclei, the minimum 
impact parameter varies. It would be interesting to study this effect in the 
future. For a related problem we refer to [82]. 

Let us mention the attempt by Benesh, Haynes and Friar [83] to include 
the effects of the i?mm-cutoff by using plane waves, a nuclear form factor, 
and some cutoff qmax (corresponding loosely to a minimum impact parameter 
Rmin)- From the above it is clear that this cannot work quantitatively, as is 
explained in detail in a comment by Baur and Bertulani [84], see also [85]. 
The equivalent photon spectrum used in [83] corresponds to the ’’point-form” 
spectrum in Fig. 2.3 of [16]. In addition the results of [83] were compared to 
experimental data for giant resonance excitations in [ 86 ]. 

For the 77 luminosity in the two-photon (TP) case for AA collisions, modi¬ 
fications due to the strong interaction are important especially for the high 
energy end of the luminosity function. Again we apply the condition that the 
minimum impact parameter between the two ions must be larger than Rmin- 
The photon-photon luminosity can then be calculated as [30,32] 


dL 


77 


dWdY 


^Jd^JdX 

b^)N2{^e-Y, h2)Q{\h + h 


Rn 


(49) 


Here, especially for non-strongly interacting hnal states, one does not neces¬ 
sarily have a condition for the individual bi to be larger then i?*. Therefore 
form factor effects might become important in contrast to the 7 A collisions 
considered above. 


This formula has also been derived ab initio within the semiclassical approx¬ 
imation in [87,88]. Effects due to the photon polarization are also included. 
They are neglected here for simplicity, as they were found to be small. A simi¬ 
lar approach for the derivation of the impact parameter dependent luminosity 
L^.y was done in [89]. There an expression not in terms of the individual bi, but 
a direct impact parameter dependent luminosity function L^^{b) or n{(jJi,uj 2 , h) 
was derived. (Of course this expression can also be transformed into a formula 
dependent on the individual &*, as explained in [16]). 

For a calculation of the 77 luminosities it is easiest to calculate first the 77 
luminosities without the restriction b < Rmin = Ri + R 2 , which can be done 
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analytically, e.g., for point particles with the restriction bi > Ri. The correct 
result will then be obtained by subtracting those cases, where b < Rmin- One 
obtains [32] 


dWdY W \2 ) \2 ) dWdY' 


(50) 


with 


d 

dWdY 


„ oo bl+Rrnin 

-Ri max(fei-Rmi„,R2) 


&2 db2 


xiv(ye^,6i^iv(ye ^,&2)0crit- (51) 

The integral over 0 goes only from 0 to and from 27r — to 27r. Only 
in these cases do we have contributions from b < Rmin- 0crit given by 

0crit = arccos ^ 26162^"" ) ’ 


For 61 —cx) we also have &2 —^ and ^^rit integral converges 

rapidly. It was also noted in [32] that for symmetric collisions the 77 luminosity 
can be written in the form 


T 


dL 


77 


dr 


Loi{z), 


with 


Lq — 


37r2 


(53) 


(54) 


and a universal function ,^(;7) of ;7 = 2MR^ = mR/'j. A parameterization 
of this universal function is given there as well. The case 2: -C 1 was first 
discussed in [27]. 

The last modihcation concerns possible hnal state interactions of the produced 
particles with the colliding particles. If the produced particles are hadrons, 
they can interact with the ions, if they are produced “within” one of the ions 
bi < Ri or 62 < -^2- Such processes are excluded in the approach given above, 
as integration over bi and 62 are starting from Ri and R 2 , respectively. Since 
the electric held strength inside the nucleus decreases with decreasing radius 
such effects are expected to be small. They will be of importance only at 
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the very high end of the 77 spectrum. A calculation taking into account the 
hnal state interaction has to our knowledge not been performed. But one can 
obtain an estimate (or an upper bound) of the size of this effect by comparing 
77 luminosities with either an integration over all bi or restricted to bi > Ri 
respectively. Of course for this one needs to use impact parameter dependent 
equivalent photon spectra with an elastic form factor. 


2.8 Effective -Luminosities and Perspectives for RHIC and LHC 


The effective 77 luminosity L^s in ion-ion collisions is dehned in terms of the 
the beam luminosity Laa as 


dLgg 

'dM 


j. d L/ 
dM ' 


(55) 


where the 77 luminosity L.y.y is given by Eq. (49). Luminosities of the heavy 
ion beams are several orders lower compared to those of light ions and pro¬ 
tons. One reason is the large cross sections of electromagnetic processes of 
heavy ions, which either disintegrate the ions or change their charge state due 
to electron capture. In [90] the influence of different beam-beam interaction 
processes on the beam lifetime is considered for the LHC. The main processes, 
which contribute to the beam-beam interactions, are hadronic nuclear interac¬ 
tions, electromagnetic dissociation, where an ion is excited and subsequently 
decays, and bound free electron-position pair production, see below. Thus, the 
maximum ion luminosity is derived from the cross sections of the beam-beam 
interaction. The luminosities for different ion species at the RHIC and the 
LHC are given in [91] and [92], respectively. In Table 1 we quote the average 



Projectile 

Z 

A 

Vi, A GeV 

Luminosity, cm ^s ^ 

L 

P 

1 

1 

14000 

1.4 • 10^^ 

H 

Ar 

18 

40 

7000 

5.2 • 10^9 

C 

Pb 

82 

208 

5500 

4.2 • 10^6 

R 

P 

1 

1 

500 

1.4 • lO^i 

H 

Cu 

29 

63 

230 

9.5 • 1027 

I 

Au 

79 

197 

200 

2.0 • 10^6 

C 







Table 1 

Average luminosities at LHC and RHIC for pp, medium and heavy ion beams, 
luminosities at RHIC and LHC for pp, medium and heavy ion collisions. The 
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specific experiments may also apply their requirement to lower the luminosities 
to satisfy the detector load [93]. 



W(yy), GeV/c^ 



Fig. 6. Effective 77 luminosity at LHC (left) and RHIC (right) for different ion 
species and protons as well as at the e"*'e“ collider LEP-II. The Laa iu the last two 
cases of course corresponds to the pp or ee beam luminosity. 

Two-photon luminosities for various ion species and protons as a function of 
77 mass are shown in Fig.6 (left) for LHC and (right) for RHIC. The 77 
luminosities are calculated by the Monte Carlo program TPHIC [54]. The 
luminosities at the ion colliders are compared with the 77 luminosity at LEP- 
II with a c.m. energy of 200-GeV and an e’''e“ luminosity of 5 ■ 10^^ cm“^s“^ 
|94], 

Although the heaviest ions generate the highest 77 luminosities due to their 
large charge, the beam luminosities are lower for these ions. As a result, the 
medium-weight ions (Argon at LHC, Copper at RHIC) are the most prominent 
source of two-photon processes. RHIC can compete with LEP in two-photon 
studies at low 77 masses while LHC will be the best machine to study two- 
photon physics at all ranges of 77 masses. 



































3 77 Physics at Hadron Colliders, General Considerations 


We will now give a general discussion of possible photon-photon physics at rel¬ 
ativistic heavy ion colliders. Two-photon interactions allow to test the main 
properties of the standard model in both the electroweak sector and in QCD, 
as well as, physics beyond the Standard Model. In contrast to hadronic pro¬ 
cesses, many of the processes here can be calculated in principle to some de¬ 
gree of accuracy. Below we review QCD interactions and also physics related 
to electroweak interactions, which are of interest to be studied in two-photon 
interactions in heavy ion collisions. This covers the range starting from exclu¬ 
sive meson production, meson pair production up to the total hadronic cross 
section. In the next chapter we will then be concerned with the potential for 
the discovery of new physics, which seems to be possible in principle due to the 
high two-photon invariant masses of up to about 100 GeV, which are avail¬ 
able at the LHC. An interesting topic in itself is the electron-positron pair 
production. The helds are strong enough to produce multiple pairs in a single 
collision. A discussion of this subject will be given in Chap. 7 below. 

Up to now photon-photon scattering has been mainly studied at e’''e“ colliders. 
Many reviews [59,95,96] as well as conference reports [97-100,57,58] on this 
subject exist. Whereas in the past the range of invariant masses has been 
the region of mesons, ranging from 7 r° (m^rO = 135 MeV) up to about rjc 
{rrir,^ = 2980 MeV), the higher invariant masses at LEP2 have allowed to study 
an invariant mass range up to about 185 GeV [101], where one interesting 
subject has been the study of the total 77 —> hadron cross-section. 

We are concerned here mainly with the invariant mass region relevant for the 
LHC; see the 77 -luminosity in Fig. 6 . At RHIC 77 physics can be done in 
the 77 invariant mass region up to several GeV. This topic has also been 
studied experimentally now by the “Peripheral Gollision Group” at STAR 
and this will be discussed below in Ghap. 8 . Apart from the production of 
(and pairs, the photons can always be considered as quasireal. 

The cross section for virtual photons deviates from the one for real photons 
only for Q^, which are much larger then the coherence limit ^ see 

also the discussion in [59]. Therefore photon-photon processes are restricted 
to quasireal photons. Tagging of the elastic scattered beams is very difficult 
in the heavy ion case, but for the pp case at the LHG such a possibility has 
been investigated in [76]. The standard detector techniques used for measuring 
very forward proton scattering will allow a reliable separation of interesting 
77 interactions 
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3.1 The Photon at High Energy 


Real photons have a complicated nature. In a first approximation the photon 
is a point-like particle, although in held theory it may huctuate also into a 
fermion pair [102], Fluctuations into a quark-antiquark pair 'j ^ qq can in¬ 
teract strongly and give a large contribution to the two-photon hadronic cross 
section, while huctuations into a lepton pair 7 —*> // interact electromagneti- 
cally only and do not inhuence the total 77 cross section very much (see also 
the rates predicted for the different contributions below for the LHC). Lepton 
huctuations can be calculated perturbatively, which is not true for the quark 
pairs, as low-virtuality huctuations need to be described by nonperturbative 
QCD. Therefore, the spectrum of the photon huctuation can be separated 
into low-virtuality and high-virtuality parts, according to [102]. To describe 
the hrst part, a phenomenological model of vector-meson dominance (VMD) 
is used, according to which the photon huctuates into a sum of vector me¬ 
son states. Quark pair with high virtuality are described by the perturbative 
theory. As a whole, the photon wave function can be written then as follows 



(b) 


Fig. 7 . The photon can fluctuate into 
other states, among them into a pair 
of fermions (leptons, quarks, (a)), but 
also into vector mesons (b). The total 
cross section in photon-photon collisions 
is dominated by the fluctuations into 
two quarks and their strong interaction. 


7 ) — Co I 7o) + Z]v=pO,a;,0, J/7,T Cy I V) + Y.q=u,d,s,c,bCq I 5'Q') + 

C( I ^ ^ )) 


(56) 


where | 70 ), | R), | qq) and | Hl~) are wave functions of the point-like photon, 
a vector meson, a quark- and a lepton pair, respectively. The coefficients c* in 
general depend on the scale n, which probes the photon. The coefficients for 
the contribution to a lepton pair are well known to be equal to 

cz ~ ^ ^ 

To separate the low- and high-virtuality quark fluctuations one introduces a 
parameter pq. For high-virtuality fluctuations one sets 

Cq ~ (aem/27r)2Qjln(/iVPo)- 

The qq fluctuations with virtuality below Pq are described by the vector meson 
dominance part, which does not depend on the scale p. As a rule of thumb 
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the scale jj, is taken equal to the transverse momentum of the parton process. 
The value of po, which provides the best description of the total cross section, 
is found to be po = 0.5 GeV [102]. The cy are dehned as 

Cy = 47raem//y, 

where the decay constants /y/ 47 r are determined from experimental data on 
V —they are fyj^'n = 2.20 for 23.6 for cy, 18.4 for 0, 11.5 for J/'^. 
Finally the cq is then dehned from unitarity, and is usually close to one. 

According to the photon wave function representation of Eq. (56), the fol¬ 
lowing partonic processes can take place in two-photon interactions, see also 
Fig. 8 : 

• Purely electromagnetic processes of order 0{al^) involving quarks q and 
leptons /, i.e., 77 — qq, 7 / — 7 /, 77 — ll, W ll', as well as, the production 
of W~^W~ pairs and of pairs of particles beyond the Standard Model. 

• Processes with one electromagnetic and one strong vertices of order O (cremes) 
such as 7 g —gq and 'jg qq. 

• Strong processes involving quarks and gluons of order O^a"^) qq' qq', 
qq q'q', qq gg, qg qg, gg qq, gg gg. 

• Nonperturbative QCD processes (elastic, diffractive scattering, normally at 
low transverse momentum). 



Fig. 8 . Diagrams showing the contribution to the 77 ^hadron reaction: direct mech¬ 
anism (a), vector meson dominance (h), single (c) and double (d) resolved photons. 

To describe those two-photon processes where quarks and gluons interact, one 
uses the structure function of the photon, which has been measured experi¬ 
mentally (see, e.g. [103,104]). The structure functions of the vector mesons are 
less well known, therefore the approximation |p°) ^ | 7 r°) ~ (Itt"'") -|- | 7 r “))/2 is 
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often used, where the structure functions of pious are known [105]. Below we 
give the cross section formulae of the elementary processes within the leading 
order of the perturbative theory and describe some nonperturbative models, 
which involve photons. 


3.2 Total 77 Cross Section 


The main contributions to the total 77 collision comes from the strong inter¬ 
actions. The total cross section 77 —hadrons can be parameterized in the 
form predicted by Regge theory 

a^:,^A{s/sor + B{s/so)-\ (57) 


where s is the square of the invariant 77 mass, sq = IGeV'^, and the exponents 
e and 7 are universal parameters, that is, they are identical to those in yp, pp, 
but also np total cross sections, corresponding to Pomeron and Reggeon ex¬ 
change, respectively. The universal values for these exponents from a combined 
£t of all cross section are e = 0.093(2) and p = 0.358(15) [106]. 

This total cross section is an important input in the study of high energy 
e’''e“ collisions, see, e.g., the review of Peskin [107]. It is also interesting to 
understand what part comes from point-like processes and what part from 
soft processes involving the hadronic constituents of the photon. Eventually 
theory should explain, as well as, be constrained by the data of both < 7 ( 77 ) 
and ct( 7 p). As also done in [ 101 ] there are in general two fits to the total cross 
section data: one with PHOJET and the other one with the PYTHIA event 
generators. They differ substantially, and this is due to the fact that about 
40% of the cross-section is unobserved at LEP and that theoretical models 
differ considerably in the size of the contribution from these very soft events. 

The best fitted values before LEP2 were as follows [102]: 

A ^ 211nb, B ^ 297nb, e 0.0808, p ^ 0.4525. (58) 


The L3 collaboration recently made a measurement of the total hadron cross- 
section for photon-photon collisions in the interval 5GeV < < 185GeR 

[108,101]. Fitting the data up to 65 GeV, it was found that the 77 —i>hadrons 
cross-section is consistent with the universal Regge behavior of total hadronic 
cross-sections. Values of 

A = 173 ± 7nb, R = 519 ± 125nb, e = 0.0790, p = 0.4678 (59) 
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were given in [108] (where the values of e and 77 had been held hxed). Their 
cross section were found to be in agreement with the universal values of rj 
and e. 

Using the larger invariant mass range up to 185 GeV in [101] a deviation from 
the universal value for e was found. Fitting the data with rj = 0.358 fixed gave 
values of 

A = 59 ± lOnb, B = 1020 ± 146nb, e = 0.225 ± 0.021, (60) 


that is a value for e more than twice as large as the universal one. This large 
exponent was also found to be independent of the Monte Carlo model used to 
correct the data. 


3.3 Charged Fermion Pair Produetion 


For the invariant mass of the 77 system above the threshold ^/s > 2mf a 
fermion (lepton or quark) pair can be produced in two-photon collision. In 
lowest order of QED this process is shown in Fig.9. The differential cross 



section of the process 77 — f~^f in the c.m. system is given by the equation 
[109,110] 

da ^ e'^PQjNc 1 + 2(3'^{1-6) - 6 

dcoso^^^ ^ ^ ^ ~ StTS (1-/32cOs20)2 ’ ^ 


where Qf is a fermion charge in units of electron charge, in a number of 
colors, (3 = is the velocity of the fermion in the 77 rest frame. 

The integral cross section of this process is 


c^(77 ^ ) 


- - —/5 

s 


3-/34 

2/3 


In 


\ (3 


-2 + /32 


(62) 
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The dependence of the integral cross section on the ratio of the fermion mass 
to the beam energy in the c.m. system, is shown in Fig.lO. 



Fig. 10. Dependence of the process 
77 ^ f^f~ cross section on 


The production of fermion pairs in two-photon collisions can be used to study 
the coupling of s-, c- and 6 -quarks to photons and to study the fragmentation of 
these quarks into iF, D and B mesons, because these processes are clean from 
hadronic background. To normalize the two-photon luminosity the process 
77 —can be used as it is easy to observe and simple to calculate 
[35,111,112]. In Table 2 production cross section of ss, cc and bb are 

shown in different invariant mass intervals in collisions of various species of 
ions at the RHIC and LHC. Figure 15 shows cross section and event rates for 
the LHC. 


/+/- 

Mass range, GeV 

1 

T 

b 

.AA+f+n fih 



AuAu 

CaCa 

PbPb 


1 < IF < 10 

280 

730 

1.1 • 10® 

ss 

1 < IF < 10 

7.6 

20 

3.1 • 10^ 

cc 

3.7 < IF < 10 

1.5 • 10-2 

6.2 

790 

bb 

10.6 < IF < 20 

— 

1.3 • 10-2 

1.2 


Table 2 

Production cross sections of fermion pairs in two-photon interactions in AuAu col¬ 
lisions at y/SAA = 100 GeV per nucleon at the RHIC, as well as in CaCa and PbPb 
collisions at ^saa = 7.2 TeV per nucleon ^saa = 5.5 TeV per nucleon at the LHC. 
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Fig. 11 . Cross section per GeV for different charged fermion pair production at 
the LHC for PbPb (a) and CaCa (b) collisions, using the lowest order QED cross 
section. Also shown is the total rate for 77 ^hadrons, see Eq. ( 57 ). Also shown are 
event rates per s and per (10®) year. 

3.4 W^W~-Pair Production 


If the energy of the two-photon collision is higher than twice the mass of 
LF^, the process 77 — W^W~ is possible. This process involves the gauge 
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W 



Fig. 12. Process 77 ^ W~^W in the lowest order. 

couplings of these gauge bosons . The lowest order graphs are shown in Fig. 12. 
The differential cross section of this process is given by 


da 


dcos9 


(77 ^ W+W-) = 


TTa^/3 19 — 6/5^(1 — + 2(8 — 3/3^)/3^ cos^ 6 + 3/3^ cos"^ 6 

(1 — cos^ oy ’ 


( 63 ) 


with the same notations as those in Eq. (61). The integral cross section is [113] 


a(77 ^ W+W-) 




In 


1 + /3 


, ^22-9/32 + 3/3^ 


(64) 


The magnetic moment, as well as, the quadrupole moment of the W are hxed in 
the standard model, for the more general case of a magnetic moment different 
from the one given by the standard model see also [114]. 

Fig. 13 shows the dependence of the integral cross section on the ratio of a 
mass to the beam energy in the c.m. system, 2mw/\/s- For example, the cross 



Fig. 13. The process 77 ^ W~^W cross 
section as function of 2 mw / \/s- 


section estimates of the W~^W production in CaCa and PbPb collisions at the 
LHC energies are a{AA AA + W~^W~) = 29 pb and 190 pb, respectively. 
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3.5 Vector Meson Pair Production 


Among the nonperturbative processes in photon-photon interactions there is 
the exclusive vector-meson pair production. There are various mechanisms to 
produce hadrons in such collisions, recall the discussion about the nature of 
the photon above. According to this idea, vector meson pair production is 
explained to happen in the following way: both photons fluctuate first into a 
vector meson (p, u, 0, .J/ijj and so on, the VMD component). These vector 
mesons then scatter elastically to produce the (real) vector meson pair in the 
hnal state. The study of these processes allows a test of the Pomeron-exchange 
factorization relation, see [115,116]. If the exchange of a single Pomeron dom¬ 
inates in these processes, the factorization relation among pp, py and 77 col¬ 
lisions predicts 


lffioi(7P)]^ = fftoi(pp) X ffioi(77)- 


(65) 


As for the available experimental data there exist measurements of the produc¬ 
tion cross section of vector meson pairs in 77 -interactions by the ARGUS 
[117] and pp by the L3 [118,119] collaboration. These data were compared with 
the cross section estimate of the reaction 77 — W (U, V = p°, 0 ) based on 
the Pomeron factorization model and all possible combinations of the exist¬ 
ing sets of data on the reactions 7 p — Vp and pp pp [ 120 , 121 ]. A strong 
discrepancy between the existing experimental data (by order of magnitude) 
and the factorization model prediction was found. These estimates are in good 
agreement with independent calculations in [102]. Thus we come here to the 
puzzling situation: why does one of the most well-grounded phenomenological 
models predict cross sections, which are larger by order of magnitude in the 
reaction 77 —> VV compared to those measured in the experiments? One 
explanation proposed in [ 121 ] is that we face here a new '^defiant phenomenon 
in the formation mechanism of the Pomeron exchange for quasi-two-hody re¬ 
actions" . 

The production of vector meson pairs can well be studied at RHIC with high 
statistics in the region of up to several GeV [9]. In connection with the above 
mentioned puzzle the results from the STAR would be quite desirable and 
interesting. For the possibilities at LHG, we refer the reader to [24] and [22], 
where also experimental details and simulations are given. 


3.6 Resonance Production 


One may say that photon-photon collisions provide an independent view on 
meson and baryon spectroscopy. They provide powerful information on both 
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the flavor and spin/angular momentum internal structure of the mesons. Much 
has already been done at colliders. Light quark spectroscopy is very 

well possible at RHIC, benefltting from the high 77 -luminosities. Detailed 
feasibility studies exist [9,46-48]. In these studies, 77 signals and backgrounds 
from grazing nuclear and beam gas collisions were simulated with both the 
FRITIOF and VENUS Monte Carlo codes. The possibilities to produce these 
mesons at the LHC have been discussed in detail in the FELIX Lol [24]. Rates 
are given and possible triggers are discussed. The general conclusion is that 
all these processes are very promising tools for meson spectroscopy. 

In two-photon collisions with real photons general symmetry requirements re¬ 
strict the possible final states, as is well known from the Landau-Yang theorem 
[122]. Especially it is impossible to produce spin 1 final states. Only resonance 
states with positive C-parity can be produced, such as = O”’*', O’*"’'', 2 ++, 
■ ■ -. Two photon collisions therefore give access to most of the C = -|-1 mesons. 
In e'''e“ annihilation on the other hand only states with = 1 , that is, 

with odd C-parity can be produced directly. 

In principle C = — 1 vector mesons can be produced by the fusion of three 
(or, even less important, five, seven, ...) equivalent photons. This cross sec¬ 
tion scales with Z®. But it is smaller than the contribution coming from 7 -A 
collisions, as discussed in Chap. 5, even for nuclei with large Z (see also [15] 
and the corresponding discussion for the positronium in Sec. 7.7). 

The production cross section of the resonance is given by 


where W is the 77 invariant mass , Mpi is the mass of the resonance R, F..j,..|, and 
Ftot its two-photon and total width. For sufficiently narrow state (Ftot ^ Tf/j) 
the expression ( 66 ) can be approximated by 

lim = + (67) 


This makes it easy to calculate the production cross-section aAA^AA+R of a 
particle in terms of its basic properties. 


In this case the production cross section of a resonance in heavy ion collisions 
factorizes, according to equation (9), to 


a(AA ^AA + R) 


47r2(2J+ 1) 


77 


dL 


77 


dW 


W=Mji 


( 68 ) 
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i.e., it is given only by its two-photon width and the value of the 77 -luniinosity 
function at its resonance mass. 

In Fig. 14 the function /dW/W"^, which is universal for a produced 

resonances, is plotted for various systems. It can be directly used to calculate 
the cross-section for the production of a resonance R with Eq. ( 68 ). 



W(GeV/c^) 

Fig. 14. The universal function An'^dL^^/dW/W'^ is plotted for different ion speeies 
at LHC. We use R = 1.2 x fm and 7 = 2950, 3750 and 7000 for PbPb, CaCa 
and pp collisions, respectively. 


3.6.1 Quarkonia 

Among the interesting processes for resonance production is the study of the 
quarkonium states cc and bb such as ?7c(b) (^>S'o), Xc{b)o Xc{b)2 with the 

aim of measuring their two-photon widths, and therefore to test quarkonium 
models, as well as, studying their decay modes. Production of light qq states, 
like 7 r°, 7 and 7 ', whose widths are well-known, can be used for calibration. 

In the nonrelativistic quarkonium model of [123], which takes into account 
first-order QCD terms, analytical expressions for gluon, lepton and photon 
widths of various quarkonium states are obtained. 

In Table 3 these expressions for states with even total angular momentum are 
given, as well as, those for the vector state ^Si for normalization. In these 
expressions Q denotes the electric charge of a quark, constituting the quarko¬ 
nium, M is the quarkonium mass, and rriq the quark mass. For the quark 
masses we use rric = 1.5 GeV, rrib = 4.8 GeV. The strong coupling constant 
ttg is derived from the relation between the total and partial widths of the 
decay of the known states. Gomparison with experimental data gives [123] 


39 




















Process 

Width 

Correction 0(as) 


77 


1 — 3.4as/7r 


gg 

STTOg 'h(0)p/3mq 

l + 4.8as/vr (7c) 

l + d.dos/vr (7fe) 


‘ e^e~ 

IGTrg^agjjj 'I'(0)p/M^ 

1 — 16as/37r 


ggg 

40(7r^ — 9)q;s 'I'(0)p/81mg 

1 — 3.7q;s/'P (d/U) 

1 — 4.9as/'P (T) 

CO 

77 

27g^aLlKp(0)lV"iq 

1 -|- 0.2Q;s/7r 

^Po- 

gg 

6a2|i?;^(0)|Vm^ 

l + 9.5as/vr (Xco) 

1 + IO.Oos/tt (xbo) 

^P2- 

77 

36g^aLl-Rnp(0)lV5m^ 

1 — 16as/37r 

^P2- 

gg 

8a2|i?(^p(0)|2/5m^ 

1 - 2.2as/vr (Xc2) 

l-O.lOs/vr (Xfe2) 


Table 3 

Quarkonia decay widths in the nonrelativistic model. Adapted from [123]. 

asirric) = 0.19, asirrih) = 0.17. The unknowns are the value of the wave func¬ 
tions at the origin 'h(O) and R'^piQ). 


Assuming that the wave function T(0) is the same for all IS states, and equal 
to those in Table 3, the width ratio of J/'ip(T) to rjc {rjb) (l^*S'o) can be 

calculated: 


e+e ) ^ ^ f2mgY A _ q 

T{RSo ^77) 3g2 V M M 


(69) 


This ratio for J/ijj and rjc is equal to 0.74, in perfect agreement with the 
experimental result 0.71 ± 0.14 [106]. Therefore, one can predict two-photon 
width of the rji, from the lepton width of the T(IS'). The ratio of the total 
and two-photon width of the rjb is obtained analogously, assuming that the 
two-gluon decay dominates: 


r(7fe ^ gg) ^ 2a^ , 7 

r(7fe ^ 77) 9a2^g2 ^ 


(70) 


Using the experimental lepton decay width r(T(lS') —e’''e ) = 1.32 keV 
[106], one obtains from Eqs. (69) and (70) the total and the two-photon widths 
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of the rjb: 


rtot(hb) = 6 MeV, ^y'riVb) = 0.43 keV. (71) 

P-states of bottomonium are not yet studied, therefore predictions of their 
total and partial widths cannot be obtained from a comparison of experimental 
data. In [124] lattice calculations give a value for the wave function of the IP- 
state of bottomonium at the origin |P(p(0)p = 0.75 GeV^. These calculations 
allows to dehne total and two-photon widths of the Xbo and the Xb 2 from the 
analytical expression of Table 3: 


rtot(xw) = 0.24 MeV, 

r 77 (Xfeo) = 25 eV, 

(72) 

rtot(X 62 ) = 65 keV, 

r 77 (X 62 ) = 6.7 eV. 

(73) 


For charmonium production, the two-photon width of the r]c (2960 MeV, 
= 0 - '+) is known from experiment [106]. But the two-photon widths of the 
P-wave charmonium states have been measured only with modest accuracy. 
Two photon widths of P-wave charmonium states can be estimated following 
the PQCD approach of [124]. Similar predictions of the bottomonium two- 
photon widths can be found in [123]. For RHIC the study of r]c is a real 
challenge [46]; the luminosities are falling strongly with increasing 77 mass 
and the branching ratios to experimentally interesting channels are small. 

In Table 4 (adapted from Table 2.6 of [24]) properties of some qq states are 
given, and their production cross sections are predicted, where possible. Sim¬ 
ilarly an overview of different rates is also given in Fig. 11 for both PbPb and 
CaCa collisions at the LHC. The results for AuAu collisions at RHIC and 
PbPb collisions at LHC are done on the basis of [31] and from calculations of 
the program Tphic [54]. Mass values and known widths are taken from [106], 
bottomonium widths are used from Eqs. (71) to (73). Also given is the number 
of events in a 10®s run. Ion luminosities of 2 x 10^®cm“^s“^ for AuAu collisions 
at RHIC and 4 x 10^°cm“^s“^ for CaCa and 10^®cm“^s“^ for PbPb collisions 
at LHC are used. Millions of C-even charmonium states will be produced in 
coherent two-photon processes during a standard 10 ®s heavy ion run at the 
LHC. The detection efficiency of the charmonium events has been estimated 
to be about 5% for the forward-backward FELIX geometry [24], i.e., one can 
expect the detection of about 5 x 10^ charmonium events in PbPb and about 
10® events in CaCa collisions. This is two to three orders of magnitude higher 
than what is expected during the hve years of LEP200 operation. Experiments 
with a well-equipped central detector like CMS on the other hand should have 
a much better efficiency. Further details — also on experimental cuts, back¬ 
grounds and the possibilities for the study of C-even bottomonium states — 
are given in [24]. 
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M[GeV] (b) 

Fig. 15. Cross section for the production of different resonances at the LHC for 
PbPb (a) and CaCa (b) collisions, using the parameters of Table 4. Also shown are 
the event rates per s and per (10®s) year. 

3.7 Exotic Mesons 

The two-photon width of a resonance is a probe of the charge of its con¬ 
stituents, so the magnitude of the two-photon coupling can serve to distinguish 
quark dominated resonances from glue-dominated resonances (“glueballs”). 
The absence of meson production via 77 fusion would therefore be one signal 
of great interest for glueball search. In 77 collisions a glueball can only be 
produced via the annihilation of a qq pair into a pair of gluons, whereas a 
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normal gg-meson can be produced directly. Therefore we expect the ratio for 
the production of a glueball G compared to a normal gg meson M to be 


a(77 ^ G) T{G 77) 2 

- = - ~ Q 

(t(77 —M) r(M — 77) * 


where o;* is the strong interaction coupling constant. On the other hand glue- 
balls are most easily produced in a glue-rich environment, for example, in 
radiative J/T decays, J/T — ■ygg. In this process we expect the ratio of the 
cross section to be 


r(J/T ^ 7 G) 1 

r(J/T ^ 7 M) 


A useful quantity to describe the gluonic character of a mesonic state X is 
therefore the so called “stickiness” [125], dehned as 


r(j/T ^ 7X) 

r(A' ^ 77) 


(76) 


One expects the stickiness of all mesons to be comparable, while for glueballs 
it should be enhanced by a factor of about l/a^ ~ 20. In the recent work of 
[126] results for the search for /j( 2220 ) production in two-photon interactions 
were presented. A very small upper limit for the product of was 

given, where Bk^Ks denotes the branching ratio of its decay into KgKg. From 
this it was concluded that this is a strong evidence that the /j( 2220 ) is a 
glueball. 

Two-photon processes can also be used as a tool to observe mesons beyond 
the quark model. The ao(980) and /o(980) could be four-quark states gggg (or 
“quarktets”) [127]. The ARGUS collaboration observed a = 2+(2++) 

peak in the reaction 77 —pp near threshold [128]. This state was called by 
PDG X(1600) and also interpreted as a four-quark state. 

For a rather general discussion of glueballs but also other (even more) exotic 
mesons like “quarktets”, “hybrids” (or “centauros”, made of a quark, an an¬ 
tiquark and a gluon) we refer to [129]. The QGD studies in 77 collisions at 
the GLEO detector were summarized recently by Savinov in [130]. Also the 
anti-search for glueballs is described. 
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particle 

M, 

MeV 

r 

L 77, 

keV 

r tot: 

MeV 

a{AA AA + R) 

Production rate per 

10®-sec run 

AuAu, 

7 = 100 

PbPb, 

7 = 3000 

CaCa, 

7 = 3700 

AuAu, 

7 = 100 

PbPb, 

7 = 3000 

CaCa, 

7 = 3700 

0-+ 

TT® 

135 

8■10-3 

8■10-® 

5.0 mb 

46 mb 

210 /ib 

1.0 - 10® 

4.6 - 10® 

8.4 - 10® 

0-+ 

V 

547 

0.46 

1.2 ■ 10-3 

0.85 mb 

20 mb 

100 /ib 

1.7-10® 

2 - 10® 

4.0 - 10® 

0-+ 

i 

958 

4.2 

0.2 

0.59 mb 

25 mb 

130 /ib 

1.2 - 10® 

2.5 - 10® 

5.2 - 10® 

2++ 

/2(1270) 

1275 

2.4 

185 

0.41 mb 

25 mb 

133 /ib 

8.2 - 10^ 

2.5 - 10® 

5.2 - 10® 

2++ 

02(1320) 

1318 

1.0 

107 

0.14 mb 

7.7 mb 

49 /ib 

2.8 - 10^ 

7.7 - 10® 

2.0 - 10® 

2++ 

/'(1525) 

1525 

0.1 

112 

6.6 /ib 

0.45 mb 

2.9 /ib 

1.3-103 

4.5 - 10^ 

1.2 - 10^ 

0-+ 

Vc 

2979 

7.4 

13.2 

1.8 yub 

0.54 mb 

3.7 /ib 

360 

5.4 - 10^ 

1.6 - 10^ 

0++ 

XcO 

3415 

4.0 

14.9 

0.38 /ib 

0.17 mb 

1.2 /ib 

76 

1.7-10^ 

4.8 - 10® 

2++ 

Xc2 

3556 

0.46 

2.0 

0.17 /ib 

85 /ib 

0.59 /ib 

34 

8.5 - 10^ 

2.4 - 10® 

0-+ 

Vb 

9366 

0.43 

6 


0.32 /ib 

2.8 nb 


32 

1.1 - 103 

0++ 

XbO 

9860 

2.5 ■ 10-2 

0.24 


15 nb 

0.15 nb 


1.5 

600 

2++ 

Xb2 

9913 

6.7-10-3 

6.5-10-2 


20 nb 

0.18 nb 


2.0 

720 


Table 4. Prediction of production cross section and rates for a 10® sec-run of some mesons at RHIC and LHC. Masses and widths for 
mesons except bottomonia are taken from [106], bottomonium widths are given from equations (71), (72) and (73). The beam luminosities 
used are 2 • 10^®cm“^s“^ for AuAu at RHIC, 10^®cm“^s“^ for PbPb and 4 x 10^®cm“^s“^ for CaCa at LHC 



4 Two-Photon Collisions as a Tool for the Search of New Physics 


The large photon-photon luminosity and the high invariant mass range make 
two-photon collisions at heavy ion colliders also of interest for the search for 
new particles and new physics. This includes the possible production of the 
Higgs-boson in the 77 -production channel (unfortunately, this possibility is 
only marginal, see below) or new physics beyond the standard model, like 
supersymmetry or compositeness. Many studies for different extensions of the 
standard model have been performed. In the following the conclusions of some 
of these studies will be summarized, further discussions can also be found in 
|16], 

Before doing this we mention the plans to build a future e’''e“ linear collider. 
Such a linear colliders will be used for e+e“, ey and 77 -collisions (PLC, “pho¬ 
ton linear collider”). The photons will be obtained by the scattering of laser 
photons (of eV energy) on high energy electrons (up to the TeV region), see 
[131]. These photons in the TeV energy range will be roughly monochromatic 
and polarized. The physics program at such future machines was discussed , 
e.g., in [132,133]; it includes Higgs boson and gauge boson physics and the 
discovery of new particles. A recent review can be found in Parts 3 and 6 of 
[134]. The physics topics reach from Higgs boson physics to supersymmetry 
and extra dimensions. Such a collider will provide a rather clean environment 
for physics in a new energy region. This is of interest for the (far) future, 
whereas the LHC and its detectors are built at present. 

In [71,72] 77 -processes at pp colliders (LHC) are studied. It is observed there 
that non-strongly interacting supersymmetric particles (sleptons, charginos, 
neutralinos, and charged Higgs bosons) are difficult to detect in hadronic col¬ 
lisions at the LHC. The Drell-Yan and g'g'-fusion mechanisms yield low pro¬ 
duction rates for such particles. Therefore the possibility of producing such 
particles in 77 interactions at hadron colliders is examined. Since photons 
can be emitted from protons which do not break up in the radiation pro¬ 
cess, clean events can be generated which should compensate for the small 
production number. In [71] it was pointed out that at the high luminosity of 
L = 10^"^cm“^s“^ at the LHC(pp), one expects about 16 minimum bias events 
per bunch crossing. Even the elastic 77 events will therefore not be free of 
hadronic debris. Clean elastic events will be detectable at luminosities below 
10^^cm“^s“^. This danger of “overlapping events” has also to be checked for 
the heavy ion runs, but it will be much reduced due to the lower luminosities. 

Similar considerations for new physics were also made in connection with 
the proposed eA collider at DESY (Hamburg). Again, the coherent held of 
a nucleus gives rise to a factor in the cross-section for photon-photon 
processes in eA collisions [135,136]. 
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4-1 Supersymmetry Partiele Pair Produetion 


Two-photon collisions allow to study the production of particles beyond the 
Standard Model. Here we consider the processes of chargino, slepton and 
charged Higgs pair production in the frame of the minimal supersymmetric 
standard model (MSSM) [137]. 

Charginos are coupled to photons as standard fermions of spin 1 / 2 , and there¬ 
fore the cross section of the process 77 —xlxi is described by Eqs. (61) 
and (62) with Qf = l and Nc = 1. Following [138] we assume i?-parity conser¬ 
vation, i.e. continuity of supersymmetric lines. We take the chargino yf as the 
lightest observable supersymmetric particle and the neutralino y)* as the light¬ 
est supersymmetric particle. Thus, it is not observable. Sleptons, squarks and 
gluinos are assumed to be much heavier than the Xi ■ All parameters of the su¬ 
persymmetry breaking can be taken as real, i.e., CP-parity violation does not 
take place in the chargino processes. All these assumptions result in the fact 
that the yf decays to neutralino y5 and a weak duplet of fermions via a virtual 
W^, or via a sfermion or Higgs xt ^ iXiW+*, f*Jp y?P+*) ^ xUifj- 

Experimental search for supersymmetric particles impose limits on their masses 
[106]. The lightest neutralino y° was not found at the mass m < 32.5 GeV. 
The chargino yf mass is bound by m > 67.7 GeV. This restricts the possibil¬ 
ity to search for the processes of the chargino pair production in two-photon 
interactions at heavy ion colliders, since the invariant mass of the two-photon 
system must be above 130 GeV, while the spectrum of the 77 -system at LHG 
is limited essentially to ~ 200 GeV (see Eq. (3)). 

Gross sections of the chargino pair production vs. their mass in two-photon 
interactions in Pb-Pb and Ca-Ca collisions at LHG are shown in Fig. 16. The 
results are obtained with the program Tphic [54]. 


Fig. 16. Production cross section of 
xf Xi in Pb-Pb and Ca-Ca collisions at 
yJSAA = 5.5 and 7.2 TeV per nucleon 
respectively. 
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Fig. 17. Process 77 ^ S~^S in the leading order, is either a slepton or 
charged Higgs 


When the mass of sleptons and charged Higgs are not too high, their 
pairs can also be produced in the two-photon interactions. Sleptons and Higgs 
are scalar particles and the differential cross sections of the processes 77 — 
/+/'- and 77 ^ H+H-: 


da{'j'y S~^S ) 
d cos( 6 ') 



2 (1-/3^) 

1-/3^ cos‘^{9) 


(1 — /32 cos2(6'))^ 


, (77) 


as well as, the total cross section 


a(77 ^ S+S ) 



2-P^- 


l-/3^ 1nl + /3 

2p ^ I-P 


(78) 


are well known, see [139,72]. Fig. 18 shows how the integral cross section 
depends on the ratio of a mass to the beam energy in c.m. system, 2ms/y/s. 
Comparing this plot with the cross section behavior of fermions (Fig. 10 for 



Fig. 18. The 77 ^ S~^S cross section 
as a function of 2 ms / \/7- 


Nc = I, Qf = 1) shows that the yield of supersymmetric scalars is lower than 
that of charginos. 
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The search for sleptons and charged Higgs imposes a lower limit on their mass 
m{H^) > 69.0 GeV, m(i^) > 67.1 GeV [106]. 

Recent (unpublished) studies done for FELIX and ALIGE show that the 
chargino pair production can be detectable, if the lightest chargino would have 
a mass below 60 GeV. Unfortunately the chargino mass limits set recently by 
LEP experiments already exclude the existence of charginos below 67.7 GeV 
with 95% conhdence level [106]. Therefore the observation of MSSM-particles 
in 77 -interactions in heavy ion collisions seems to be hard to achieve. GMS 
on the other hand should be more suited for these observations. 


4-2 Higgs Search 


One of the most important tasks in present high energy physics is the discovery 
and the study of the properties of the Higgs boson. In this respect, the 77 
production of the Higgs boson is of special interest. 

There are a number of calculations of the 77 production of a medium heavy 
standard model Higgs [73,140,31,141,142]. For masses mn < 2,mw± the Higgs 
boson decays dominantly into bb. In [141] a comparison of gg ^ H and 77 — 
H emphasized the favorable signal to background ratio for the latter process. 
Unfortunately, at the LHG a heavy ion year consists only of the order of 10® 
s instead of the assumed 10^ s. Ghances of Ending the standard model Higgs 
in this case are marginal [ 22 ]. 

The search for anomalous Higgs couplings in peripheral heavy ion collisions 
at the LHG was studied by Lietti et al. [143]. They consider corrections to the 
standard model from new physics associated with a high energy scale A. It is 
concluded that 'limits for anomalous Higgs couplings which can he obtained 
in peripheral heavy ion collisions at the LHC via electromagnetic processes 
are one order of magnitude tighter than the limits that can he obtained in the 
upgraded Tevatron and comparable to limits coming from the pp mode of the 
LHC. For further details we refer to this paper. 

An alternative scenario with a light Higgs boson was, e.g., given in [144] in 
the framework of the “general two Higgs doublet model”. Such a model allows 
for a very light particle in the few GeV region. With a mass of 10 GeV, the 
77 -width is about 0.1 keV. The authors of [144] proposed to look for such a 
light neutral Higgs boson at the proposed low energy 77 -collider. We want to 
point out that the LHG GaGa heavy ion mode would also be very suitable for 
such a search. A systematic parameter study of the production of Higgs bosons 
in the minimal supersymmetric extension of the standard model (MSSM) was 
done in [145]. For certain values of the free parameters in enhancement of the 
cross section of up to a factor 10 was found. 
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In [76] it was shown how to tag two-photon production in pp collisions at the 
LHC, see Fig. 4 there. It was concluded that “the significant luminosity of 
the tagged two-photon collisions opens an exciting possibility of studying the 
exclusive production of the Higgs boson as well as search for new phenomena.” 

Last year all 4 LEP experiments have reported on their Higgs search in the 
data collected during the LEP-H run. All of them set their new lower limits 
on the Higgs mass which are Mh > 107.0 GeV (L3), 109.7 GeV (OPAL), 

110.6 GeV (ALEPH), 114.3 GeV (DELPHI) at 95% CL [146-149]. However, 
the likelihood analysis shows a preferences for a Higgs boson with a mass of 

115.6 GeV. 

In the Standard Model the arguments of self-consistency of the theory can 
be used to obtain upper and lower limits on the Higgs mass. The upper 
limit is obtained from arguments of perturbative continuation to the GUT 
scale up to Aqut = 10^® GeV [150], i.e., from the requirement that the 
electroweak interaction would be still ’’weak” up to this energy scale. The 
high-energy behavior of the weak interaction was discussed also in a more 
general framework [151,152]. A lower limit for the Higgs mass is obtained 
from quantum corrections involving top-loops to the Higgs interaction po¬ 
tential [153-157]. Thus it appears that the requirement of self-consistency of 
the SM up to 10^® GeV leads to the theoretical limits for the Higgs mass of 
130 GeV < Mho < 190 GeV (to be compared with the lower mass limits of 
Mh > 107.0 — 114.3 GeV at LEP, see above). This range is just the mass 
range acceptable for the Higgs production in two-photon processes of CaCa 
collisions at the LHC. 


4-3 Search for Extra Dimensions 


Graviton production in very peripheral heavy-ion collisions was recently stud¬ 
ied theoretically in [158]. Such a possibility has recently become of interest 
due to a new class of theoretical models, for further explanations see, e.g.. 
Chap. 6 of [107]. In these models quantum gravity and string physics may 
become accessible to experiment; they may even appear directly in the realm 
of LHC and the linear colliders. Graviton production in 77 collisions can be 
calculated in such models. It depends on a fundamental gravity scale M which 
is related to the scale R at which the Newtonian inverse-square law is expected 
to fail (see, e.g., Eq. (25) of [107]). Cross sections for graviton production in 
heavy ion collisions were found to be substantially greater than for graviton 
production in e+e“ collisions [158]. A value of M = ITeV was assumed in 
those calculations. The signature for such a process would be a large missing 
mass in the reaction. In the heavy ion case this is hard to realize experimen¬ 
tally and the authors conclude that “a definite experimental signature cannot 
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be predicted”. For further details we refer to this paper. 


4-4 Search for Heavy Point-like Dirac Monopoles 


In [34] it was proposed to search for heavy magnetic monopoles in 77 collisions 
at hadron colliders like the Tevatron and LHC. The idea is that photon-photon 
scattering 77 —> 77 below the monopole production threshold is enhanced due 
to the strong coupling of magnetic monopoles to photons. The magnetic cou¬ 
pling strength g is given hj g = where n = ±1, ± 2 ,.... Since = 1/137 the 
magnetic coupling strength is indeed quite large. In this reference differential 
cross sections for 77 scattering via the monopole loop are calculated for en¬ 
ergies below the monopole production threshold. The result depends strongly 
on the assumed value of the spin of the monopole. With this elementary cross 
section as an input, the cross section for the process pp 77-I- anything is 
calculated. Elastic, i.e., anything= pp and inelastic contributions are taken 
into account. The signature of such a process is the production of two photons 
where the transverse momentum of the pair is much smaller than the trans¬ 
verse momentum of the individual photons. At the Tevatron such a search was 
performed. They looked at a pair of photons with high transverse energies. No 
excess of events above background was found [159]. Thus a lower limit on 
the mass of the magnetic monopole could be given. A mass of 610, 870, or 
1580 GeV was obtained, for the assumed values of a monopole spin of 0, 1 / 2 , 
or 1 respectively. 


4-5 Tightly Bound States 


One can also speculate about new particles with strong coupling to the 77 - 
channel. As seen in Sec. 3.6 above, large r.y.y-widths will directly lead to large 
77 production cross sections. The two-photon width of quarkonia, for example, 
is proportional to the wave function squared in the center of the system, 
see Table 3. Thus we can expect, that if a system is very tightly bound it 
should have a sufficiently large two-photon width due to the factor |4/(0)p 
which is large in these cases. Examples for such tightly bound systems are 
discussed, e.g., in [160,161]. Composite scalar bosons at ~ 50 GeV are 
expected to have 77 -widths of several MeV [160,161]. The search for such 
kind of resonances in the 77 -production channel will be possible at LHC. In 
Part 3 (p. llOff) of the TESLA Design Report [134] the reader can End some 
interesting remarks about the ’’agnostic” approach to compositeness. From 
Eq. (68 and Fig. 14 one can easily obtain a value for the production cross- 
sections of such states and the corresponding rates in the various collider 
modes. Of course, such ideas are quite speculative. However due to the high 
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flux of equivalent photons such searches seem worth- while, and a possible 
discovery would be quite spectacular. 

Certainly, many aspects in the present section are quite speculative, and one 
will have to wait for the future e’''e“ and 77 colliders to do the physics in the 
region of several 100 GeV. However, it may be possible to have a glimpse into 
this region with the very peripheral collisions at LHC which will be taking 
data in a few years from now. 
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5 Photon-Nucleus and Photon-Hadron Interactions 


Let us start this section with a few rather qualitative remarks which may 
serve as a guideline. Two-photon processes in relativistic ion collisions have a 
relatively simple nature: both ions interact by means of the quasireal photons, 
which they both emit. Theses photons can be thought of as being away, well 
separated from the nuclei (h > Ri + R2 , bi > Ri). The nuclei remain almost 
intact after such a collision. The coherence of the ion interaction is taken into 
account in terms of ion form factors, see Chap. 2. And thus we have virtually 
pure 77-interactions and the rich 77-physics as discussed above. 

Photon-nucleus processes in relativistic ion collisions are in general more com¬ 
plicated. The photon radiated by one (“spectator”) ion interacts with the 
other (“target”) ion in a wide range of photon energies. Besides interacting 
electromagnetically it can also have a hadronic component, due to its fluc¬ 
tuation into a vector meson (“vector meson dominance”, etc.). This leads to 
quite different interaction mechanisms in photon-ion reactions. If the photon 
is of low energy we observe the excitation of the giant dipole or other multi¬ 
pole resonances. As these are relatively low energy excitations, they can occur 
at rather large distances. When the photon (being a non-strongly interacting 
particle) interacts inside the target ion with a single nucleon we observe the 
excitation of nucleon resonances and related phenomena. 

When the photon fluctuates either to a gg-pair or a vector meson before inter¬ 
acting with the nucleons of the target, diffractive processes are possible. These 
can also be called “photon-Pomeron” processes. This signature is similar to 
the photon-photon or Pomeron-Pomeron processes. An important case is the 
diffractive vector meson production. Due to the production with a nuclear tar¬ 
get, both the coherent and incoherent production can take place. As the vector 
meson produced inside the nucleus has to go through the nuclear medium, new 
phenomena will appear. They are widely discussed in the literature in terms of 
color transparency of the ions, of formation- and coherence-lengths of hadrons 
etc., see, e.g., the review [ 162 ]. 

Below we give an overview of all these quite interesting phenomena, we present 
basic formulae and discuss also practical aspects, which are — as it appears 
now — important for experiments at ion colliders. 


5.1 Photonuclear Processes: Giant Dipole Excitation, Beam Loss and Lumi¬ 
nosity Monitor 


At relativistic heavy ion accelerators in the region up to several AGeV the 
study of the electromagnetic excitation of the giant dipole resonance (GDR) 
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and other multipole resonances has been a subject of considerable physics 
interest. Even the excitation of the double phonon giant dipole resonance 
(DGDR) was observed and studied. This is made possible due to the strong 
electromagnetic fields in these collisions, for recent reviews see [163-165]. In 
the hxed target experiments at SPS/CERN with energies in the 100 AGeV 
region, fragmentation processes have been investigated like the hssion of Pb in 
PbPb interactions observed by the NA50 collaboration [166]. Such experiments 
also help to extrapolate to the collider energies (A similar situation is also 
present in the case of bound-free electron-positron production, see Sec. 7.6 
below). While in general a good agreement with theoretical expectations is 
found [167], a puzzling discrepancy still exists. In [168] a measurement of the 
total dissociation cross section of Pb for different targets was done. Beside the 
dominant electromagnetic interaction (proportional to Z^) and the nuclear 
interaction, a component proportional to was found. This contribution is 
in addition to the contribution due to the target electrons, which is sometimes 
called ’’incoherent”. The size of this extra contribution in the case of a Pb 
target is 9.74 barn, larger than the hadronic/nuclear cross section of 7.86 barn. 
This large value precludes explanations in terms of purely nuclear effects. A 
possible interpretation in terms of incoherent electromagnetic processes on 
protons in the target nucleus was investigated in [15]. But the size of this 
effect was found to be too small to explain the results. 

The dominant contribution to the electromagnetic excitation of the ions comes 
from the excitation of the giant dipole resonance (GDR). This is a strongly 
collective nuclear state which exhausts the Thomas-Reiche-Kuhn sum rule 
to a very large extent. The impact parameter dependent probability of the 
one-phonon GDR excitation is given to a good approximation by 

PGDR{b) = S/b\xp[-S/b^) (79) 


for b > Rmin = Ri + R2- The quantity S denotes an area which is given by 

[ 10 ] 


S = 5.45 X lQ-^ZlN2Z2Af^^{m^. 


(80) 


The total cross section is obtained by integrating over the impact parame¬ 
ter from a minimum value of Rmin up to the adiabatic cutoff radius Rmax ~ 
'^ianc/oj. The energy of the GDR is given phenomenologically by hu = 80 MeV/A^P. 
This and the assumption that the classical Thomas-Reiche-Kuhn sum rule is 
exhausted to 100% by the GDR enter into the expression for S given above. 

To lowest order in S/RP^ the cross-section is given by 


(^GDR = 271 Sin • 

\UjlhyYiiffl / 


( 81 ) 


53 



The interest in these (low energy) excitations at the heavy ion colliders is 
mainly motivated by more practical considerations: the huge cross section for 
electromagnetic interaction (of the order of lOO’s of barn) leads to beam loss 
and — as was only noticed very recently [169] — to a local beam pipe heating, 
with the possible danger of quenching of the superconducting magnets. This is 
due to the effect that most of the nuclear excitations (to the GDR) are followed 
by the emission of a few neutrons. The remaining nucleus (with a change in 
the ZjA ratio) will hit the beam pipe and deposit its energy in a rather small 
area. An even more serious problem is the formation of one-electron atoms in 
bound-free electron-positron pair creation (see below, Sec. 7.6). 

On the other hand it was noticed also that the neutrons coming from the GDR 
decay can conveniently be detected in the zero degree calorimeter (ZDG). 
These low energy neutrons in the rest frame of the nucleus are high energy 
neutrons in the collider frame (with energies in the TeV range in the case of 
the LHG). This leads to the possibility of using them as a luminosity monitor. 
Especially the use of mutual GDR excitation has been studied in detail as a 
possible luminosity monitor for RHIG, as well as, LHG [170]. 

The mutual excitation probability is to a good approximation given by the 
product of the probabilities for the single excitations [30]. For symmetric col¬ 
lisions {Zi = Z 21 Ai = A 2 ) we have 

= PoBRl.b) = |r {-2S/A (82) 


and the total cross section is given by the integration over the impact param¬ 
eter from Rmin to 00 : 


ttS 


^mutual 


1 — exp 



TlS^ 

d 2 ’ 
■^min 


(83) 


where the last equation is valid if is small. 

This formalism can be extended to include also the excitation of the double 
giant dipole resonance (DGDR) and even higher phonon states with phonon 
numbers Ni and N 2 respectively. Taking into account the Bose character of 
these excitations one obtains 


PNi,N2iP) 


1 

Ni\N2\ 



N1+N2 

exp 




(84) 


The total cross section is found again by integration over the impact parameter 
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from Rmin to cx): 


<^Ni,N2 — 


{Ni +N2-2 )\tiS 

iVi!iV2!2^i+^2-i 

nS 


1 - exp(-2^/i?^i„) 


S 


A^iH-A^2 —1 


A^i+A^ 2~2 

E 


fc =0 



Ni\N2\{N^ + N2-l) 


p 2 

-^min 




where again the last equation is valid to lowest order in Even for the 

heaviest systems is less than one, but not very much: with Rmin = 

2Ra) we obtain S/R^^^ = 0.53 and 0.48 for the system PbPb {A = 208) 
and AuAu {A = 197), respectively. The probability to excite a one phonon 
GDR in close collisions is thus of the order of 30%. Detailed calculations using 
photonuclear data (instead of the gross properties of the GDR as discussed 
above) as an input were performed by Vidovic et ah [171] and Baltz et al. 
[172]. An improved value of A = (17.4 fm)^, including photonuclear excitation 
beyond the GDR, was given in [173] for the PbPb collisions at LHG, based on 
1172], 

Glose collisions of ions contribute most to the mutual excitation. This leads to 
a stronger dependence of the mutual excitation cross section on the detailed 
form of the cutoff, as compared to the single neutron production, with its 
rather weak logarithmic dependence on the cutoff radius. On the other hand 
the measurement of the mutual excitation process is less affected by possi¬ 
ble background processes — like beam-gas interactions — than the single 
GDR excitation. This makes it more suitable as a luminosity monitor. Mutual 
excitation is sensitive to some extent to nuclear effects. The theoretical calcu¬ 
lations [170] and especially [174,167] show that nevertheless a good accuracy 
can be achieved. In [167] it is concluded that “Good description of CERN SPS 
experimental data on Au and Pb dissociation gives confidence in the predictive 
power of the model for AuAu and PbPb collisions at RHIC and LHC”. This 
method to monitor the luminosity is currently used at RHIG and plans are 
also underway to use it at LHG/GERN. In [175] it was concluded that this 
method is possible also for the ALIGE experiment, see also the discussion in 
Ghap. 8 below. 

Very recently the hrst experimental results have become available from the 
ZDG at RHIG [176]. The ratio of the total cross section (mutual Goulomb plus 
nuclear interaction) and nuclear interaction alone was measured to be 0.64 in 
good agreement with the expected theoretical result of 0.66. Also ratios of 
cross sections for the emission of different numbers of neutrons from each ion 
were compared with theoretical predictions. The energy spectrum in the ZDG 
measured in the experiment at RHIG is shown in Fig. 19, taken from [177], 
where the dominant peak is caused by high-energy neutrons in the ZDG. 
At LHG energies a clear structure due to multiple-neutron emission of the 
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Fig. 19. Online ZDC energy distribution of neutrons obtained during RHIC colliding 
beam operation with beam energy 65 GeV per nucleon. Extracted from [177]. 

ions is predicted, see Fig. 20 [178], where the expected energy distribntion for 
nentrons in the ZDC for the ALICE experiment is shown. Another application 



Fig. 20. Expected ZDC energy distribution of neutrons from e.m. of giant dipole 
resonance in the ALICE experiment at LHC. The ZDC energy resolution is included 
in the calculations. Extracted from [178]. 

of the neutron emission processes in heavy ion collisions is to use them as a 
trigger for the electromagnetic interaction (only) of the ions in the collisions 
by observing neutrons in both of the ZDC’s. This opens the possibility to 
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select 77 -processes and suppress other peripheral processes in ion collisions 
— like photon-Pomeron and Pomeron-Pomeron processes in experiments with 
an open geometry, i.e., where it is impossible to identify the 77 -processes by 
measuring directly the very small Pt of the system produced in 77 -processes. 
This will be discussed in more detail in Chap. 8 below. 

For both of these purposes the cross section needs to be known very well. Espe¬ 
cially the different decay channels resulting in a different number of neutrons 
emitted have therefore been studied in a detailed investigation in [174,167]. 
For the simulation the internuclear cascade in the excited nucleus was used. 


5.2 Excitation of Nucleon Resonances 


The electromagnetic excitation of the A resonance in peripheral heavy ion 
collisions was already observed at the SPS [179]. Cross sections for the elec¬ 
tromagnetic excitation of the nucleon resonances are sizeable at the collider 
energies. This was hrst studied in [41], where the following values are given: 
17 b for UU collisions at RHIC and 25 b for PbPb collisions at LHC. These 
numbers are based on the folding of the equivalent photon spectra with an 
experimental input for the 7 -nucleon cross-sections up to E.y = 2 GeV (in the 
nucleon rest frame). For further details we refer to that reference. These reso¬ 
nances predominantly decay by the emission of mesons, mainly pions. These 
pions will be emitted in a rather narrow cone in the beam direction. Maybe 
they are useful? The cross section for nucleon resonances is 25 b for PbPb 
collisions at LHC, compared to 8 b for the total nuclear cross section. With 
the branching ratio of about 33% for the decay of the A to mr~^ and of about 
66 % for the decay to pn^, and relatively small decay energy, a beam of pions 
with a narrow energy range will be produced. For a recent reference of pho¬ 
toproduction of mesons we refer to [180]. Some time ago the production of 
beams of pions and other strongly interacting particles by photons was con¬ 
sidered by Drell in [181]. A detailed study of the pion emission in heavy ion 
collisions is performed in [182]. Various mechanisms are considered along with 
their characteristics. The electromagnetic excitation of nucleon resonances is 
found to be one of the important channels. 


5.3 Photonuclear Reactions to the Continuum above the Nucleon Resonances 


The continuum above the nucleon resonances is also strongly excited at collider 
energies. A hrst estimate of this effect was given in [41]. In this exploratory 
calculation a constant value (independent of the photon energy) of a = 100 /ib 
was assumed for photon energies above 2 GeV (in the nucleon rest frame). This 
gives a value of 2.6 b for UU at RHIG and 19 b for PbPb at LHG. More detailed 
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calculations were done in [43] and [171]. In this last work experimental values 
of photonuclear cross-sections for Pb up to 80 GeV were used as an input into 
the calculations. For PbPb collisions at LHC a cross section of 53.7 b was 
given for the region between 40 MeV and 2 GeV. From 2 GeV to 80 GeV the 
cross section is 18.7 b in good agreement with the estimate above. For AuAu 
collisions at RHIC the corresponding experimental photonuclear data points 
only exist up to 9.5 GeV. Scaling the experimentally known cross section for 
Pb to the case of Au a value of 5.6 b was given for the region of 2-80 GeV. For 
the region above 80 GeV, where no experimental data were available, several 
parameterizations have been used. Only a small cross section of 0.7 barn was 
found for AuAu at RHIC, in contrast to 16.3 barn for PbPb collisions at LHC. 
This is due to the much harder equivalent photon spectrum at the higher beam 
energies at the LHC. 

As has been mentioned already before, the electromagnetic breakup of the 
nuclei is one of the main loss processes in PbPb collisions at LHC. Therefore 
a precise knowledge of these cross sections also at higher photon energies is of 
interest. 


5-4 Coherent Vector Meson Production in Photon Diffraction 


A lot of theoretical studies of coherent, as well as, incoherent photoproduction 
of vector mesons in nuclei exist. With the hard equivalent photon spectra 
present at the relativistic heavy ion colliders new experimental perspectives are 
opened up for these studies. As was emphasized above the photon virtuality 
is essentially restricted to zero, as opposed to the situation in electron 
scattering, where this quantity can be varied. Let us mention here some basic 
concepts of this vast held, where we essentially stick to the case of quasireal 
photons. 

Among the many possible channels, the coherent production of vector mesons 
(that is, when the nucleus stays in the ground state), as well as, other G-odd 
mesons in photon-Pomeron processes is a very interesting one to be studied 
in peripheral heavy ion collisions. From an experimental point of view these 
processes are quite similar to the two-photon processes discussed above. The 
main difference lies in the quantum number of the produced system, which 
in a photon-Pomeron collisions corresponds to C = —1. The low value of the 
transverse momentum of the mesons is a very clear experimental signature 
for this, comparable to the similar one in the photon-photon case. In fact 
the production of mesons was already observed at STAR/RHIG [50] by 
looking for its decay into 7r’''7r“ pairs, see Sec. 8.3.2 below. The invariant mass 
distribution for pairs with P^ < 0.1 GeV shows a clear peak at the p mass. 
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Let us describe the main ideas to treat coherent vector meson production, 
where we follow essentially the work of [183] and [19]. The calculation is done 
in the rest frame of one of the nuclei. The minimum momentum transfer 
is given to a good accuracy by tmin = {My/2kY, where My is the vector 
meson mass and k the photon momentum. This quantity (it corresponds to 
the minimal q\ dehned in Eq. (92) below) is very small (essentially zero) for 
high photon momenta. The coherence length lc = I/q'l is larger than the size 
of the nucleus. This means that the photon will fluctuate into a vector meson 
pair already outside the nucleus and this virtual meson is then put on-shell in 
the interaction with the nucleus. 


We are interested in the total cross section for the reaction yA — VA. To 
get to this several steps have to be done. One starts with the experimentally 
known cross section for jp ^ Vp as an input. This cross section can usually 
be parameterized in the Regge form [184] 


daipfp Vp) 


dt 


i=0 


bv (X\V‘ + YW-"') , 


(85) 


where W is the center-of-mass energy and the constants by, X, Y, e and p are 
determined from hts to the data (compare also with the parameterization of 
the total 77 cross section Eq. (57)). A major simplihcation comes from the use 
of vector meson dominance [185], which allows to relate this photoproduction 
cross section to the cross-section for forward elastic Vp Vp scattering: 


dai^pp Vp) 


dt 


t=Q 


2na d(j{Vp Vp) 
If dt 




t=0 


( 86 ) 


where fv denotes the vector-meson-photon coupling. The further calculations 
are simplified if one assumes that Re/(0)/Im/(0) 1, i.e., the Vp scattering 

is mainly absorptive. This assumption is, e.g., relaxed in [186] at the expense 
of having to introduce a new parameter. This assumption is certainly better 
fulhlled for the p meson than for the J/T. With this the forward scattering 
cross section is proportional to the square of the imaginary part of the forward 
scattering amplitude. The optical theorem can be used to relate this forward 
elastic scattering cross section to the total cross section avp for vector meson 
nucleon scattering: 



IGtt 


daiVp Vp) 
dt 


t=o 


(87) 


Using this elementary cross section ayp and the nuclear density distribution 
p^(r) we can obtain in a Glauber calculation [183] the total elastic scattering 
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cross section as well as the total cross section for VA scattering: 


o-eiastic(V^A) = J d% [1 - exp {-avpTAib)/2)f 
o-totai(V^A) = 2 J d%[l - exp {-avpTA{b)/2 )], 


(88) 


where Ta is the nuclear thickness function 


= j dzpAVb^ + z^-). 


(89) 


VMD allows us now to relate the cross section for elastic VA scattering back 
to the coherent production cross section 

OiTTCy. 

a{^A ^ VA) = ^^(Telastic(V^A). (90) 

Jv 


This expression allows us to discuss the dependence of the cross section on A 
in two limiting cases. In the transparent limit, we expand the exponent to first 
order. Using a nuclear density distribution corresponding to a homogeneous 
charged sphere, integration of Eq. ( 88 ) gives us < 7 ( 7 ^ VA) AA3. In the 
black disc limit on the other hand we hnd a( 7 A —> UA) ^ 712/3. If only 

a hnite range of t needs to be taken into account, the differential scattering 
cross section da/dt, see [183] can be used. 


A slightly different path is used in [19]. The optical theorem for nucleus A and 
vector meson dominance is used to get 

d(j{'yA UA) 
dt 


aa. 


total 


(UA) 


t=n 


4/. 


(91) 


V 


From this equation one can directly see the A-dependence of the coherent for¬ 
ward photoproduction amplitude for the two limiting cases: in the transparent 
limit there is a A^ behavior (typical for coherent processes) and in the black 
disc limit we have an AA3 rise with nucleon number A. 

The t-dependence of the coherent photoproduction is modeled by the nuclear 
form factor Fei{t). A convenient parameterization of this form factor is given 
in Eq. (14) of [19]. This form factor falls off very fast with f, the angular distri¬ 
bution is peaked very much in the forward direction (up to a tmax of the order 
of 1/i?^). Integration over t in order to obtain the total coherent cross section 
yields a factor of A“2/3 Thus the total coherent photoproduction scales as 
A 4/3 transparent and as A2/3 for the black nucleus case, in agreement 

with the discussion above. The above description illustrates an important idea. 
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which goes back to the work of Drell and Trefil [187]: using the nucleus as a 
laboratory, one can determine otherwise unaccessible vector-meson-nucleon 
scattering cross sections. Several groups studied such processes already in the 
1960s and 1970s. This is reviewed in [185]. More recent experiments using 
virtual photons from muon or electron scattering are less comparable because 
the photons are more virtual (compared to 1/Ra)- There are some experi¬ 
mental results on coherent cross sections: older data exist from DESY and 
Cornell. In Figs. 5,6,7 and 9 of [186] there are good fits to the Cornell data. 
Photons of = 6.1 and 8.8 GeV were scattered on complex nuclei like C, 
Mg, Cu, Ag, Pb and U. At RHIC and LHC one has the possibility to study 
these processes in great detail, for vector mesons up to the Upsilon and over 
a much wider energy range. With the forthcoming new detailed results, the 
theoretical analysis can be improved, for a detailed formalism see, e.g., [183] 
or [186]. 

Further interesting ideas are advanced in [19,20,188]. Since either nucleus can 
emit the photon, there is an interference between these two indistinguishable 
processes. Such interference effects are discussed theoretically in [20,188]. We 
refer the reader to this work for further details. Interference effects are sensitive 
to amplitudes. Therefore it will be interesting to see whether such effects can 
tell us more about the real and imaginary parts of the amplitudes involved. 

Another possibility is multiple vector meson production. There is an analogy 
to the double phonon giant dipole excitation, see, e.g., [10,163]: in both cases 
identical bosons are produced and symmetrisation has to be done. This is very 
important in the case of the double phonon states, since the dipole phonons 
are only labeled by their spin-projection quantum numbers M = 0, ±1. On the 
other hand for the produced vector mesons, there is in addition a continuous 
quantum number — namely their momenta — which can be quite different 
in general. Therefore one can expect that boson effects are less important in 
this case. However, it would be most interesting to observe the production of 
two identical vector mesons with close enough momenta and study such Boson 
effects. For details and further references, see [20]. 

It is noted in [20] and [188] that coherent vector meson production cross 
sections are huge. E.g., the p production in AuAu collisions at RHIC is 10 % 
of the total hadronic cross-section, in PbPb collisions at LHC it is about 
equal to the total hadronic cross-section. So heavy ion colliders can act as 
vector meson factories with rates comparable to those of e’''e“ vector meson 
machines. Up to 10^° $ mesons will be produced in 10® s with Ca beams at 
LHC. This can be compared to the expected rates at the dedicated <l>-factory 
DA<hNE [189]: at the beginning, a number of 2.2 x 10^® yearly produced $ is 
expected. This will allow for searches of rare decay modes, CP violation and 
the like. Also there is the possibility to do vector meson spectroscopy: mesons 
like the p(1450), p(1700) and $(1680) will be copiously produced. The current 
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situation of vector meson spectroscopy is outlined in [190]. It is concluded there 
that many interesting questions are still not understood. Hopefully relativistic 
heavy ion colliders can contribute to these QCD studies. Of course, it will have 
to be seen how well such experiments can be done in the hostile environment 
of the violent central collisions with their extremely high multiplicities. 


5.5 Incoherent Production of Vector Mesons 


In the incoherent vector meson production process the photon interacts with 
a nucleon to produce a vector meson leaving the nucleus in an excited state. 
We follow here mainly the formalism of [191]. 

The coherence length is dehned as /c = I/ql, where qi is the difference between 
the longitudinal momenta of the photon and the produced vector meson. It is 
given by: 


My + Pt 


(92) 


where pt denotes the transverse momentum of the produced vector meson 
with mass My, and uj is the photon energy. The incoherent cross section for 
the production of a vector meson in a reaction, where the nuclear state changes 
from |0) to I/) (/ ^ 0) in Glauber theory is given by: 



where the transition operator (b) is dehned as 


ry(6; {?,,*-,}) = 




k(M) 


1 - - Sk) &{zk - Zj) 


(94) 


It consists of the vector-meson production amplitude and an amplitude 
, which describes the elastic scattering of the vector meson on the nucleon. 
By taking |/) = |0), we can treat also the elastic/coherent case, above. Of 
course the equations there are related to the one here by making use of the 
“cumulant expansion” [192]. 

Because of the phase factor in this equation, the photon production 

amplitude on two nucleons with positions \zi — Zf \ < Ic adds up coherently. 
E.g., with My = 0.8 GeV and u = 8 GeV we have /c = 4 fm (for pr ~ 0). 
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The equivalent photon spectrum at RHIC extends up to several hundred GeV 
(as viewed from an orbiting nucleus), therefore this coherence condition will 
be fulhlled in many cases. 

The low and high energy limits of incoherent photoproduction are given in 
Eq. (12) of [191] as 


aJnc = -^VN) j (fb J dz p{h,z) x 

—oo 

f (low energy) 1 

I (high energy) I 


(95) 


This formula contains the elementary vector meson production cross section 
and a shadow correction factor. In the low energy limit, with Ic -C R, the atten¬ 
uation of the outgoing vector meson is governed by the thickness Tz{b), which 
is the thickness experienced by the vector meson from its point of creation 
{b, z) till its exit from the nucleus. In the high energy limit the attenuation 
is governed by T{b) = Tz^-ooib), the thickness of the nucleus along the to¬ 
tal path. This corresponds to the interpretation already given above, that at 
high energies the photon has already converted to a virtual vector meson long 
before hitting the nucleus. The meson is then put on the mass shell by the in¬ 
teraction with a target nucleon. This interpretation is familiar from the vector 
dominance model. 

An interesting quantity is the nuclear transparency. It is dehned as the ratio 
of the incoherent prodnction cross section on a nncleus A to A times the 
elementary photoprodnction cross section on a nncleon. In Fig. 21 we show the 
cross sections for incoherent vector meson prodnction at the LHC. A nnclear 
transparency of 1 was assnmed. This should (of course) be scaled with some 
realistic value of the transparency factor, which can, e.g., be read off from 
Fig. 1 of [191]. For a Pb nncleus those authors find a valne of abont 0.1 for 
cu > 10 GeV (g2 = 0). 

Both the coherent and incoherent prodnction of vector mesons have been 
stndied at HERMES/HERA in eA collisions [193]. The photon virtnality was 
0.4 GeV^ < < 5 GeV^which is mnch larger than the < 1/R^ in 

AA collisions. A theoretical discussion is given in [194] where a systematic 
mnltiple scattering formalism is developed for vector meson electroprodnction 
from nnclei. Also formation-, propagation-, and hadronization-scales for qnark- 
glnon flnctnations of the virtnal high-energy photon are considered. 
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Fig. 21. The cross sections for the diffractive vector meson production of p,uj,J/'^, 
and T, as well as, the total diffractive cross section is shown for the collision of a 
virtual photon coming from a Pb ion at LHC energies with a proton as a function of 
the invariant mass M of the photon-proton system. For the completely incoherent 
cross section, this cross section needs to be multiplied by A. Parameter for the 
different cross section are taken from [184]. 

5.6 Photon-Gluon Fusion and the Gluon Structure Function in a Nuclear 
Environment 


The production of cc and bb pairs in photon-gluon processes has been sug¬ 
gested as a method to measure the gluon distribution in a nucleon by Frixione 
et ah [195]. Such processes are under study experimentally at HERA [57] (es¬ 
pecially [196]), and [58] (especially [197]). Also there has been much progress 
in the theoretical study of these processes and the investigation of higher order 
effects. A recent reference is [198]. With the high photon fluxes at relativis¬ 
tic heavy ion colliders such studies will also be possible for nuclei instead of 
nucleons. This is reviewed in [16] and we give only a brief update here. 

In these photon-gluon fusion processes, one of the ions provides the photon, 
which hits the gluon of the other ion to produce a qq pair. Such processes are 
coherent for one ion and incoherent for the other. This leads to a special kind of 
event topology: there is one nucleus, which remains intact (or at least almost, 
the “spectator”, which provides the equivalent photons), and the “target” 
(which provides the gluon), will be broken up (often called the “remnant”) 
plus a qq jet. 

The first estimate of the production of heavy quarks in such collisions was 
done in [199]. Central collisions were not excluded in their approach. More 
refined calculations focusing on the production of cc and bb quarks were done 
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in [42-44,16]. 


The flavour content of the quark jets can be determined by detecting the 
corresponding D or B mesons. This of course needs to be included in the 
specification of the triggering. In this context the study of inclusive processes 
with heavy flavour production is also of interest. We mention the measurement 
of inclusive Df photoproduction at HERA [200]. The photon-proton c.m. en¬ 
ergies were in the range of 130 < W < 280 GeV. Photoproduction events 
were selected by the requirement that no scattered positron was identified. 
The c.m. energies of the photon-proton system can therefore not be deter¬ 
mined directly from the momentum of the scattered positron, but needs to be 
extracted from the transverse momentum of the jet from the two produced 
quarks. They were determined with the Jacquet-Blondel estimator of W, see 
also Sec. 7.4 and especially Eq. (116). 

Modification of parton distribution functions inside the nucleus is one of the 
interesting topics of QCD [201-203]. These modifications are expected to be 
present especially at small value of x, where the partons distributions of the 
individual nucleons start to overlap with each other. In a recent work [204] 
the gluon distribution in this low-a; region was assumed to be given by the 
“color glass condensate”. In this model, the photoproduction of a qq pair 
in very peripheral AA collisions was calculated. The authors conclude that 
“R appears that this process is very sensitive to the properties of the gluon 
distribution in the region where saturation might play an important role". 

Quite recently the production of top quarks at the LHC was studied in [45]. It 
was found there that 210 ti pairs will be produced in a 10®s 00 run at LHC 
(with an assumed luminosity of 1.4 x 10^^cm“^s“^). Top quark production has 
been studied before in [199]. In [45] the photon spectrum excluding central 
collisions and a more modern version of the gluon distribution function is 
used. The results are sensitive to the large behavior of the distribution 
function. Shadowing is also included. Eor more details we refer to this work. 
It will be an interesting channel to be studied at the forthcoming LHC. 
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6 Diffractive Processes 


Diffractive processes are of interest in pp and AA collisions for their own sake 
and as a possible background to 77 collisions. At HERA diffractive photon- 
proton interactions have been studied in great detail, see [205-207]. The pro¬ 
cesses of diffractive photoproduction of vector mesons on nuclei were discussed 
in subsections 5.4 and 3.5. Elastic proton proton scattering will be studied at 
RHIC [208]. 

In this subsection we discuss double diffractive (or double Pomeron) processes 
in ion-ion collisions. This case is certainly even more difficult than the diffrac¬ 
tion in the 77 and 7 p interactions, as discussed above. The double diffractive 
processes have been studied already at the ISR in pp-collisions and at the 
TEVATRON in pp collisions [209-211]. These processes are also one of the 
main topic in the current COMPASS experiment at the SPS [ 212 ]. 

Generally one could dehne diffractive processes as hadronic or ion interactions 
by means of exchange of one or several Pomerons. The Pomeron itself is re¬ 
garded in QCD as a colorless object consisting of 2 (correlated) gluons gg. 
Also discussed now are three-gluon objects ggg, i.e., Odderons. The simplest 
diffractive process meeting this dehnition is the elastic scattering and there 
are several special experimental programs devoted to this question in pp colli¬ 
sions [213,214,36]. For the elastic scattering in AA collisions, see also Chap. 2 
in [15]. 

We are interested in diffraction mainly as a possible background for photon- 
photon and photon-hadron processes. We will therefore look only at those 
processes where particles are produced in the central region together with 
no breakup of the two ions, as they correspond to the tagging conditions 
for photon-photon and photon-hadron processes, see Chap. 8 . Experimentally 
these double diffractive processes (“double Pomeron exchange”) have been 
studied for alpha-alpha collisions already at the ISR in 1985 [215]. Results for 
double diffraction in pp collisions at the TEVATRON have also been studied 
recently [210]. Diffractive photon-nucleus processes could also be studied with 
nuclear targets in a future option at HERA [136]. For an overview of diffractive 
processes we refer to [216,217,24]. 

These double diffractive processes are characterized by hadrons produced in 
a region of rapidity y well separated by so-called “rapidity gaps” from the 
rapidity of the initial ions, see Fig. 22. If a colored object (e.g., a single quark 
or gluon) would be transfered from one of the ions, the probability for the 
existence of such a gap is suppressed exponentially with y [216,36], therefore 
these events will be mainly due to the exchange of colorless objects. This 
characteristics coincides also with the condition for Ending the two ions in 
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the ground state (“coherent case”). This process will therefore be due to the 
exchange of either neutral mesons or Pomerons, where these two contributions 
were already present in the total cross section 77 —hadrons, discussed in 
Chap. 3 above. At high energies, as it takes place in relativistic ion collisions, 
only the Pomeron exchanges are dominant in the cross sections. 

rapidity gap central region 


• • 


y 

Fig. 22. Coherent double diffractive processes are characterized by particles pro¬ 
duced in the central rapidity region and the two ions leaving the interaction with 
almost the initial rapidity. Between them is a large empty region, defining the “ra¬ 
pidity gap”. Adapted from [216] 


In the phenomenological approach of Donnachie and Landshoff [218-221] the 
Pomeron is described as a single Regge pole with a “propagator” 


Dp{t,s) = - j- -TTV^xp 

sm (2 Trap 


-7rap(t) 


(96) 


where s is the square of the total c.m. energy, t = —q^ the square of the 
momentum transfer and the Regge trajectory of the Pomeron given by 

*^f(^) = 1 + e -|- (Upf, (97) 


with e = 0.085 and a'p = 0.25 GeV ^ (see also the parameterization mentioned 
above for 77 —y hadron). The coupling of the Pomeron to the proton is given 
by 




(98) 


with /do = 1.8 GeV“^ and Fi is taken as the isoscalar electromagnetic form fac¬ 
tor, for which the usual dipole parameterization is assumed. How the Pomeron 
couples with the whole nucleus is an interesting question. Very little seems to 
be known about this at present. This will introduce some uncertainty into the 
calculations discussed here. In general it is assumed to be of the form 

^Apit) = ^^PoFAi—t), (99) 


that is, a coherent coupling to all nucleons ~ A with the elastic form factor 
guaranteeing that the Pomeron emission does not lead to a breakup of the 
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ion. The coherence will be destroyed rather quickly with increasing t, t being 
restricted essentially by \t\ ^ ~ (60 MeV)^, compared to the “natural 

range” given by Regge theory |t| ^ 1/rg ^ (500 MeV)^ with Tq = a'p ln(s/m^). 

The use of the elastic form factor corresponds to taking into account the hnite 
size of the nucleus, see the discussion in Chap. 2 . The additional hadronic “ini¬ 
tial state interaction” between the other nuclei need to be taken into account. 
Whereas their effect was small in most cases in the photon-photon collisions, it 
will be much more severe here, due to the short range nature of the Pomeron. 
Also the “hnal state interaction”, the third effect discussed in Chap. 2 , will 
be more important in this case. Whereas the electromagnetic held decreases 
inside the nucleus and therefore the number of equivalent photons there is 
small, this is not the case for the equivalent Pomerons inside the nucleus. 
Both the hnite size and the initial state interaction ehects are assumed to lead 
to a considerable decrease of the (ehective) coupling so that we hnally expect 
only an increase of the cross section with where <5 < 1 . 

This model for the Pomeron has been used as the basis in a number of studies 
within the “equivalent Pomeron approximation” [222]. The inclusive produc¬ 
tion of the Higgs boson (PP —TT-I-anything) was studied for pp collisions 
in [223], see Fig. 23(a). It was assumed that the major contribution for the 
Higgs boson production comes from the gluon-gluon fusion via a heavy-quark 
loop into the Higgs. The gluon distribution function inside the Pomeron is not 
fully settled, but it was assumed that the Pomeron consists mainly of gluons. 
A gluon distribution of the form 

G^{x) = Q/x{l — x)^ (100) 


was used. 

This approach was extended to AA collisions in [ 222 ]. Also an estimate for 
the exclusive Higgs boson production (PP — H, where the Higgs boson alone 
without any other hadronic particles is produced, see Fig. 23(b)) was given. 
The Pomeron-Pomeron-Higgs coupling was modeled by assuming that the 
Pomeron is a correlated two gluon exchange [218], and the coupling to the 
Higgs boson was then done again via a heavy quark loop, see Fig. 23(b). 
In contrast to the photon the Pomeron, being not a real particle but a phe¬ 
nomenological description, has a hnite size. Therefore the exclusive production 
of some heavy mass resonances (corresponding to an object with a small size) 
is highly suppressed. 

Central collisions were excluded in [ 222 ] by using only the impact parameter 
range b > Ri+R 2 . This exclusion of the central collisions leads to a suppression 
of the inclusive production cross section by six orders of magnitude. It was 
concluded that both the exclusive and inclusive production are small compared 



to the photon-photon case. The exclusive production of rj mesons was studied 
similarly in [224], Again it was found to be only a small background compared 
to the dominant photon-photon process. 



(b) 


Pomeron-Pomeron interaction. The dominant process in both cases comes from the 
coupling of the Higgs boson with the two Pomerons or gluons through a heavy quark 
loop. 


The calculation of the production of a central cluster of hadrons with a given 
invariant mass was done in [21] within the Dual Parton Model [225]. By as¬ 
suming that only one of the nucleons produces a central cluster and all other 
nucleons scatter only elastically, central collisions were removed. Within a 
Glauber-like calculation of this, the dependence of the cross section on A is 
given by a ~ This can be understood by assuming that there is only 
a “ring zone” surrounding each ion, which is active for Pomeron-Pomeron in¬ 
teractions, see Fig. 24, corresponding to the black disc limit. It was concluded 
that the cross-section for Photon-Pomeron and Pomeron-Pomeron interac¬ 
tions for this case are larger than the cross section for photon-photon events 
for almost all ion species except for the very heavy ones (Pb), were they are 
comparable. 



Fig. 24. Due to the strong interaction of the nucleons in both ions and the short 
range of the Pomeron, essentially only a narrow ring around the nuclear radius, 
contributes to the diffractive interactions without nuclear breakup. As the nuclear 
radius increases with R ~ A^A this explains the increase of the Pomeron-Pomeron 
cross section with A^A_ 

Whereas the central collisions were effectively removed in [21], the coherence 
condition, which essentially guarantees that the ion does not break up due 
to the Pomeron emission, was not included in this approach. As mentioned 
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above the Pomeron is a short ranged object, it would normally give a large 
momentum transfer of the order of 500 MeV to the nucleus. Therefore it is 
very likely that the nucleus will break up. An elastic form factor for the whole 
ion was not taken into account in this approach. Some estimates were made 
[226], indicating that including this would lead to a further suppression of the 
diffractive events, favoring again the photon-photon case. 

This further decrease of the diffractive production cross section due to co¬ 
herence effects was confirmed also in [227] within the equivalent Pomeron 
approach, see above. Three different types of hnal states were considered and 
compared to the photon-photon case: the exclusive production of resonances 
(mainly different rj mesons) was studied, as well as, the (continuum) produc¬ 
tion of a TT pair. The contribution of Pomeron-Pomeron interaction compared 
to photon-photon collisions in these cases is found to be only 1% or less for 
both RHIC and LHC for all ions. This is again due to the fact that the exclu¬ 
sive process is highly suppressed. As a third type of reaction the production 
of a central hadronic cluster is studied. Here it is found that for light ions 
(Ca,Ag) at LHC diffractive processes can become of the same size and even 
dominate at an invariant mass of the cluster below 5 GeV. Also at RHIC 
the diffractive processes will dominate over the photon-photon processes. No 
results are given for the photon-Pomeron interactions in [227]. 

A useful expression for the production of a central cluster of hadrons in 
Pomeron-Pomeron processes is given in the framework of the simple Regge 
model, see, i.e., [228,24]. The cross section of these processes are given ap¬ 
proximately by 


dapp/dm\ = + 5) In 


snn 


m 
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{mxmNRy 


( 101 ) 


where do = 1.3 j^b, mx is the mass of the produced hadron cluster in a central 
rapidity region, satat is the square of the c.m. energy of the nucleon-nucleon 
collisions, R the radius of the nucleus, nix is the nucleon mass and 6 is the 
width of the nuclear surface. After integration of Eq. (101) over dm'j^ an es¬ 
timate of the cross section for Pomeron-Pomeron processes can be obtained. 
One sees again that this cross section roughly increases with that is, 

also with This should be compared with the factor in the cross sec¬ 
tion for two-photon processes. Thus in heavy ion collisions we can expect the 
predominance of two-photon processes, whereas in light ion- and especially 
in pp-collisions the contribution of Pomeron-Pomeron processes will be much 
more essential. These more qualitative cross section estimates are fully con- 
hrmed by detailed calculations, see above and Sec. 8.2.1. 

Due to strong shadowing of the Pomeron inside the nucleus only the surface 
region will contribute to peripheral diffractive processes. Photon-photon and 
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Pomeron-Pomeron processes will therefore scale approximately with the ion 
species as and respectively [24], Thus for heavy ions, like Pb, we may 
expect dominance of the photon-photon processes whereas, say in pp collisions, 
the Pomeron-Pomeron processes will dominate in coherent collisions. 

Diffractive processes are an additional held of studies in peripheral ion col¬ 
lisions. They essentially have the same triggering conditions and therefore 
one should be able to record them at the same time as photon-photon and 
photon-hadron events. One may think also about going beyond the coherent 
case (“no breakup”). Processes with a large rapidity gap have been proposed 
[229,230] as an interesting physics topic and many interesting questions, like 
the probability for a gap as a function of its width y have been posed. 

As a potential background for coherent photon-photon and photon-hadron 
collisions in very peripheral collisions we conclude that they are in most cases 
only a small background for the heaviest ions. Without the triggering on “no 
breakup” for the two ions, they dominate the production of hadrons in the 
central region. The same is also the case for pp collisions at the LHC. 
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7 Special Aspects of Dilepton Pair Production 


Due to their small mass, electrons (positrons) and to some extent also muons 
play a special role in very peripheral heavy ion collisions. They are produced 
much more easily than other heavier particles, and — especially in the case 
of e'''e“ pair production — some new phenomena, like multiple pair produc¬ 
tion, will occur. A small mass means that the Compton wave length — for 
electrons this is Ac = — ~ 386 fm — is large, larger than the nuclear ra¬ 
dius R of several fm. The equivalent photon approximation has to be modihed 
when applied in this case. The cutoff radius is not given by the nuclear radius 
anymore but by the Compton wavelength of the particle. The reason is that 
in the cross section for the 77 —>■ e’''e“ subprocess the dependence on the 
virtualities Q\ and Q 2 of the (quasireal) photons can no longer be neglected. 
As was explained in Chapter 2 the photon virtuality in heavy ion collisions is 
limited by the nuclear radius up to ^ R^i which is much larger than 
Especially for invariant masses of the order of several rrie — which contribute 
mostly to the total cross section — the photon-photon cross section deviates 
from the one for real photons for virtualities Qf and/or Q 2 ^ m^, (see App. E 
and Eq. (6.25) of [59]). For the muon — with a Compton wavelength of about 
2 fm — we expect the standard equivalent photon approximation (EPA) to 
be applicable with smaller corrections. For pair production using EPA 
together with the cutoff radius equal to the electron Compton wave length 
gives a fair approximation for the total cross section. On the other hand EPA 
completely fails for certain regions of the phase space and also for the calcula¬ 
tion of pair production probabilities at small impact parameters (as discussed 
in detail below). 

The calculation of pair production in heavy ion collisions dates back to the 
work of Landau & Lifschitz [25] and Racah [26] in the thirties, as already 
mentioned in the introduction. In connection with the relativistic heavy ion 
colliders there has been some renewed interest and cross section calculations 
have been done using the semiclassical approximation in [10] and using Feyn¬ 
man Monte Carlo techniques in [231]. For lower energies (about 1-2 GeV per 
nucleon) coupled channel calculations have been done also by a number of 
groups [232,233]. 

The total cross section for pair production is huge (about 200 kbarn for PbPb 
at LHC, 30 kbarn for AuAu at RHIC [234]). Electron-positron pairs therefore 
present a possible background. The differential cross section show that most 
of the particles are produced at low invariant masses (below 10 MeV) and into 
the very forward direction, see Fig. 25. The highly energetic electrons and po¬ 
sitrons tend to be even more concentrated along the beam pipe, therefore most 
of them remain unobserved. On the other hand, the cross section for highly 
energetic electrons and positrons is still of the order of barns. These QED 
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pairs constitute a potential hazard for the detectors, see below in Chap. 8. On 
the other hand they can also be useful as a possible luminosity monitor, as is 
discussed, e.g., in [24,35]. 




e(deg) (b) 


Fig. 25. Single differential cross sections for the e'’'e“ pair production at LHC as a 
function of the energy (a) of either electron or positron and as a function of the angle 
of the electron or positron with the beam axis (b). Most pairs are produced with 
energies between 2-5 MeV in the Lab frame and in the very forward or backward 
direction. 


Due to the strong electromagnetic helds of short duration new phenomena, 
especially multiple pair production, but also Coulomb-correct ions, will appear 
in relativistic heavy ion collisions. They are of interest for the study of QED of 
strong helds. This will be discussed in the next two sections. In addition there 
are also a number of processes related to pair production, like bremsstrahlung 
from produced pairs, or the production of electrons and muons not only as free 
particles but either produced into an atomic state bound to one of the ions 
(“bound-free pair production”) or as a bound state of the pair (“positronium”, 
“dimuonium”, or even “ditauonium” production). This will be discussed then 
in the rest of this chapter. For an introduction to aspects of atomic physics of 
relativistic heavy ion collisions, we refer the reader to [235,236]. 


7. 1 Strong Field Effects in Electron Pair Production: Multiple Pair Produc¬ 
tion 


The special situation of the electron pairs can already be seen from the formula 
for the impact parameter dependent probability for (single) pair production 
in lowest order. (For heavy nuclei the semiclassical approximation is valid 
and the impact parameter b can be considered to be an observable quantity, 
see, e.g.. Chapter 2 of [15] and further references given there.) Using the double 
equivalent photon approximation (valid for l/mg 5 'j/rne) with the cutoff 
radius set to l/rrie one obtains [237,238] 
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P^^\b) ^ [2 - 3 \n{meb)] ln{meb) for 1 < rrieb < jiab 

971^ mi¥ L J 

for 7 zab < me 6 < 7L (102) 

mjb^^ [ TTieO 

As this calculation underestimates the probability in the range of small impact 
parameters, one might expect that also probabilities larger than one are possi¬ 
ble. As mentioned already in the introduction, the use of the equivalent photon 
approximation is not justified here. Therefore (more) exact calculations need 
to be done. 

The impact parameter dependent probability in lowest order was calculated 
numerically in [239,240]. The quantity P^P{b) found there is larger than the 
one given by Eq. (102), and values larger than one are possible at RHIC and 
LHC. As an example P^'^^b = 0) = 3.9 (1.6) [79], and P‘'^\l/me) = 1.5 (0.6) 
[239] for LHC (RHIC) are found. It means that a new class of phenomena 
appears here. 

The fact that the cross section for the electron-pair production in lowest order 
rises too quickly with beam energy and will eventually violate unitarity at 
very high energies was already mentioned by Heitler in [241], and speculations 
about possible remedies were made there. Whereas it was assumed to be an 
“academic problem” at that time, the fact that unitarity is violated in the 
lowest order calculations at RHIC and LHC energies for the heavy systems 
has generated some renewed interest [10]. It led to a series of studies starting 
almost a decade ago. It was found that the production of multiple pairs in a 
single collision restores unitarity [242]. 

For Za ~ 1 the interaction of the leptons with several photons coming from 
the external held will be much more likely than the exchange of photons 
among the singly charged leptons, see Fig. 26. Therefore we will neglect the 
lepton-lepton interaction, treating it as a pure external held problem. 

This was hrst studied in [242,243]. Hereby, the e+e“ pair was treated as a 
“quasiboson”. Neglecting rescattering terms a closed expression was obtained 
to all orders. A similar approach was also followed in [244]. In [245] another 
approach was used, starting from a description within the Dirac sea picture. 
All studies essentially make use of the “quasiboson” approximation, that is, 
treating the pair as an “unbreakable unit”. All hnd the probability to produce 
N pairs P{N, b) to be given by a Poisson distribution: 

P(V()) = ^^hALexp[-p('>(6)], (103) 

where PP\b)^ the probability for pair production calculated in lowest order 
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Fig. 26. The dominant two pair production process at the energies of RHIC and 
LHC is shown in (a). Due to Za ~ 1 this process is favored, giving a cross sec¬ 
tion of ~ In^ On the other hand the process of (b) is proportional to 

~ ln^ 7 ;^^^. It will therefore take over at extremely high ')iah- 

(which, as shown before, can become larger than one), is actually the “average 
number of e’''e“ pairs” produced in a single ion collision: 

{*'((>)>=i: iVP(iV,&) = (104) 

N 

One should keep in mind that the two-pair production process as shown in 
Fig. 26(a) is rising approximately with ~ ln^( 7 ;^^^), whereas the one 

shown in Fig. 26(b) rises with ~ ln"^( 7 f^^) [246,247,59] and will hnally 

take over at extremely high values of 'jiab- We are still far away from this 
regime. On the other hand, as mentioned in [59] the multiple pair production 
at large invariant masses is again dominated by Fig. 26(b). 

The use of the Poisson distribution Eq. (103) allows us to interpret the lowest 
order calculation as a calculation of multiple pair production. Multiple pair 
production cross sections can be dehned by integrating Eq. (104) over b, for 
example: 


(Tj' 


d% {N(b)} = Na{N). 

N 


(105) 


It was shown in [248,234] that the matrix element for iV-pair production S'at 
factorizes quite generally into an antisymmetrized product of individual pair 
production amplitudes s’*"” (corresponding to individual fermion lines) and 
the vacuum-vacuum amplitude (corresponding to all closed fermion loops), 
see also Fig. 27: 


Sn = (0|^|0)^sgn(a)4^;^^^j (106) 

(7 


where /c*,/* are the quantum numbers (momenta and spin projection) of elec¬ 
tron and positron, respectively, and a denotes a permutation of {1, • • •, N}. 



The vacuum-vacuum amplitude is present in all QED calculations, but whereas 
it is of absolute value one without external helds or for static helds (and is 
therefore factored out and dropped in the calculation), it is < 1 here, as pair 
creation occurs out of the initial vacuum state and the probability for no-pair 
production |(0| 5* |0)|^ has to be smaller than one. This result was conhrmed 
recently in [249]. The amplitude corresponds to the one of single pair 
production neglecting the vacuum-vacuum amplitude; it is the one that is 
calculated, for example, in lowest order Born approximation. 

b 

ki 
ki 


h 
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Fig. 27. Graphical illustration of the general form of the iV-pair production process. 
The interaction with the external field is shown as crosses. The production of a pair 
is described by a fermion line coming from and leaving to the future, interacting an 
arbitrary number of times with the external field. Such a line corresponds to 
in Eq. (106). The vacuum-vacuum amplitude (0| S |0) corresponds to the sum of all 
closed fermion loops. For details we refer to [248]. 


Calculating the total N pair probability 

P(A'.(<) = 7^ E \S«\" (107) 


one gets two different types of contributions. In the first type — see Fig. 28(a) 
— the same permutation a is present in and S^, whereas different ones 
are used in the second type — see Fig. 28(b). Neglecting this second class — 
and therefore the antisymmetrization in the hnal state — one recovers again 
the Poisson distribution Fq. (103). 

The fact that the second class of diagrams gives in general a small contribution 
as compared to the hrst one can be understood due to the correlation of the 
electron and positron momenta. They are not produced completely indepen¬ 
dently of each other. Especially their transverse momenta are correlated; this 
is obvious in the equivalent photon approximation where the snm of the trans¬ 
verse momenta is close to zero, which is also true in the exact calculation in 
lowest order. Therefore the “quasiboson” approximation, that has been used 
in the early studies is justihed. In addition the deviation from the Poisson 
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Fig. 28. The two pair production probability consists of two types of diagrams. The 
class of diagrams (a), where electron and positron are always matched in Sat and 
leads to the Poisson distribution. Diagrams of type (b) are assumed to be small 
due to the correlation of the momenta of electrons and positrons. 

distribution in the case of two pair production was calculated explicitly for 
small impact parameters 6 0 [248,250]. Figure 29 shows the contribution 

from the two diagrams of Figure 28. At LHC energies the deviation is found 
to be only of the order of 1%. 



Y 


Fig. 29. The contribution of diagram (a) and (b) of Fig. 28 to the two-pair produc¬ 
tion probability at impact parameter 6 ~ 0 is shown as a function of the Lorentz 
factor jiab of the collision. The dotted line shows contributions from the diagram of 
Fig. 28(a), the numerical data points and solid line those of Fig. 28(b). Taken from 
[250]. 

Using Eq. (103) one can use the impact parameter dependent probability 
P{b) of [239,240] to obtain the probabilities for Al-pair production P{N,b), 
see Fig. 30. One can see that for impact parameters b ~ 2R up to about 1/me 
on the average 3-4 pairs will be produced in PbPb collisions at the LHC. This 
means that each photon-photon event — especially those with high invariant 
mass, which occur predominantly at impact parameters close to b ^ 2R — 
is accompanied by the production of several (low-energy) e’''e“ pairs (most of 
them however will remain unobserved experimentally). Integrating over the 
impact parameter the total multiple pair production cross section were given 
in [234]. The single pair production probability P{b) falls off essentially as 
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Fig. 30. The impact parameter dependent probability to produce N e+e“-pairs 
{N = 1,2, 3,4) in one collision is shown for the LHC = 2950,PbPb). Also shown 

is the total probability to produce at least one e+e“-pair P{N-, h) = 1 — P{0, b). 
One sees that at small impact parameters multiple pair production dominates over 
single pair production. 

1/6^. The total cross section is dominated by large impact parameters and 
not very sensitive to the multiple pair production effects at small b. In [238] 
the reduction of the exclusive one-pair production cross section (7{N = 1 ) 
was estimated to be -6.4% for RHIC and -4.7% for LHC, using the numerical 
results of [248,234]. Multiple pair production cross section on the other hand 
are dominated by impact parameters, around 6 ~ Ac = l/rrie- Measuring them 
would be of interest. At the SPS (CERN) the effect was looked for, see [251], 
but only an upper bound could be given, which is still above the theoretical 
prediction. 


1.2 Strong Field Effects in Electron Pair Production: Coulomb Corrections 


The calculation of the impact parameter dependent probabilities described 
above have been done in lowest order perturbation theory, even though Za ~ 
0.7 is not small. One could therefore ask how reliable these calculations are. 
Higher order effects (in addition to multiple pair production discussed in the 
previous section) could become important. This will be discussed now. 

For low beam energies (around 1-2 GeV per nucleon) coupled channel calcu¬ 
lations have been made by a number of people [232,233]. They found a rather 
large increase of the cross section. For the higher energies at RHIC and LHC 
such an approach is no longer viable due to the large number of channels to 
consider [252] (For a calculation at 7 = 200, see [253]). Since the review in [15] 
a number of papers appeared treating higher order effects in ultrarelativistic 
heavy ion collisions in the high energy limit. 

In the following we are concerned only with higher order corrections to the 
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single fermion line, that is to s’*"” of Eq. (106), assuming that the Poisson 
distribution can then be used to account for multiple pair production. These 
effects are generally called “Coulomb corrections”, as they mostly deal with 
the Coulomb distortion of both the electron and positron wave functions due 
to additionally exchanged photons. Formally two kinds of classes can be dis¬ 
tinguished. In the hrst one only Coulomb rescattering will be considered, that 
is, only one pair will be present at all intermediate time steps. The second one, 
which we call “multiple particle corrections”, are those processes, where more 
than one pair is present at an intermediate time step, with an electron and a 
positron from two different pairs annihilating at a later step. Such a distinction 
is less useful in the usual Feynman approach, as intermediate lines describe 
both electrons and positrons, but it is well dehned in a retarded boundary 
condition approach, corresponding to the Dirac sea picture. 

The second class of diagrams was studied in [248] for small impact parameters 
6^0 within the framework of the Magnus theory and neglecting all Coulomb 
rescattering terms. The effect was found to be rather small, at most 5% for 
PbPb collisions at LHC. As these effects are expected to be largest for the 
small impact parameters, one has to conclude that their effects will be rather 
small in the cross section. We therefore consider the hrst class in the following. 


A classical result for this hrst class of higher-order ehects can be found in the 
Bethe-Heitler formula for the process 'y + Z ^ + e~ + Z (corresponding to 

the highly asymmetric case Zia 0 , Z 2 a ~ 1 ): one obtains for an unscreened 
nucleus 


28 iZaf r 2uj 109 

-- In- 

9 rrig L rrie 42 


f{Za) 


( 108 ) 


with the higher-order term given by 
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(109) 


with the Euler constant 7 ~ 0.57721 and ip the Psi (or Digamma) function. As 
far as total cross sections are concerned the higher-order contributions tends 
to a constant for a; —> cx). 

This is used as the basis for a detailed calculation of higher order ehects on 
the (single pair) cross section in [254]. The higher order ehects were classi- 
hed according to the number of photons n and n' exchanged with each ion, 
see Fig. 31. The dominant term is given by the n = n' = 1 contribution 
(Fig. 31(a)), due to the 1/6^ behavior of the single photon exchange. It leads 
to the famous In^'y rise of the (total one pair production) cross section. The 
next important terms are those with n = 1 , n' > 1 and n > 1 , n' = 1 (see 
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Fig. 31(b) and (c)). They are enhanced proportional to ~ In^y. The last class 
of diagrams will only be ~ In 7 and is therefore neglected in [254]. A system¬ 
atic way to take the leading terms into acconnt in e+e” pair prodnction was 
used in [255,254], leading to a result which corresponds to the Bethe-Heitler 
result above and is therefore in accord with the Bethe-Maximon corrections 
[256-258]. The authors hnd corrections to the total Born cross sections, which 
are negative and equal to -25% and -14% for RHIC and LHC. Therefore the 
Coulomb corrections in the single pair production cross section (better ut 
as dehned in Eq. (105), as this is what is really calculated) seem to be well 
understood. A negative correction was also found in the work of [10]. These 



Fig. 31. Contributions to the single pair production cross section can be classified 
by the number of photons that is emitted by each ion. If each ion emits only one 
photon (a) the cross section has a In^ 7 dependence. If only one of the ions emits 
only one photon ((b) and (c)), only a In^y dependence is found. Finally with both 
ions emitting several photons (d), only a In 7 dependence is expected. 

theoretical developments are further supported in [237]. 

Please note that such a classihcation according to n and n', resulting in dif¬ 
ferent powers of In 7 is not helpful for multiple pair production, as at least 
two photons must be emitted from each nucleus in this case (corresponding 
always to the case n > l,n' > 1 above) and one will not have a cross section 
proportional to In^ 7 from the beginning. A different approach, which is good 
in the small impact parameter region, is therefore needed. 

In a series of papers the exact solution of the Dirac equation of an electron 
in the held of the two nuclei in the limit of the Lorentz factor 'yiab —> 00 
was studied. The interest started after an article showed how the summation 
to all orders in the high energy limit can be done in the related problem of 
bound-free pair production [259], see also below. In a hrst article [260] it was 
shown that also in the free pair production case the summation to all orders 
can be done analytically. Due to the form of the electromagnetic held of the 
two ions in the high energy limit — both of them being essentially localized 
on two diherent sheets corresponding to z = ±f, where is the direction of 
the beams — only a certain class of diagrams was found to be dominant in 
this case, see Fig. 32. The sum of the interactions with one of the ions in the 


80 



highly energetic case gives the typical eikonal phase. It was shown that this 
leads to a matrix element very similar to the one in lowest order, but where 
the photon propagator is replaced with 


^2 _i_ ^2 ( 2 I ^ 

<li. + Qi (gi + qt) 

The authors of [261] come to the same conclusion; they also show that by 
integrating over the impact parameter the cross section becomes identical 
to the lowest order Born result, as only the absolute values squared of the 
modihed photon propagators Eq. (110) are needed. The same result is also 
found in [262]. In [263] the scattering of electrons is studied, which is then 
related via crossing invariance to the pair production process. 

A numerical calculation of the impact parameter dependent probability was 
done in [264] and also in [265] (see also the two comments following [265]). 
Calculating the impact parameter dependent probability in this approach a 
reduction of up to 50% was found for small impact parameters, leading to 
the same reduction for the multiple pair cross section. Therefore Coulomb 
corrections are expected to be important also at small impact parameters. 

The fact that the cross section to all orders corresponds to the lowest order 
Born result is obviously in contradiction to the Bethe Maximon corrections, as 
described above. This can be seen from Eq. (7.4.3) of [10] and it was pointed 
out in [254] and later also in [237,266]. In [266] it was noted that in the limit 
of 7 —cx), different types of diagrams are of importance in the electron scat¬ 
tering case and in the pair production case, that is, taking only the dominant 
diagrams in both cases will violate crossing symmetry. In [237] the question of 
how to do this limit is studied in more detail. By making a careful transition to 
7 —cx) the Bethe Maximon corrections were found to be present also in this 
case. The failure of the eikonal technique to reproduce the Bethe Maximon 
results was already noted in [267]. 

The equality of the nonperturbative cross section to all orders with lowest 
order result can therefore be seen as an artifact of the regularization of the 
modihed photon propagator. Eq. (110) is found by taking the limit 7 —cx) 
hrst; in this limit the modihed photon propagator is Taking this 

limit — which essentially corresponds to the sudden approximation — is jus- 
tihed only for small impact parameters. In fact integrating this expression 
over b would lead to a diverging cross section, essentially due to the singular 
behavior of the propagator for g_L —0. Therefore in [261] the form of the 
propagator Eq. (110) was suggested. It is mainly this form of the propagator 
which in the end leads to the equality of the two cross sections. A more care¬ 
ful regularisation of the propagator would lead to an expression in agreement 
with the Bethe Maximon theory [268]. On the other hand this means that 


81 




z 


Fig. 32. The special structure of the electromagnetic interaction can be seen in this 
t-z plot. Due to the Lorentz contraction the electromagnetic fields are localized in 
two sheets corresponding lo z = ±t. In the “retarded” or “Dirac sea” approach (left) 
an electron with negative energies comes from t = —oo, crosses the field of each ions 
only once before leaving as an electron with positive energy. In the “Feynman” 
approach, the positron comes from t = +oo and can go forward and backward in 
time therefore interacting a number of times with the ions. 

the probabilities for small impact parameter, which get a large part of their 
contribution from large should be expected to be reliable. 

In all derivations of the pair production in the high energy limit [260,261,263] 
the pair production are not calculated using Feynman boundary conditions, 
that is, treating the positrons as particles going backward in time [269-272]. 
Instead a retarded approach, that is, essentially the “Dirac sea” picture is 
used: the problem is treated as the scattering of a negative energy electron 
into a positive energy state, see, e.g., [269,272]. In this approach the structure 
of the interaction with the external held can be exploited easily, see Fig. 32. 

This has two consequences: the total cross section as well as the differential 
cross section with respect to either the electron or the positron should be reli¬ 
able, but not the differential cross section with respect to both leptons, as the 
“rescattering of the hole” is not taken into account [237]. In addition in [249] 
it is shown that the cross section calculated in this way does not correspond to 
the exclusive one-pair production ai but to ar, as dehned in Eq. (105). This 
result was already found there to be true for the calculation of the single¬ 
pair production cross section, neglecting the vacuum-vacuum amplitude. But 
whereas the Poisson distribution was needed there, it is found in [249] to be 
exact in the general case as well. 

Single and multiple pair production in charged particle (heavy ion) collisions 
without nuclear contact is a QED process and should in principle be calcu¬ 
lable precisely with the methods of QED. Although the calculations in detail 


82 



are quite complicated, the theoretical concepts are well understood by now. 
Especially the importance of higher order corrections due to the multiple pair 
production cross section would be an interesting test. First experimental ob¬ 
servations of pairs has been reported already by STAR at RHIC [273]. 
A comparison with experiments will hopefully be possible soon. 


1.3 Electromagnetic Production of Electron Pairs and Other Particles in 
Central Collisions 


Differential production probabilities for dileptons from 77 processes in central 
relativistic heavy ion collisions are calculated using the equivalent photon 
approximation in an impact parameter formulation and are compared to Drell- 
Yan and thermal ones in [87,274,275]. The very low p± values and the angular 
distribution of the pairs give a handle for their discrimination to those coming 
from other sources, e.g., from meson decay. 


In [80] it is shown that particle production in central collisions via the 77- 
mechanism is so small that it can be neglected in practice. The probability to 
produce a hnal state with invariant mass M and rapidity Y is given by 


d^P 

dMdY 
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where ct^t is the corresponding 'yy ^ f fusion cross-section and a = {R 
is the larger one of the nuclear radii). This probability is very small, since we 
typically have utt ^ R^- Electromagnetically produced electrons were also 
measured in distant S — Pt collisions in the CERES/NA45 experiment [276]. 
The authors also conclude that “the rate of e+e“ pairs from QED produc¬ 
tion in collisions with nuclear overlap for the CERES acceptance is negligible 
compared to e+e“ pairs from meson decay f. For further details we refer to 
[80,276]. 


7.4 Electrons and Muons as Eguivalent Particles, Eguivalent Electron (Muon) 
Approximation 


The equivalent photon approximation can be generalized to other particles 
[277,278]. A light or zero mass particle moving with high energy splits into 
a pair of light or zero mass particles. The constituents of this daughter pair 
subsequently react as ’’equivalent beams” with the target system. This sub¬ 
process is easier to describe than the original reaction. Here we are interested 
in those processes, where the (equivalent) photon splits into a lepton pair and 
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the (equivalent) leptons induce reactions with the target, e.g., deep inelas¬ 
tic lepton scattering. (The discussion here has some similarity to the more 
general one in Chap. 3 above, where also vector mesons and quark pairs are 
considered.) 

The production of lepton pairs from (real or virtual) photons producing a 
pair by inelastic scattering from a nuclear target was first studied in [279], see 
also [69], where it was mostly regarded as an additional contribution to the 
(elastic) pair production. Now we suggest that the equivalent photons may 
also be used as a lepton beam to study deep inelastic lepton scattering from 
a nucleus. 


Up to now only the production of dileptons in heavy ion collisions was con¬ 
sidered, for which the four-momentum of the photons was less than about 
1/7?^ (coherent interactions). There is another class of processes, where one 
of the interactions is coherent {Q^ < 1/7?^) and the other one involves a 
deep inelastic interaction (Q^ 3> 1/77^), see Fig. 33. These processes are much 
rarer than the events with two coherent interactions. But they may also be of 
interest in the future. They are readily described theoretically using the equiv¬ 
alent electron- (or muon-, or tau-) approximation, as given, e.g., in [277,278]. 
They are characterized by the fact that the lepton is almost on-shell (i.e. the 
propagator of the virtual particle will become very large, since its momentum 
squared |p^| ~ w?). In the following we will speak of muons, but the same 
considerations apply also to electrons. 

The equivalent photon can be considered as containing muons as partons. The 
equivalent muon number is given by [277]: 
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where denotes the muon mass. The energy of the equivalent muon 
is given by 77^ = xu, where u is the energy of the equivalent photon. This 
spectrum has to be folded with the equivalent photon spectrum given approx¬ 
imately by, see Eq. (11): 
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we then get: 


Rmax 

//i/z(a^i) = J duf^/z{u)f^/^{uEA,Xi/u). ( 114 ) 
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The energy of the heavy ion is denoted by Ea-, the photon energy is given by 
ui = Mi?A,and the muon energy by = XiEa = uxEa- This expression can 
be calculated analytically. 


The deep inelastic lepton-nucleon scattering can now be calculated in terms 
of the structure functions Ei and E 2 of the nucleon. It would be interesting 
to see how the nucleon structure functions Ei and E 2 are modihed in the 
nuclear medium. (It remains however unclear at present how accurate such 
measurements could be.) The inclusive cross section for the deep-inelastic 
scattering of the equivalent muons is given by 
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where can be calculated from the usual invariant variables in deep 

inelastic lepton scattering (see, e.g., Eq. (35.2) of [280]). Here the scattering 
angle of the lepton is 6 and its energy E'. The accompanying muon of opposite 
charge, as well as the remnants of the struck nucleus, will scatter to small 
angles and generally remain unobserved, see Fig. 33. The hadrons scattered 




Fig. 33. Events where the “equivalent muon” is scattered deep inelastically on the 
other ion are characterized by one muon being emitted with a large transverse mo¬ 
mentum, whereas the other muon escapes unobserved (together with the remnants 
of the struck ion) in the forward direction. Energy and momentum of the struck 
parton can be used to infer the energy of the (initial) equivalent muon. 

to large angles can be observed, with total energy Eh and momentum in the 
beam direction of Pzh- One can use the Jacquet-Blondel variable i/jb dehned 
as 


yjB 


Eh - Pzh 

2E^ 


(116) 


The missing hadrons which do not enter the detector (which we assume to 
be almost 47r) have a small pt, and thus their contribution to the numerator 
Eh — Pzh is small (ultrarelativistic particles are assumed). One gets a good 
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estimate of y using this method. Combining this with Eq. (3.27) of [207] the 
initial energy of the equivalent muon can in principle be reconstructed as 

^ Eh - Pzh + cos 9) . (117) 

This is quite similar to the situation at HERA, with the difference that the 
energy of the lepton beam is continuous, and its energy has to be reconstructed 
from the kinematics (How well this can be done in practice remains to be seen). 
Of course the struck parton also has to go through the nucleus, and there is 
the possibility of ESI with it. This could introduce some uncertainty in the 
reconstruction. 

Quite similar also the excitation of individual nuclear states is possible in this 
way. As these states can decay electromagnetically, they are a possible source 
of high energetic photons. Recently this excitation was discussed [281] for a 
related process, where pairs are produced first and then excite the nnclear 
levels in a second step. It is astonishing that the anthors fonnd a qnite large 
cross section of 5.1 barn for this process in the case of CaCa collisions at the 
LHC. 

7.5 Radiation from e^e~-Pairs 

The bremsstrahlung in very peripheral relativistic heavy ion collisions was 
found to be small, both for real [10] and virtnal [282] bremsstrahlnng photons. 
This is due to the large mass of the heavy ions. Since the cross section for 
e'''e“ pair prodnction is so large, one may expect to see sizeable effects from 
the radiation of these light mass particles. For e’''e“ colliders it was fonnd that 
“at large values of s the photon emission proeess in electroproduction of the 
e~^e~ pair ... becomes very important (and in a number of cases, decisive)’’ 
[283]. 

In the soft photon limit (see, e.g., [284]) one can calculate the cross section 
for pair prodnction accompanied by soft photon emission 

Z + Z—i>Z + Z + e’'' + e“+7 (118) 

as 

da{k,p_,p+) = -e'^ -^dao(p+,p_), (119) 

P-K P+k\ 9:71^00 

where dao denotes the cross section for the pair production in heavy ion 
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collisions (without soft photon emission). The corresponding Feynman graphs 
are given in Fig 34(a) + (b). An alternative approach is to use the double equiv¬ 
alent photon approximation and calculate the (lowest order) matrix element 
for the QED subprocess 


7 + 7 




+ e +'y. 



Fig. 34. Emission of a bremsstrahlung photon from pairs produced electromagneti- 
cally. Process (a) and (b) are dominant in the infrared regime. For the calculation 
either the IR approximation together with a full calculation of the pair production 
was used or a full calculation of all three diagrams together with photon spectra 
from double equivalent photon approximation (DEPA). Similar diagrams, where the 
photon line is attached to the ions are generally small due to the heavy mass of the 
ions. They are neglected here. 


In Fig. 35 we show results of calculations for soft photon emission. The cross 
section for the QED process (Eq. (7.5) above) in lowest order was used and 
folded with the corresponding (double) equivalent photon spectra [285]. 



Fig. 35. The energy-dependence of bremsstrahlung-photons from e'’'e pair pro¬ 
duction is shown for different angles. We show results for PbPb collisions at LHC. 


These low energy photons might constitute a background for the detectors. 
Unlike the low energy electrons and positrons, they are not bent away by the 
magnets. The angular distribution of the photons also peak at small angles, 
but again a substantial amount is still left at larger angles, even at 90°. The 


87 
























typical energy of these low energy photons is of the order of several MeV, i. 
e., much smaller than the expected level of the energy equivalent noise in the 
CMS ECALs [286], 

7.6 Bound-free Pair Produetion 


The bound-free pair production, also known as “electron-pair production with 
capture”, is a process which is also of practical importance for relativistic 
heavy ion colliders. An electron-positron pair is produced, where the electron 
is not free but bound in an atomic state of one of the ions: 

Zi -\- Z 2 ^ {Zi e + e"'' -|- Z 2 - (120) 


As this changes the charge state of the nucleus, it is lost from the beam in 
the collider. Together with the electromagnetic dissociation of the nuclei (see 
Sec. 5.1 above) these two processes are the dominant loss processes for heavy 
ion colliders. It has only been realized recently [169], see also [92,287], that this 
process can also result in a localized beam-pipe heating: the atomic states are 
produced with a perpendicular momentum of the order of rrie, which is rather 
small and therefore leads to a narrow singly-charged ion beam. These beams 
with altered magnetic rigidity will deposit their energy in a localized region of 
the beam pipe and cause a localized heating. This can lead to the quenching 
of the superconducting magnets. The energy deposited per unit time is of the 
order of Watts. This limits the luminosity of the PbPb collider at LHC [92]. 
Due to the lower beam energy at RHIC, the energy deposit there is smaller 
and therefore less important. 


From a simple kind of coalescence model, a scaling law for the cross-section 
an for bound-free pair production of the form 


O'n ~ 


z\zl 


n'" 


( 121 ) 


can be derived, where n is the principal quantum number of the bound state 
and only s-states are populated [288,289]. It is easiest to see this in a system 
where the nucleus Zi is at rest. 

A typical graph is shown in Fig. 36. An atom {Zi + e~) with momentum k 
is produced. The bound state wave function is a strongly peaked function of 
the relative momentum Krei = I — {.Mzi/Mf^Zi+e-))k ^ I — k where I is the 
three-momentum of the photon exchanged with the target. This momentum 
is of the order of rUe, while the momentum Krei is of the order of ZiarUe 
(i.e. Zi times the Bohr momentum psohr = crnie-)- Since this momentum I is 




Fig. 36. A typical graph describing the bound-free pair production process, see 
[288]. 

usually much larger than Krei the loop integration over the relative momentum 
factors out and the factor 'h(O) ^3/2 

is found, where d' is the bound state 
s-wave function in coordinate space. The cross section to produce free pairs is 
proportional to ZfZ^; so the scaling law Eq. (121) above is obtained directly 
from the scaling properties of the hydrogen-like wave-functions. 

The calculations to be discussed in the following are not done along these 
lines. Instead the matrix elements involving the full wave functions of the 
bound and free lepton in the held of nucleus Zi are evaluated directly. As will 
be seen below, these scaling rules still hold to a good degree of accuracy. Not 
unexpectedly deviations are found for the large values of Zi , where Zia . 
In these cases the argument given above is not so well fulhlled {urei ~ but 
not TUe). But still, the 1/n^ scaling is also very well fulhlled in this case, 
see below (and Fig. 37). 



0 10 20 30 40 50 60 70 80 90 
Z 


Fig. 37. The cross section for the bound-free pair production process are shown for 
LHC energies (using 7 = 3400) as a function of the charge of the ion Z. The cross 
sections for capture to a given state with quantum numbers n/, Kf are scaled by 
the factor , see Eq. (121). 


Recently, the cross-section for the process Eq. (120) was calculated in PWBA 
using exact Dirac wave functions for a point nucleus [290]. The formalism 
of [291] was extended in a straightforward way to the heavy ion case. The 
cross-section for capture to the n, k state can be written as 
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Here /Sjon and 7 jon are the velocity and the Lorentz factor of the projectile ion 
in the rest frame of the target ion, k± denotes the perpendicular momentum 
of the exchanged photon, kz = uj/P ion is its component in the beam direction 
and the energy is given by u. The quantities Tj_ and Ti are functions of u and 
Q^, independent of 7 . For the integrations over k± and u a logarithmic grid 
in u and k±, with energies and momenta up to 500 rrie and angular momenta 
up to / = 200 was used. 

In addition to lsi /2 capture the capture to the 2 si/ 2 , 2 pi/ 2 , 2 ^ 3/2 and 3 si /2 
states was also calculated explicitly. The only approximation in these calcula¬ 
tions is the neglection of higher order effects (in Z 2 ). In Fig. 37 we show the 
dependence of the cross-section (divided by the approximate scaling factor Z'^) 
as a function of Z for various bound states (We consider the symmetric case 
Zi = Z 2 = Z; a.s the dependence on the projectile charge (in lowest order) is 
given by — see the scaling law Eq. (121) — the asymmetric cases can be 
found as well.) We especially note the rise of the cross-section for pi /2 states. 
This is due to the s-wave character of the small component of the bound-state 
wave-function. We also see a deviation from the Z"^ scaling, which is due to 
the behavior of the Dirac wave functions for small values of r [292]. While the 
bound state wave function is further increased as compared to the nonrela- 
tivistic Z^ rise, the Coulomb repulsion of the positron causes the Si /2 curves 
in Fig. 37 to decrease with increasing Z. The 1/n^ scaling with the principal 
quantum number n for the s-states is very well fulhlled in these calculations 
for all values of Z. 

The cross section per target in Eq. (122) can be shown to be of the form 

a = Aln'yiab +B = Cln'yion +D. (123) 


This form has been found to be a universal one at sufficient high values of 7 . 
The constant A and B (and alternatively C and D) then only depend on the 
Z value of the target. Parameterizations for A and B [290] for typical cases 
are given in Table 5. 

The calculation of bound-free pair production has a long history. Let us men¬ 
tion here the papers [293-298]. Most of these calculations have either not been 
done at the high beam energies relevant for RHIC and LHC or are based on 
the EPA. In [10] an analytical formula for bound-free pair production was 
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Ion 

A 

B 

(^,llab ~ 100 

0 -, Jlab ~ 3400 

Pb 

35.5barn 

— 22 . 1 barn 

146 barn 

272 barn 

An 

28.7barn 

— 17.8barn 

114 barn 

212 barn 

Ca 

4.46mbarn 

— 2 . 88 mbarn 

18.7 mbarn 

33.9 mbarn 


Table 5 

Parameters A and B (see Eq. (123)) as well as total cross sections for the bound- 
free pair production per target for RHIC and LHC. The parameters are taken from 
[290]. Capture to the Is-, 2s-, 2 pl/2-, 2p3/2-, and 3s-states are taken into account. 

given, starting from approximate wave functions for electron and positron. 
For large values of Z this formula underpredicts the cross-section. The reason 
for this was recently discussed in [299]. All other results were found to be in 
agreement with each other [290]. 

Most calculations were done in hrst order in the interaction with the projectile 
only. For a long time the effect of higher order processes have been under inves¬ 
tigation. Such higher order processes would violate the scaling in Eq. (121). 
At lower beam energies, in the region of a few GeV per nucleon, coupled chan¬ 
nel calculations have indicated, that these give large contributions, especially 
at small impact parameters, increasing the cross section by orders of magni¬ 
tude. Newer calculations tend to predict considerably smaller values at higher 
energies, indicating that the number of channels used in the original calcula¬ 
tions was not sufficiently large enough [252]. In a recent calculation such effects 
were studied at high beam energies [259] in the high energy limit 7 —cx). In 
these calculations the higher orders were found to be quite small (less than 
1%) and even tend to reduce the cross section. Higher order effects are more 
important for the smaller impact parameters b ~ l/rrie, where the helds are 
largest. It seems interesting to note that very large (almost macroscopic) im¬ 
pact parameters contribute: for a; = 1 MeV ~ 2me and yjon = 2 x 10^ (relevant 
for PbPb collisions at LHC) one has bmax = ~ 4/im. 

There are experimental results of bound-free pair production for hxed targets 
at the Bevalac, AGS, and SPS/CERN. At Bevalac energies, the experimental 
results [300-302] are not quantitatively reproduced by lowest order theories. 
Higher order effects are present for the systems with Za ^ 1, see the conclusion 
of [302]. At the higher AGS energies [303,304], higher order effects become 
smaller and there is good overall agreement with theory. A cross-section of 
8.8 b was found for AuAu collisions at = 12 . 6 , in agreement with the 
theoretical results of [290,295,298]. Screening effects due to the target electrons 
will be small, since only small impact parameters contribute appreciably. At 
the even higher SPS energies, the effects of the target electrons (screening, as 
well as, “antiscreening”) should be taken into account [305-307]; for further 
details see the references given above. The bound-free pair production was 
measured at the SPS for 7 = 168 with a hxed target [308,309]. There is an 
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overall agreement with theoretical expectations, for details we refer to these 
papers. 

In principle the same process is also possible for muon-production. Due to 
their larger mass and therefore smaller Bohr radii, muon-capture is much 
more sensitive to the hnite size of the nucleus. But the cross section for this 
process are rather small. 

The process corresponding to Eq. (120) with one of the ions replaced by an 
antiproton p was used to produce fast antihydrogen [288,310,311]. In [311] 
a value of 1.12 ± 0.14 ± 0.09 pb was given for the cross section for Z 2 = 
1. The momentum of the antiproton beam was in the range between 5203 
and 6323 MeV. Due to experimental limitations all capture into the Is-state 
and 98% of those into the 2s-state but no higher states were measured. This 
experimental value is in agreement with the theoretical ones of 0.89 pb [291] 
and 1.02 pb [312] for 7 ^ 0 ^ = 6 . (The contribution of the capture to the 2 si/ 2 - 
state is given by the 1 /v? scaling law, as discussed above.) 

Finally we discuss a practical aspect. When passing through the magnetic 
helds in the collider, it is possible in principle that the atom formed in the 
interaction region will be subsequently ionized in the magnetic helds of the 
system. Since this would alleviate the problems related to beam pipe heating 
(see the discussion above) we give now a simple estimate, which shows that 
this effect will be negligible in practice. Let us assume that the {Z + e~) atom 
moves through a constant magnetic held of strength B. In its rest-frame this 
corresponds to an electric held of strength E = 'yiabl3iabB, i.e., a linear term 
eEx is added to the Coulomb potential of the electron in the held of nucleus Z. 
This could cause ionization by tunneling through this barrier, see, e.g., [292]. 
For a given binding energy E = IRy with IRy = 13.6eR, this escape of 
the electron from the Coulomb barrier is even allowed classically for a certain 
critical held. For a given Z and electric held E a critical value Ucrit is found, 
above which the atom ionizes in this classical approach. Such held ionization is 
well known for the case of the loosely bound Rydberg atoms. This ionization 
will happen very fast, it is of the order of the classical orbiting time of an 
electron in the corresponding state with principal quantum number n. We 
hnd 


nV, = 3.2 X lO'zyB, 


(124) 


where E is the held strength given in units of R/cm. See also Eq. (54.2) of 
[292]. As an example we take Z = 82, E = 10^'^R/cm which corresponds to 
a magnetic held of R = 1 Tesla and 7 = 3000. We have Ucrit ~ 12, i.e. only 
states with very high principal quantum number will be ionized. According to 
the 1 /n^ scaling law, their contribution can safely be neglected. 
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7.7 Formation of Dilepton Bound States: Positronium and Muonium 


Finally we note that the electron and the positron can also form a bound 
state: positronium. This is in analogy to the 77 -production of mesons {qq 
states) discussed in Sec. 3.6. With the known width of the parapositronium, 

= (125) 

the photon-photon production of this bound state was calculated in [313]. The 
production of orthopositronium, n = 1 ^Si was calculated recently [314] and 
[315], see also [316] where the three-photon production of vector mesons is 
treated. In [317] an effective parameter p is introduced, which controls the 
relative importance of two- and three-photon- processes (see also p. 1677 of 
[15]): 


P = 



(126) 


with 1 /hf = 1/6 (r^) and (r)^ is the mean square radius of the charge distri¬ 
bution and is the mass of dimuonium (or any other produced particle). 
This factor p is always much smaller than 1 except for positronium, where p 
is given by Za (in this case K k, 1 /R should be replaced by A ^/rne, see 
the discussion at the beginning of this chapter). Therefore the production of 
orthopositronium is only suppressed by the factor (Za)^ as compared to the 
two-photon production of parapositronium. Since (Za)^ ~ 0.5 for the heavy 
systems like An or Pb, one expects that both kinds of positronium are pro¬ 
duced in similar numbers. Detailed calculations show that the three-photon 
process is indeed not much smaller than the two-photon process [314,318]. 

The bound pFp,~ system (“dimuonium”) is an interesting example of a pure 
electromagnetic bound state, which can be effectively studied in heavy-ion col¬ 
lisions. Dimuonium has a mass of about two muon masses, ^211 MeV 
and the binding energy in the ground state is of the order of 1.4 keV. Although 
this system is well known theoretically, see [319-324], dimuonium has not been 
observed experimentally up to now. Like in the case of positronium, one should 
distinguish para- and ortho- dimuonium with = O’*" and 1 “ respectively, 
which can be produced in electromagnetic decays of light mesons, i.e. ortho- 
dimuonium in radiative decays of 7 - [322,325] and K^- mesons [323], while 
paradimuonium in decays of uj- and 0- mesons. The hyperhne structure and 
the decay rate of dimuonium are influenced by the electronic vacuum polariza¬ 
tion in the far time-like asymptotic region, which does not yield any contribu¬ 
tion in any other bound states [326,327]. In particular, a relative contribution 
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of 1.6 % to the decay rate of orthodimuonium is predicted in [320,321] and 
this would allow a test of QED in a previously unexplored kinematic region. 

Dimuonium production in heavy ion collisions is studied in [328] in detail. Be¬ 
low we will mainly follow this paper. Paradimuonium is produced dominantly 
in two-photon processes, orthodimuonium in three-photon processes while its 
production by bremsstrahlung is suppressed, see Fig.38. The paradimuonium 



(a) (b) (c) 


Fig. 38. Diagrams for two-, one- (bremsstrahlungs) and three- photon production 
mechanism of a dimuon bound state in heavy ion collisions [328]. 

production cross-section is determined according to Eqs. (67) and ( 68 ) by the 
77 -width, which in the leading order approximation is given by Eq. (125) with 
rUu instead of rUe- Orthodimuonium cross section is suppressed compared with 
that of paradimuonium by more than an order of magnitude see Eq. Al of 
[328]. The production cross sections are given in Table 6 . The main decay 



a{AA^ AA + PDM), pb 

a{AA AA + ODM), pb 

AuAu, RHIC 

0.15 

0.021 

PbPb, LHC 

1.35 

0.089 

CaCa, LHC 

6.6x10-3 

6.9x10-3 


Table 6 

Production cross sections for paradimuonium (PDM) and orthodimuonium (ODM) 
in heavy ion collisions, as predicted in [328]. 

channels are the annihilation processes into two photons for paradimuonium 
and e’''e“ for orthodimuonium. Dimuonium can therefore be observed via these 
decay channels. Further details, especially concerning background processes, 
can also be found in [328]. 

One might even think of the production of ditauonium (DT). With the di- 
tauonium mass uidt = 3554 MeV, its width (see Eq. (125)) is given by 
0.018 eV. From these values we can immediately obtain the cross section of 
ditauonium production in 77 -interactions in heavy ion collisions. For exam¬ 
ple, in PbPb collisions at the LHC we obtain from Eq. 68 and Fig. 14 the 
production cross section = 0.7 nb. We note that the lifetime of the r of 
291 X 10“^^s corresponds to the width of r-lepton P,- = 0.00226eV, i.e., the 
ditauonium 77 -width is about 5 times larger than the width of its weak decay 
given by 2 Vt- and thus the ditauonium can be really produced as a bound state 
of two r-leptons. 
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8 Methods of Detecting Very Peripheral Collisions in Experiments 
at Heavy Ion Colliders 


In this chapter we discuss the study of 77 -events in heavy ion collisions from 
an experimental point of view. We discuss the general characteristics of these 
kind of events and how they can be used as basic signatures for the triggering 
and off-line selection of these events. We discuss further various background 
sources for 77 -processes, including single-photon, photon-Pomeron and double 
Pomeron processes, their corresponding cross sections and expected trigger 
rates. We also give an overview of the experimental methods used for the 
selection of different hnal states produced in 77 interactions, based on the 
experience of current e’''e“-colhder experiments, as well as, on the hrst studies 
of very peripheral events in heavy ion collisions at RHIC and Monte-Carlo 
simulations done for the future heavy-ion experiments at LHC. 


8.1 Signatures of -Processes in Heavy Ion Collisions 


The experimental technique for the detection of 77 -processes is based on spe¬ 
cial observables (signatures) which make the 77 -processes different from other 
very peripheral processes, as well as, from strong peripheral interactions of the 
ions. In particular one of the main features of coherent electromagnetic inter¬ 
actions of relativistic ions is the small transverse momentum transfer from the 
ions. As a result the ions remain intact after such an interaction and in general 
they will be kept for some time in the beam pipes and thus cannot be detected 
in an experimental setup. Therefore very clean events will be detected with 
particles produced mostly in the central rapidity region. (In Sec. 2.5 we al¬ 
ready mention the interesting possibility to tag the scattered protons in pp 
collisions studied in [76]. In the case of ions this is probably impossible to 
achieve, due to the low limits on x, see Sec. 1 above.) 

It is useful to mention in this context also the excitation of the Giant Dipole 
(GDR) and other higher Multipole and/or Multiphonon Resonances (MR), 
which lead in their decay to the emission of one or several low-energy neutrons 
in the nucleus rest frame, see Sec. 5.1. In the laboratory frame — due to the 
very high Lorentz factor of ions and the low momenta of the neutrons in the 
ion frame — one has almost monoenergetic neutrons, which can be detected 
in the very forward hadron calorimeters commonly called the Zero Degree 
Galorimeters (ZDG) of the experimental setup by the deposit of multiples 
of the energy 'jiabmn for each detected neutron. These almost monoenergetic 
signals in the ZDGs, see Figs. 19 and 20, and their small number are therefore 
a good signature of the electromagnetic nature of the interaction between the 
two ions. 
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Another feature of 77 -processes in heavy ion collisions is a relatively low mul¬ 
tiplicity of secondary particles compared to those in strong ion interactions. 
Fig. 39 shows the average charged particle multiplicity for the process 77 —X 
as a function of the effective 77 -mass as predicted by PYTHIA [329]. This 
should be compared with the multiplicity expected in strong peripheral ion 
collisions. As an example the expected average charged particle multiplicity 
distribution vs. pseudorapidity dNch/dr] in PbPb collisions with a large im¬ 
pact parameters b between 1.8 x Rpi, < b < 2 x Rpi, is shown in Fig. 40 for 
LHC energy. This multiplicity distribution has been calculated using HIJING 
[330]. The total charged particle multiplicity in hadronic peripheral collisions 
is much higher than that in 77 -processes; even if one is using a restricted 
aperture, let’s say, |? 7 | < 1 , the multiplicity is still expected to be of the order 
of hundreds and therefore remains very high. The most important signature 



Fig. 39. The average charged parti¬ 
cle multiplicity in 77 -interactions vs. 
the effective 77 -mass predicted by the 
PYTHIA model [329]. 



Fig. 40. Charged particle multiplic¬ 
ity dNch/dr] in hadronic peripheral 
PbPb collisions with impact parame¬ 
ter 1.8 X Rpf) < b < 2 X Rpi, calcnlated 
with the HIJING generator [330]. 


of 77 processes in coherent ion-ion collisions is the quite low total transverse 
momentum Pj_ of the produced system. As discussed in Chap. 1, see Eq. (4), 
the 77 system has a transverse momentum less than about \/2/R, where R is 
the nuclear radius. But the average transverse momenta is even smaller than 
this, namely 


Pj_ ~ M^^/'yidh. 


(127) 


For the effective 77 mass range at RHIC < 3 GeV), as well as, at 

LHC (My.y < 100 GeV) we have P± ~ 30 MeV. It should be remembered 
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however that transverse momentum distributions have not yet been calculated 
explicitly within the semiclassical approximation. The problem is formulated 
in [ 88 ]. In order to use this feature for selecting 77 -processes the experimental 
setup has to provide a large enough aperture for the detection of both charged 
and neutral particles, including also secondary photons. It should provide also 
a sufficiently good transverse momentum resolution. 

Many of the 77 reactions, which are planned to be studied, are exclusive 
processes. In this case the net charge of the detected particles or the number 
of charged tracks can be taken into account. This leads to one more selection 
criteria for 77 -processes, requiring an even number of charged track of opposite 
signs. 

We therefore may formulate the criteria for the selection of the two possible 
classes of two-photon processes in collisions of relativistic ions: 

• The pure jj-processes are characterized by a small multiplicity of the charged 
particles in the central detectors, a small total transverse momentum of the 
produced system and the absence of any signal in both ZDCs; 

• The '-y'^-processes with EM excitation of the secondary ions are character¬ 
ized as well by a small multiplicity of the charged particles in the central 
detectors, a small total transverse momentum of the produced system and 
by neutrons coming from the GDR or MR decay in one or both ZDCs. 

Actually these are ideal criteria for the event identihcation, which in real 
experiments are often less clear due to different limitations of the experimental 
setup that lead to distortions and mixtures between the different event classes, 
as well as, to the contamination of them by background processes. 


8.2 Background Processes 


In this section we discuss different sources of background for very periph¬ 
eral processes in relativistic heavy ion collisions. As mentioned already in 
Chap. 6 , one should distinguish between photon-photon, photon-Pomeron and 
Pomeron-Pomeron processes. All of them are quite interesting in themselves 
and they also have similar signatures. The physical nature of them is quite 
different, and each of these processes should therefore be considered as a back¬ 
ground for the others. All those processes belong to the class of very peripheral 
collisions, taking place at impact parameter b > 2R. The strong interaction 
of the ions with impact parameter b ~ 2R is a common source of background 
for all very peripheral processes and should be taken into account. This is a 
physical background, as it results from the interaction of the ions with each 
other. 
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On the other hand there is background coming from the apparatus, which is 
then specihc for each experimental setup. As potentially important sources 
of such a background we discuss in the following one-photon processes, the 
photonuclear excitation of the ions and multiple e’''e“ pair production. Cosmic 
ray events and beam-gas interaction should be mentioned also in this respect 
as possible background sources of the apparatus. 


8.2.1 Physical Background 

Two-photon processes in relativistic ion collisions are discussed in detail in 
Chap. 3 including the corresponding cross sections and rates. Pomeron-Pomeron 
processes as a possible background have been discussed in Chap. 6. 

As was mentioned there the Pomeron-Pomeron cross section increases with 
^1/3 Qj, equivalently with compared to the rise of the photon-photon 
processes. Therefore in heavy ion collisions we can expect the predominance of 
two-photon processes, whereas in light ion and especially in pp-collisions the 
contribution of Pomeron-Pomeron processes will be much more essential. As 
shown in Chap. 6 this is conhrmed in detailed calculations, see, e.g., [21,227]. 

The cross section for the photon-Pomeron processes in ion-ion collisions can 
then be estimated with the help of the factorization relations [115,116], as 
explained in Sec. 3.5, using the cross section for photon-photon and Pomeron- 
Pomeron interactions: 

aU-fP ^X) = atoth-f ^ X) aUPP ^ (128) 


which is valid also for the corresponding differential cross sections. This can 
be used to estimate the total cross section of the photon-Pomeron processes 
in relativistic ion-ion collisions using the known cross sections of the photon- 
photon and Pomeron-Pomeron processes as = ^a^^a-p-p. 

The cross sections of the three different processes for PbPb, ArAr and pp 
collisions at the LHC, as well as, for AuAu, CuCu and pp collisions at RHIC 
are shown in Table 7. 

From this it is seen that at the LHC in the very peripheral PbPb collisions 
we have dominantly two-photon processes, in pp collisions on the other hand 
double Pomeron processes dominate, while in the ArAr collisions cross sections 
of all three very peripheral processes are comparable. Similar cross sections 
ratios for different ion species are valid also at RHIC. 

Apart from the different A-dependence of the very peripheral processes the se¬ 
lection of one of the processes can be enforced by using the following additional 



signatures: 

(1) A different total Pt distributions for very peripheral processes is expected 
because in the case of Pomeron-Pomeron processes the Pt are in general 
higher than in the case of photon-Pomeron processes, and these are again 
higher than those in photon-photon processes, see [24], 

(2) The different processes will lead to a different C-parity of the (exclu¬ 
sive) system produced in the different very peripheral processes. C-parity 
should be even in the two-photon and double-Pomeron processes, while 
it is odd in the case of photon-Pomeron processes, 

(3) The probabilities for the excitation of the GDR and MR of the ions, that 
is, of the additional electromagnetic excitation of the nuclei should be 
a good sign for a very peripheral collision, discriminating them against, 
e.g., grazing collisions or diffractive processes. The probability to excite a 
GDR or MR in the ions is much higher for an additional photon-induced 
reaction than for a hadronic collision. The GDR favors neutron over pro¬ 
ton emission, whereas a direct hadronic nucleon emission is insensitive 
to proton/neutron, and thus the emission of a soft proton in the ion rest 
frame and detection of the almost monoenergetic proton in the charged 
ZDG could be regarded as a sign of a Pomeron-ion interaction. 



Reaction 

AA^X 

AA AAX 

AA AAX 

AA AAX 



(nuclear) 

(FF^ hadrons) 

( 77 ^ hadrons) 

( 7 F^ hadrons) 

L 

cr(PbPb), barn 

7.8 

00 

h-^ 

0 

1 

1.5- 10-1 

1.1 • 10-2 

H 

O'(ArAr), barn 

2.1 
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6.1 • 10 -^ 

C 

cr{pp), barn 

0.07 

5.2 • 10“^ 
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2 . 8 - 10-3 

R 

cj(AuAu), barn 

7.5 

2.2 • 10 "^ 

1 . 8 - 10-3 

6.3 • 10-^ 

H 

cj(CuCu), barn 

3.1 

2.1 • 10 "^ 

9.3- 10-3 

1.4- 10-^ 

I 

cr{pp), barn 

0.07 

to 

h-^ 

0 

1 

1.7- 10-1* 

6.1 • 10-'^ 

C 






Tab 

le 7 


Cross sections of the nuclear, Pomeron-Pomeron, photon-photon and photon- 
Pomeron processes (M > 1 GeV) in PbPb, ArAr, pp collisions at the LHC and 
in AuAu, CuCu, pp collisions at the RHIC. 


Please note that in the case of heavy ion collisions, due to the large probabil¬ 
ities for electromagnetic excitations, see Sec. 5.1, additional electromagnetic 
processes, which we illustrate for the case of p'^-production in AuAu collisions 
at RHIG in Fig.41, are very likely. In heavy ions, the probability for the excita¬ 
tion of the GDR is large for small impact parameter, but this is not essential 
for the selection of 77 -processes because for heavy ions 77 -processes domi¬ 
nate over the processes involving the Pomeron. The signature 3 above could 
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Fig. 41. Diagrams of production in heavy ion collisions when neither, one or 
both nuclei are excited electromagnetically. Diagrams are show up to vertex permu¬ 
tations. 

become important in the case of intermediate ion collisions, where on the 
one side the cross sections of photon-photon, photon-Pomeron and Pomeron- 
Pomeron processes are comparable, bnt on the other side the cross section of 
GDR excitation is low as well, and thns the analyzing power of signatnre 3 
becomes very high. 

Concerning the strong interaction of the two ions, these events in general 
are qnite different in their properties from the very peripheral ones discnsses 
above, compare, for example, Figs. 39 and 40. Nevertheless in the case of 
detectors covering only a small apertnre range, like the central detector and 
especially the ZDC, the backgronnd from these processes shonld be taken into 
acconnt. 


8.2.2 Apparatus Backgrounds 

Since the cross sections of the interesting coherent processes in relativistic very 
peripheral ion-ion collisions are mnch smaller than the total cross sections, the 
key point in all experiments to stndy very peripheral ion-ion collisions is the 
organization of the primary or “Level Zero” (LO) trigger, which is needed in 
order to start the Data Acqnisition System (DAQ) of the experiment. 

The Higher Level (HL) triggers are nsed then fnrther only to stop or continne 
the data taking by the DAQ of the triggered event from the detectors. The 
main aim of the HL-triggers is to prevent storage of the physically nninterest- 
ing events. Dne to the hnite dead time of the DAQ the LO trigger rate cannot 
be too high, otherwise the data taking in the experiment will be fnlly blocked. 
Thns processes, which give a too high a load on the LO trigger detectors, 
should be considered as a most important source of background. 

The hrst trigger on coherent processes in relativistic heavy ion collisions was 
developed for the STAR experiment at the RHIC [331]. The LO-trigger is based 
on the charged particle multiplicity in the pseudorapidity interval \ri\ < 1, 
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measured in the Central Trigger Barrel (CTB) and Multi-Wire Chambers 
(MWC) at the endcaps of the TPC. The information from the trigger detectors 
is divided into 4 bins in azimuthal angle which can be used to apply cuts on 
event topology to select events with tracks lying in opposite 0-sectors. At the 
LI and L2 trigger levels the improved multiplicity information is available 
after rehning kinematics. At the L3-level the tracking information is available 
to apply momentum cuts and use the charge of the tracks. The decision time 
of the trigger levels 0, 1/2 and 3 were given by the authors as 1.7 ps, 100 ps 
and 10 ms respectively. Besides these triggers STAR also used a minimum bias 
trigger to select events when both nuclei were dissociated [50]. This trigger is 
based on zero-degree calorimeters (ZDC). The presence of a signal in each of 
the ZDCs corresponding to one or more neutrons with the full beam energy 
defines events when the nuclei are electromagnetically dissociated. 

The positioning of the main trigger detectors of an experiment is highly sen¬ 
sitive to the background processes. In this connection the low mass single and 
multiple e'''e“-pair production in heavy ion collisions has to be considered 
as one of the major backgronnd processes, because its cross section is really 
huge (200 kbarn for PbPb at LHC and 30 kbarn for AuAn at RHIC), the 
rapidity distribntion is rather flat (and thus covers all detectors) and all other 
signatures are quite similar to those of the photon-photon processes, because 
the e'''e“-pair production is also just a photon-photon process, see Sec. 7.1 
above. Electrons and positrons are produced preferably in the very forward 
or backward directions, see Fig. 25(b), but these pairs do not constitute a 
background for the ZDCs as their energy is rather small. Fig. 25(a), compared 
to the typical hadron energy in this pseudorapidity interval. But the central 
trigger detectors (i.e., the CTB in the case of the STAR experiment) can be 
essentially loaded by these pairs. Dne to their low transverse momentum 
an essential reduction of this backgronnd can be done, if one uses the central 
trigger detector at a snfficiently large radius, a strong solenoidal magnetic held, 
or both. (In the case of the CTB at STAR its radius, 2.2 m and B = 0.2 T, 
is sensitive to charged particle with Pi > 130 MeV [331], which reduces the 
background from soft e’''e“-pair signihcantly). On the other hand if the central 
trigger detectors are too large, we may expect some background from cosmic 
ray events. As was shown in [331] the rate of cosmic rays entering the TPC is 
due to the large size of the TPC (16 m^) almost 3 kHz, which is rather high. 
In [49] the cosmic ray hux through the STAR trigger detectors was estimated 
to give an important percentage of all background processes. 

Another possible trigger detector, the ZDC, can be contaminated by the elec¬ 
tromagnetic excitation of the ions alone to the Giant Dipole and other reso¬ 
nances, see Sec. 5.1. We are studying here the following two types of processes: 
one-photon exchange with the electromagnetic excitation of one of the ions and 
no particle detected in the central rapidity region and the production of one 
or more low mass e’''e“-pair together with the electromagnetic excitation of 
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one of the ions. The cross sections of these processes in heavy ion collisions at 
RHIC and LHC are epual to o' auAu^AuAu* 58b, o'AuAu^AuAu*{e+e~)ffyQ 13b, 
and opbPb^PbPb* = 114b, opbPb^PbPb*{e+e-)i,yo = 22b, respectively, using for 
the GDR excitation Eq. (79) and for the pair production Eq. (103). The 
similar process but with the mutual excitation of both ions have almost the 
same order of magnitude, as the probability for excite the GDR is not much 
smaller than one, see Sec. 5.1. For intermediate ions (CuCu at RHIC and ArAr 
at LHC) the exchange of additional soft photons leading to electromagnetic 
excitation is suppressed by a factor of about see Eq. (79). Therefore the 

background load on the ZDC as a potential trigger detector by electromag¬ 
netic excitation of the ions alone is important only in the case of heavy ion 
collisions, but it is negligible small for the intermediate and low mass ions. 


8.2.3 Trigger for Very Peripheral Processes 

In this subsection we briefly summarize the most essential aspects, which 
should be taken into account for a trigger to study the very peripheral ion 
collisions. 

The central trigger barrel (like the CTB in STAR) with endcups have to be 
considered as the major detectors for the LO-trigger to select very peripheral 
ion collisions. The trigger requirement is a low charged particle multiplicity 
rich in the central trigger detectors, let’s say 1 < rich < 20. The solenoidal 
magnet is an obligatory part in such type of experiments, because it prevents 
the overloading of the trigger barrel by background processes. Two backgro¬ 
und processes should be taken into account for an optimization of the trigger 
detector parameters: the processes of low mass multiple e+e“ production and 
the background from cosmic events. The basic parameters are the magnetic 
field in the solenoid B and the trigger barrel radius Rbar or the inner radius 
of endcups, Rinn- For the optimization of the trigger one increases the value 
of BR (where R = min{i?f,ar, Rinn}) to reduce the background from low mass 
multiple e’''e“-pair production up to a value, at which the background from 
cosmic ray events is still not too high. The optimal detector parameters corre¬ 
spond probably to a detector conhguration, when the background rates from 
the e'''e“-pair production and from cosmic events are comparable. The sig¬ 
nals from the ZDCs are not so important, due to their contamination with 
processes of the electromagnetic excitation of the heavy ions. But it could 
be useful for the suppression of Pomeron-Pomeron processes in the case of 
intermediate ion collisions. 

Another possibility for a very peripheral trigger in ion-ion collisions is the 
barrel Electromagnetic Calorimeter (ECAL). This is one of the most popular 
detectors used in experiments at e’''e“ and hadron colliders, see, for example, 
[286]. The requirement that the total deposited energy in the ECAL is suf- 
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ficiently low compared with the one in strong interaction of the ions, allows 
an effective selection of very peripheral events at a LO trigger level. Because 
the energy in the ECAL from the low mass multiple e’''e“-pairs is quite small, 
these events do not provide any background. The strong ion interactions in 
peripheral collisions as a background process can be effectively suppressed by 
using additional signals from the ZDCs, which for these events are very high 
as was discussed in Chap.5. Only cosmic ray events provide in this case some 
background, which can be suppressed at the HL trigger level, based on general 
event topology, see [50]. 


8.3 Detection of Some 'j'y-Processes 


Two-photon processes in heavy ion collisions have not yet been widely studied 
experimentally with the first heavy ion collider RHIC starting its operation 
only a year ago. However, a group of experimentalists at the STAR detector 
have studied the possibility to observe 77 -processes during the last six years. 
Earlier 77 -physics has been studied at e’''e“-colliders by almost all collabora¬ 
tions. The experience gained of these e’''e“-experiments to detect 77 processes 
can be applied now to the heavy ions experiments. Below we give an overview 
of the major experimental observations of 77 -processes at e’''e“-colhders and 
the studies done for heavy ion experiments. 


8.3.1 Experiments at e’''e -Colliders 

All e’''e“-collider experiments which have been running so far, have accumu¬ 
lated a large experience for the detection of 77 processes. Although the signa¬ 
tures of two-photon processes in electron-position collisions differ from those 
in heavy-ion collisions, they still have many common features which allows to 
take over their experience to the future heavy ion collider experiments. Many 
recent experimental data from e’''e“-colliders (LEP-II, CLEO, KERB) are for 
quasi-real photon interactions, like those in heavy-ion collisions. Here we give 
a review of the 77 processes detected by those experiments. 


77 —> hadrons. This reaction was studied by the L3 collaboration at LEP 
[ 101 ] to measure the total cross section for the scattering of two real photons 
into hadrons. The two-photon cross section < 7(77 —hadrons) was measured in 
the invariant mass range of 5 < Hdy.y < 185 GeV. Hadronic two-photon events 
were selected by the following criteria: 

• Events with the scattered beam electrons detected were rejected. This re¬ 
stricts the virtuality of the interacting photons to < 8 GeV^ with an 
average value of ~ 1.5 x 10“^ GeV^. 
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• The total energy in the electromagnetic calorimeter was required to be be¬ 
tween 0.5 and 50 GeV to suppress beam backgrounds and exclude radiative 
events. 

• The detected multiplicity should be greater than six to exclude events con¬ 
taining r-pairs. 

After the selection the visible effective mass of the event, ITvis was measured. 
An unfolding procedure based on the Monte Carlo event generators PHOJET 
and PYTHIA was used to calculate the two-photon mass The system¬ 
atic uncertainties due to the Monte Carlo generators were estimated to be 
negligible for < 65 GeV, but are important for higher 77 -masses. 


Production of light resonances. The L3 collaboration studied also light res¬ 
onance production in 77 -colhsions [332]. The resonances were studied in the 
and final states. The 7(1440), /i(1420), /i(1285) were ob¬ 

served. The mass region between 1200 and 1500 MeV is expected to contain 
several states. The pseudoscalar state 7(1440) was observed in hadron colli¬ 
sions and in radiative J/T decay, but not in 77 -colhsions. Only a very low 
upper limit of its two-photon width was measured, which supports the in¬ 
terpretation of 7 ( 1440 ) as a glueball candidate. Axial-vector states /i(1420), 
/i(1285) are also present in this mass region. 

Since these are exclusive processes, their selection criteria differ from those for 
the process 77 —X above: 

• The TT^TT" was identified by the secondary vertex separated from the 

interaction point, the effective invariant mass of the two pions and the angle 
correlation between the momentum and the two hnal tracks. 

was identihed by the two tracks coming from the interaction point 
and their dE/dx measurement being consistent with a Xtt hypothesis. 

• The final state 77 r+ 7 r“ was identified by the presence of two photons and 
two tracks of opposite charge and the 77 invariant mass. 


77 —> K^K~, 77 — KgKg. These reactions were studied by the L3 collab¬ 
oration at LEP [333], by the BELLE collaboration at KEKB [334] and by the 
CLEO at CESR [335]. The events were selected according to usual criteria 
for exclusive reactions, i.e., the requirement of two or four charged tracks for 
the K^K~ and hnal states respectively, with a net charge of zero, 

mesons were identihed by secondary vertex reconstruction and the invariant 
mass of the secondary track pairs. The vector sum of the transverse momenta 
of all tracks should be small. These experiments aimed to study the resonances 
in the mass region from 1.3 to 2.3 GeV, where glueballs are expected to be 
found. A tensor state /j(2220) was observed in the radiative decays of J/T 
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and was considered a good candidate for a glueball. Since gluons do not cou¬ 
ple to photons they cannot be produced directly in 77 collisions. The 
experiments did not observe any resonance around the mass 2220 MeV, thus 
they support the glueball interpretation of the /j( 2220 ) state. 


77 — ccX, 77 —> bbX. Open charm and beauty production was measured 
by the L3 collaboration [336,337]. The event selection was performed in two 
steps. The hrst one selects hadronic states and the second one identihes c 
and b quarks by their semileptonic decays. The hadronic hnal states were 
selected by requiring at least 5 tracks and a visible energy below y/s/3. The 
visible mass of the events was restricted to ITvis > 3 GeV. Then the unfolding 
procedure based on either PYTHIA or PHOJET was applied to obtain the 77 
c.m. energy. 


77 — D*^X. The inclusive production of i7*-mesons was measured by L3 
[338] and OPAL [339]. D*^ mesons were detected via their decay chains 
D*^ —and > K~tt^, . The reconstruction of these 

decay chains requires a sample of events containing hadronic hnal states. The 
hadronic events were selected by cuts on the energy measured in the electro¬ 
magnetic and hadronic calorimeters and using tracing information. 


11 ^ Vci XcOi Xc2- Charmomium production was studied by CLEO [340,341] 
and L3 [342]. The 77 and total widths of these states were measured. CLEO 
searched for rjc meson in the decay mode, and Xco and Xc 2 via 

their decay into 7 r+ 7 r“ 7 r+ 7 r“. The L3 collaboration detected rjc in various de¬ 
cay modes with tt^, K^, K^, tj' and in the hnal state. Therefore the 

selection criteria were the usual ones for exclusive events with such type of 
particles. 

Summarizing the experience collected by e’''e“ experiments to study 77 in¬ 
teractions, one can see typical event selection criteria, which can be adapted 
to heavy ion collisions. The 77 events are characterized by the absence of 
the scattered beam particles, a low (or hxed in the case of exclusive reactions) 
multiplicity in the hnal state, low transverse momentum of the hnal state, zero 
net charge of the tracks, low (compared to the annihilation processes) energy 
deposited in the detectors. Depending on the reaction to be studied, diherent 
topology cuts can be applied. In the case of inclusive reactions a Monte Carlo 
based unfolding procedure is used to obtain the 77 mass from the visible event 
mass. 
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8.3.2 Experiments at Heavy Ion Colliders 


The STAR collaboration is the first one to study 77 -processes in heavy ion col¬ 
lisions in the AuAu collisions at 100-250 A GeV at RHIC [47,46,343,331]. The 
scope of their interest is lepton pair production, meson pair production, meson 
spectroscopy, baryon pair production, and possibly charmonium production. 
They studied background processes, trigger and experimental techniques to 
select 77 -processes. 


Background. STAR studied several major sources of background for the two- 
photon processes: peripheral hadronic nuclear collisions, beam-gas interac¬ 
tions, cosmic rays, photon-nucleus collisions and Pomeron processes. Cosmic 
rays and beam gas reactions are a problem mainly at the trigger stage, be¬ 
cause the more precise data analysis at the later stage rejects such kind of 
events via track and vertex information. Detailed Monte Carlo simulation has 
been performed to separate signal from background [343]. Peripheral hadronic 
nuclear collisions were simulated with two standard nuclear event generators, 
FRITIOF [344] and VENUS [345] with impact parameters from 12 to 20 fm. 
Beam-gas interactions were simulated by the same event generators as Aup 
and AuN collisions where hydrogen p and nitrogen N represent the beam-gas 
content. Photonuclear reactions were studied by DTUNUC [346]. Only yAu 
collisions with photon energies larger than 10 GeV were simulated. Cosmic 
rays were simulated with a flux of 1.8 x 10 ^ m“^s“^ and the trigger rates in 
the TPC were estimated. The most important cuts, which discriminate 77 
events from background, were found to be: 

• Multiplicity: many two-photon reactions, which can be detected at RHIC, 
have two or four charged particles in the hnal state; 

• the sum of the transverse momentum of the hnal state particles should be 

small, of the order \^hc/R; in the analysis a cut < 40 MeV was used; 

• the c.m. rapidity distribution of the 77 system is centered around zero 
with a narrow width; the cut < 1.0 was applied; this cut is realized 
automatically by the detector acceptance and reduces also the beam-gas and 
photonuclear interactions since these are characterized by an asymmetric 
particle emission. 

Monte Carlo simulation performed by the STAR collaboration [343] showed 
that these cuts are able to separate the signals from background. 


First experimental results from STAR. RHIC had a hrst AuAu run in 2000, 
and the hrst reports on observation of coherent interactions in nuclear colli¬ 
sions appeared at the beginning of 2001. In [50] the hrst observation of exclu¬ 
sive production is presented. The are produced with small perpendicular 
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momentum, as expected, if they couple coherently to both nuclei. Here we pay 
special attention to the event selection used. 


The production of p°s were studied in the decay channel p° —>• . The data 

on AuAu collisions at y/SNN = 130 GeV was collected during the Summer 2000 
run. The p production was studied with two separate triggers. The topology 
trigger was designed to trigger p° decays detected in the central trigger barrel 
(CTB). The topology trigger divided the CTB into 4 azimuthal quadrants as 
shown in Fig. 42. This trigger selects events with at least one hit in the North 


Transverse Momentum 




Fig. 42. Topology trigger on coherent 
interactions at STAR. Figure taken 
from [50]. 


Fig. 43. The transvers momentum P±_ of the 
two pions triggered with the topological trig¬ 
ger. The enhancement at small transverse 
momentum, coming from the very periph¬ 
eral collisions, is clearly seen. Figure taken 
from [50]. 


and the South sectors. The top and bottom quadrants were used to reject most 
cosmic rays. The selected events were then passed through a level-3 trigger, 
which rejects events with more that 15 tracks or a vertex far outside the 
interaction point. These cuts rejected 95% of the events and rejected events 
from central collisions, beam-gas interactions and cosmic rays. 


A minimum bias trigger used both zero degree calorimeters (ZDC) to select 
events, where both nuclei dissociated. This trigger required a coincidence be¬ 
tween the two ZDCs with a single (or more) neutron depositing their energy 
in each ZDC. 


In the analysis STAR selected events with exactly two tracks with a vertex 
within 2 cm of the center of the TPC in the transverse plane, and within 2 m 
along the beams. A slightly acoplanarity of the track pairs was required to 
reject the remaining cosmic rays. The sum of the transverse momentum of 
the track pairs has a peak at low pj--, which shows that they are coming from 
coherent collisions, see Fig. 43. The invariant mass spectra of the 7r'''7r“-pairs 
in events selected by both kind of triggers shows a peak around the mass. 
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The STAR collaboration was the first one who studied experimentally the 
coherent interaction of heavy ions. In 2001, STAR has upgraded the detectors 
and the trigger, which will allow to increase the data taken by several orders 
of magnitude. Also it allows to study other reactions in coherent heavy ion 
interactions. 

8.4 Perspectives of ion experiments 


As was discussed before, 77 processes have been studied thoroughly at 
experiments up to the 77 -invariant masses of about 3 GeV. Recently the 
R-factory at KEK, the e’''e“-collider KERB started operating, and a few pa¬ 
pers have already appeared concerning the detection of 77 -processes in the 
experiment BELLE at KERB. These recent reports have attracted some at¬ 
tention of the physics community towards RERB as a new facility to study 
two-photon physics. Actually RERB provides an '^unprecedented luminosit]/' 
[347] of 4.49 X 10^^ cm“^s“^. Thus the 77 -luminosity achieved at RERB is 
several orders of magnitude higher than that at RHIC in the whole range of 
77 invariant masses possible in both colliders. So far the BELLE collaboration 
has reported only on a few measurements of 77 processes [334,348], based on 
the hrst run stored in 2000 , but one expects that a lot of high-statistics results 
will appear from BELLE soon and the low-mass 77 -physics will be exhausted 
there on the level of femtobarn cross sections. However the physics domain 
of photon-Pomeron and Pomeron-Pomeron processes at RHIC is out of the 
competition with any lepton collider. 

For the experiments at LHC, the possibilities to study very peripheral ion 
collisions are currently discussed by the CMS and ALICE collaborations. The 
effective 77 -luminosity at LHC is large enough due to the high collision energy. 
In addition the 77 mass range, that can be achieved in ion collisions there, 
allows to study 77 processes well beyond the mass range studied at LEP, see 
Fig. 6 . 

CMS considers an extensive heavy ion program, which includes also the physics 
of very peripheral ion collisions [22,349]. This program is quite attractive, 
especially as CMS provides a detector with a rather wide aperture and allowing 
thus to study various 77 -reactions. The physics of very peripheral ion collisions 
at ALICE has been discussed in a number of papers, see [350,55,351]. Although 
ALICE is a quite promising setup for studies of this kind of physics a trigger 
conhguration for very peripheral events has not been developed up to now. 
Nevertheless some possible options were proposed, see, i.e., [352], and problem 
is still under discussion. 
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9 Conclusion 


In this report we describe the basic properties of very peripheral ion-ion col¬ 
lisions. Due to the very strong electromagnetic helds of short duration, new 
possibilities for interesting physics arise. The study of these electromagnetic 
processes, that is, photon-photon and photon-ion collisions is an interesting 
option, complementing the program for central collisions. It is the study of 
events characterized by relatively small multiplicities and a small backgro¬ 
und (especially when compared with the central collisions). These are good 
conditions to search for new physics. The method of equivalent photons is 
a well established tool to describe these kinds of reactions. Reliable results 
of quasi-real photon fluxes and 77 -luminosities are available. Unlike electrons 
and positrons heavy ions and protons are particles with an internal structure. 
Effects arising from this structure are well under control and minor uncer¬ 
tainties coming from the exclusion of central collisions and triggering can be 
eliminated by using a luminosity monitor from jj,- or e-pairs. 

The high photon fluxes open up possibilities for photon-photon as well as 
photon-nucleus interaction studies up to energies hitherto unexplored at the 
forthcoming colliders RHIC and LHC. Heavy ion colliders are a unique tool 
to study (quasi-real) photon-nucleus collisions in a hitherto inaccessible en¬ 
ergy range. First experimental results at RHIC on coherent p production on 
An nuclei are forthcoming, as well as, pion- and electron-pair production in 
very peripheral collisions. A wealth of new data is to be expected. Interesting 
physics can be explored at the high invariant 77 -masses, where detecting new 
particles could be within range. Also very interesting studies within the stan¬ 
dard model, i.e., mainly QCD studies will be possible. This ranges from the 
study of the total 77 -cross section into hadronic hnal states up to invariant 
masses of about 100 GeV to the spectroscopy of light and heavy mesons. The 
production via photon-photon fusion complements the production from single 
timelike photons in e+-e“ collider and also in hadronic collisions via other 
partonic processes. 

A good trigger for very peripheral collisions is essential in order to select 
photon-photon events. As it was shown in this report, such a trigger will be 
possible based on the survival of the nuclei after the collision and the use 
of the small transverse momenta of the produced system. A problem, which 
is difficult to judge quantitatively at the moment, is the influence of strong 
interactions in grazing collisions, i.e., effects arising from the nuclear strato¬ 
sphere and Pomeron interactions. With the pioneering experiments of the (Ul¬ 
tra) Peripheral Collisions Group at STAR at RHIC this held has dehnitely left 
the area of theoretical speculations and entered the area to be experimental 
feasible. 
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Very peripheral collisions with Photon-Pomeron and Pomeron-Pomeron in¬ 
teractions, that is, diffractive processes are an additional interesting option. 
They use essentially the same triggering conditions and therefore one should 
be able to record them at the same time as photon-photon events. 

This review is written in a time of rapid progress, both from an experimental 
and theoretical point of view. A few theorists are studying various aspects of 
very peripheral collisions, many study the held of 7 -hadron processes with, 
e.g., the interactions of vector mesons in the nuclear medium and the vast 
held of 77 collisions ranging from QCD studies to possible new physics. On 
the experimental side, experiments on very peripheral collisions are presently 
performed and analyzed at RHIC. With these experiments the held has left 
the area of purely theoretical ideas and entered into the reality of experiments. 
The hrst heavy ion beams are expected at LHC in 2007 and experiments are 
in a construction stage. Certainly, the experience at RHIC will help to plan 
(even) better for LHC. The future is coming and one thing is clear: with all 
the new forthcoming experimental results the next review of the held of very 
peripheral collisions will be very diherent from this one. 
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